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ABSTRACT 



A simplified pol} energetic cell theory is formulated 
to determine spatially averaged energy dependent thermal 
fluxes in the moderator,, cladding, and fuel regions within 
the unit cell of a reactor lattice. The derived spectra 
are then utilized in the calculations of the thermal integral 
parameters and average cross sections required for reactor 
computations . 

The cell theory, as formulated, postulates an infinite 
moderator region with the absorption cross section of this 
region appropriately modified to account for the neutron 
leakage into and absorption by the fuel element. The mod- 
ifications to the moderator absorption cross section are 
formulated both in terms of the net current at the fuel 
element -moderator interface and in terms of energy dependent 
moderator and fuel element escape probabilities, the latter 
approach offering physical transparency and ease of cal- 
culation . 

Analytic expressions for the escape probabilities are 
presented, integral transport theory being applied to the 
fuel element region, while diffusion theory is utilized in 
the moderator region. Using these analytic expressions, 
the theory is applied to actual lattices in the form of the 
light water moderated and uranium dioxide fueled cores of 
the Cornell University Zero Power Reactor. Room temperature 
parameters and their temperature coefficients are determined 
using both the monatomic gas model and the Nelkin water 
kernel to describe the energy transfer process in the moder- 
ator. Calculations are made with PROGRAM COUTH, a Fortran-63 
program written for use with the Control Data Corporation 
1604 digital computer. A typical lattice calculation 
including the computation of the spatially averaged fuel, 
cladding, and moderator spectra and the thermal integral 
properties and average cross sections takes approximately 
thirty-five seconds of computer time. This figure is exclusive 
of the compiling time and the time required to calculate the 
moderator scattering kernel. 

In an attempt to estimate the accuracy of the calculational 
results, the method is applied to the Brookhaven National 
Laboratory uranium dioxide cores and the results are then com- 
pared with those predicted by Honeck's THERMOS code. Dis- 
advantage factors agree to within 1.0% while the thermal 
utilizations agree to within 0.5%. A study of the sensitivity 
of the calculated integral parameters to variations in the 
input data leads to the assignment of rather small uncertainties 
in the results calculated with the simplified cell theory. 
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CHAPTER 1 



INTRODUCTION 



The analysis and design of any thermal reactor requires 
an accurate knowledge of the various and competing neutron 
reaction rates within the reactor assembly. This, in turn, 
implies a knowledge of the distribution of thermal neutrons 
in space and energy. In the calculation of the thermal neu- 
tron spectra one must consider the two processes which affect 
these distributions in a heterogeneous reactor lattice. The 
first of these is the moderation process wherein a neutron 
is scattered by a moderator atom or molecule with a subsequent 
change in energy. The second process is the transport or 
diffusion mechanism whereby the neutrons are carried into 
regions of relatively different absorption characteristics . 
Although the time independent Boltzmann's equation or the 
integral transport equation— ^ correctly describes both the 
moderation and transport phenomena, the difficulty in 
obtaining solutions to these equations have led previous 
investigators to develop their work along one of two lines. 

Initial attempts to characterize the distribution of thermal 

2 3/ 

neutrons were done for infinite homogeneous media—* — where 



— A. M. Weinberg and E. P. Wigner, The Physical Theory of 
Neutron Chain Reactors , pp. 221-234, The University 
of Chicago Press , Chicago, 1958. 

— ^ M. J. Poole, M. S. Nelkin, and R. S. Stone, "The Measure- 
ment and Theory of Reactor Spectra, " Progress in Nuclear 
Energy , Series I, Vol. 2, pp . 91-164, Pergamon Press, 

New York, London, Paris, and Los Angeles, 1958. 
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there is no net transport or diffusion and the spatial 
variable is of no consequence . This approach allows con- 
centration on the moderation or energy transfer process 
and sophisticated mathematical models may be incorporated 
into the calculation. 

The other approach, as exemplified by multigroup 
diffusion calculations attempts to consider spatial 
effects at the expense of the details of the moderation 
process. In these calculations, the scattering kernel 
is often characterized by a simple model which allows 
ease of mathematical description but which may lack 
accuracy in its approximation of the energy transfer 
process . 

More recently, interest has been focused on the 
use of transport theory— ^ for dealing with the heterogeneities 
present in a reactor lattice. This approach is embodied 
in the work of Honec k--* — ^ where both the spatial and 



3 / 

— Neutron Spectra in Lattices and Infinite Media , 

Proceedings of the Brookhaven Conference on Neutron 
Thermalization, BNL 719, Vol. II, Brookhaven National 
Laboratory (1962) . 

4 / 

— A. Hassitt, "Methods of Calculation for Heterogeneous 

Reactors," Progress in Nuclear Energy , Series I, Vol. 2, 
pp . 271-313, Pergamon Press, New York, London, Paris, 
and Los Angeles, 1958. 

— ^ Neutron Spectra in Lattices and Infinite Media, see 
Footnote 3 . 

— ! H. C. Honeck, "THERMOS -- A Thermalization Transport 

Theory Code for Reactor Lattice Calculations" BNL 5826, 
Brookhaven National Laboratory (1961) . 
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moderation effects on the spectra are considered through a 

space and energy calculation utilizing integral transport 

theory with only slight restrictions on the complexity of 

the scattering model used. 

Still another approach is that found in the work of 
8 9/ 

Leslie.— 2 — In this formulation the effect of moderation 

and heterogeneity are both considered, but the spatial 
effects are mathematically subordinated by consideration 
of only the spatially averaged fluxes in the moderator and 
the fuel. This method of attack leads to physical trans- 
parency while allowing one to maintain and concentrate on 
the importance of the moderator scattering model. 

This work is an extension of Leslie's ideas, culmina- 
ting in a computer code which provides for the rapid cal- 
culation of the spatially averaged thermal neutron spectra 
in the moderator, fuel, and cladding regions of a reactor 
lattice. The simplified cell theory utilized in these cal- 
culations is formulated in Chap. 2. A general formulation 
is presented as well as a formulation in terms of escape 
probabilities defined for the fuel element and the moder- 
ator regions . The escape probabilities utilized are a 
polyenergetic generalization of those defined by Amouyal, 

— ! H. C. Honeck, "The Calculation of the Thermal Utilization 
and Disadvantage Factor in Uranium/Water Lattices," 

Nucl . Sci. Eng. 18, 49-68 (1964). 

— ^ D. C. Leslie, "The Spectrox Method for Thermal Spectra 
in Lattice Cells," AEEW-M211, United Kingdom Atomic 
Energy Authority (1962) . 
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Benoist, and Horowitz— ^ in their one velocity thermal 
utilization theory. Chapter 3 deals with the escape prob- 
abilities defined for the moderator region,, while Chap. 4 
details the calculational methods available for the deter- 
mination of the fuel element escape probability. Mathe- 
matical models for the energy transfer or scattering kernel 
are discussed in Chap. 5. Emphasis is placed on the models 

actually used for calculational purposes; the monatomic gas 

1110 / 11 / 
kernel — 2 and the Nelkin water kernel. — ' Chapter 6 

defines the calculations necessary in the determination of 

spatially averaged, energy dependent thermal fluxes for the 

moderator, cladding, and fuel regions in a lattice utilizing 

the simplified cell theory. Integral thermal properties are 

defined as are the averaged cross sections necessary for the 

accurate calculation of reaction rates. In Chap. 7 the 

formulated theory is applied to actual lattices . A study 

is made of the low enrichment, high density, uranium dioxide 

fueled and light water moderated Cornell University Zero 

Q / 

— D. C. Leslie, "The Calculation of Thermal Spectra in 
Lattices," Proceedings of the Brookhaven Conference 
on Neutron Thermalization, BNL 719, Vol. II, pp . 592- 
609, Brookhaven National Laboratory (1962) . 

— ^ A. Amouyal, P. Benoist, and J. Horowitz, "Nouvelle 
Mdthode de Determination du Facteur d 'Utilization 
Thermique d'une Cellule," J. Nucl. Energy 79 (1957). 

— ^ E. P. Wigner and J. E. Wilkins, Jr., "Effect of the 

Temperature of the Moderator on the Velocity Distri- 
bution of Neutrons with Numerical Calculations for H 
as the Moderator," AECD-2275 (1944). 
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14 / 

Power Reactor cores . — Temperature effects are predicted 
through the use of temperature dependent models for the 
scattering kernel. The accuracy of these results are dis- 
cussed in Chap. 8. This chapter also includes a comparison 
of the results predicted with the simplified cell theory to 
those calculated by Honeck with THERMOS— ^ for the Brookhaven 
light water moderated, and uranium dioxide fueled lattices . 
Chapter 9 contains a summary of the work undertaken as well 
as the conclusions which may be drawn from it. 



12 / 

— ' J. E. Wilkins, Jr., "Effect of Temperature of the Moder- 
ator on the Velocity Distribution of Neutrons for a 
Heavy Moderator," CP-2481 (1944). 

— ^ M. S. Nelkin, "Scattering of Slow Neutrons by Water," 
Phys . Rev. 119 , 741-746 (1960). 

— ^ D. D. Clark, "The Cornell University Nuclear Reactor 
Laboratory," CURL-1, Cornell University (1961). 

— ^ H. C. Honeck, see Footnote 7. 



CHAPTER 2 



A SIMPLIFIED CELL THEORY 
2.1 General Formulation 

Recognizing that the neglect or oversimplification of 
the spatial effects in a heterogeneous lattice leads to 
answers which have little relationship to the physical 
processes taking place, one is prompted to overemphasize 
the spatial variations by resorting to multigroup diffusion— 
or transport—^ codes. However, these codes consume a con- 
siderable amount of computer time, particularly if the more 
complex models of the scattering kernel are used. 

Thus, there is a need for a simple polyenergetic 
formulation for determining the thermal distribution of 
neutrons in the various regions of a heterogeneous lattice. 
Hopefully one could retain the calculational ease of the 
homogeneous, infinite medium thermal spectrum computation, 
thereby offering only slight limitations on the complexity 
of the scattering model used, while still accounting for 
the heterogeneities of the lattice in a meaningful way. 

This section formulates a simplified cell theory for 
use in the calculation of thermal spectra and parameters 
in a lattice composed of uniformly arrayed unit cells . The 

— / R. S. Varga, "Numerical Methods for Solving Multi- 
Dimensional Multigroup Diffusion Equations," 

Proceedings of Symposia in Applied Mathematics , 

Vol. XI, pp .~T64-lo9, NucTear Reactor Theory, 

American Mathematical Society, Providence (1961). 
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heterogeneous effects due to neutron transport within the 
cell are formulated in a physically transparent manner, 
thereby allowing sufficient emphasis to be placed on the 
moderation process in any calculation. 

Consider first the neutron balance equation to be 
solved in the case of an infinite, homogeneous, moderating 
medium:—^ 



2®(E) + Z“(E) 



cp m (E) 



S(E) + 



c 



2(E HE)cp m (E ' )dE ' . 



( 2 . 1 ) 



E c is the thermal cutoff energy above which the moderator 

flux per unit energy, cp m (E), is assumed to be 1/E, while 

2(E'-*-E) is the energy transfer or scattering kernel 

expressed in (cm-eV)""*". £ m (E) and £ m (E) are the macroscopic 

absorption and scattering cross sections respectively of 

the moderator, X m (E) being defined by 

s 



2g(E) = 



2(E — 3-E ' )dE 1 



( 2 . 2 ) 



o 



17 / 

— H. C. Honeck, "THERMOS -- A Thermalization Transport 
Theory Code for Reactor Lattice Calculations," 

BNL-5826, Brookhaven National Laboratory (1961) . 

— / M. J. Poole, M. S. Nelkin, and R. S. Stone, "The Measure- 
ment and Theory of Reactor Spectra," Progress in 
Nuclear Energy , Series I, Vol. 2, pp . 91-164, Pergamon 
Press, New York, London, Paris, and Los Angeles, 1958. 
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The source of thermal neutrons, S(E), due to down scattering 

from energies above E c is given by 

oo 



S(E) = 



2(E'— *~E)cp™(E' )dE' 



(2.3) 



where cp™(E) is proportional to 1/E. 

Neutron conservation requires that 



S (E)dE = 



2?(E) cp m (E)dE 



(2.4) 



which is obtained by integrating Eq . (2.1) over energy from 



0 to E . 
c 

The problem as stated places no restrictions on the form 
of the scattering kernel, 2(E' — »-E), however, we require that 
it satisfy detailed balance: 



2(E ' — *-E)M(E ' ) = 2(E — ^E ' )M(E) , (2.5) 

where the Maxwellian distribution, M(E), is given by 

M(E) = E/(kT) 2 exp(-E/kT) . (2.6) 

The moderator temperature in electron volts is equal to kT, 

T being the moderator physical temperature while k is 
Boltzmann's constant. 

There are of course practical limitations which restrict 
the complexity of the kernel and the number of energy points 
to be utilized in the computer solution of Eq . (2.1). 
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The infinite medium calculation is now extended to the 

heterogeneous lattice utilizing the ideas embodied in Leslie's 

19 / 

SPECTROX method. — ' Since the only effect which the fuel 
element has upon the spectrum in the moderator of a lattice 
cell is to modify the transport process in the vicinity of 
the fuel element, while contributing very little to the modera- 
tion process itself; one may postulate an infinite moderator 
region with a contrived energy dependent absorption cross 
section which includes the heterogeneous properties of the 
lattice. In other words, if the absorption cross section of 
the moderator can be appropriately modified to account for the 
neutron leakage into and absorption by the fuel element, one 
can then turn his attention to the choice of a proper 
scattering kernel. The calculational ease of Eq . (2.1) is 
still available along with the fact that the non moderating 
part of the reactor lattice has been accounted for. 

Considering the energy dependent absorption rates, one 
defines 



absorption rate in the fuel element 

f /r.s _ at energy E, per eV 

^ ' absorption rate in the lattice cell 
at energy E, per eV 



(2.7) 



The "fuel element" includes the fuel rod and its cladding, 
and in the case where there is no neutron absorption in the 



cladding, f(E) can be considered to be the energy dependent 
thermal utilization. In terms of f(E) the neutron 



— / D. C. Leslie, "The SPECTROX Method for Thermal Spectra 
in Lattice Cells," AEEW-M211, United Kingdom Atomic 
Energy Authority (1962) . 
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loss rate per unit volume of the moderator due to leakage 
into and absorption by the fuel element may be expressed as 
2™(E)V n (E)f (E)/ [1-f (E) ] , where cp^CE) is the spatially 

a 

averaged flux in the moderator region having a volume of 



9(E) 




cp m (r,E)dV 



( 2 . 8 ) 



Equation (2.1) may now be rewritten for the cell in a 
heterogeneous lattice as 



2g(E) + Z“ l (E) 









m 



^(E) = S (E) + 



2(E'-^E)9 m (E' )dE' . 



o 



(2.9) 



S^E) is now the spatially averaged slowing source of neutrons 
which appear in dE at E due to scattering collisions at 
energies greater than E c . 



S(E) 




S(r,E)dV . 



( 2 . 10 ) 



As it stands, Eq . (2.9) is correct for an infinite 
reactor lattice and it could be modified to include the 
effects of thermal leakage out of the reactor assembly. 
However, since one is usually interested in the intensive 
properties of the lattice, conventional approaches neglect 
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this thermal leakage calculation. In addition, Eq . (2.9) 
becomes approximate as soon as one postulates appropriate 
mathematical models for the scattering kernel, Z(E' — ^E), 



in the quest for a practical and economic computer solution. 

The most general formulation for the relative fuel 
element absorption may be obtained by considering a unit 
cell in an infinite reactor lattice. Assume a source of 
neutrons in dE at energy E and position r, N(r,E), so that 
the spatially average source of neutrons, N(E), per unit 
volume, V, is expressed by: 



Given N(r,E), one can in principle calculate the net current 

into the fuel element, J n (E), at each energy. The total 

number of neutrons absorbed by the fuel element at energy 

E is then AJ (E), where A is the surface area of the fuel 
n 

element. Since the total number of neutrons produced at E 
must equal the total number absorbed, the number of neutrons 
absorbed by the moderator is VN(E) -AJ n (E) . The relative 
fuel element absorption is then expressed as : 



and for the relative fuel element absorption, f (E) / [1-f (E) ] , 




( 2 . 11 ) 




( 2 . 12 ) 



In the assumption of a spatially dependent neutron 
source, N(r,E), and the subsequent calculation of the net 
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current into the fuel element , any of several means may be 
utilized including the use of multigroup transport codes , 
integral transport theory, or energy dependent diffusion 
theory. The accuracy of the results is directly related 
to the effort expended on these calculations, and this fact 
must be reconciled with the desire to produce economy in 
the use of the computer along with the knowledge that 
certain errors will be introduced with the use of mathe- 
matical models for the scattering kernel. 

It is to be noted that by incorporating and concen- 
trating all of the heterogeneous effects into the function 
f(E), the calculational simplicity of the infinite, 
homogeneous case has been retained as is seen in Eq . (2.9). 
A discussion of practical approximations for f(E) is found 
in Secs. (2.2) and (2.3). 

2.2 A Practical Approximation of Heterogeneous Effects 

In formulating an analytic approximation for f(E), the 
function which takes into account all of the heterogeneous 
effects in the lattice cell, one is forced to compromise 
between accuracy and computational ease. As an initial 
approach, the following simple model is postulated in an 
attempt to investigate the feasibility of this method in 
determining the thermal spectra within a lattice cell . In 
this work we have restricted ourselves to systems with 
cylindrical symmetry, recognizing that analogous results 
for spheres and slabs are obtainable. 
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Consider a cylindrical unit cell situated in the 

regular, rectangular or hexagonal, array of a lattice. 

It is characterized by the fuel element radius, r , and 

the moderator radius, r^ . In determining r^ from the actual 

20 / 

lattice parameters, the Wigner-Seitz criterion — is applied. 

A spatially flat source of neutrons in the moderator, N(E), 

is assumed. The validity of this assumption for the 

calculation of integral properties has been discussed by 
21 / 22 / 

Pomraning . — Leslie — has indicated that multigroup 
transport calculations yield moderator flux shapes which 
compare favorably with those predicted by one group 
diffusion theory. In light of this it is proposed to 
calculate the net current into the fuel element, J (E), 
through the extension of monoenergetic diffusion theory to 
the polyenergetic case. The flux shape in the cell 
moderator at any energy E is assumed to be spatially dis- 
tributed in a manner given by solution of the constant 
cross section diffusion theory approximation, using cross 
sections evaluated at the energy E: 

D m (E)V 2 cp in (r,E)-2“(E)cp in (r,E) + N(E) = 0 . (2.13) 



20 / 

— S. Glasstone and M. C. Edlund, The Elements of Nuclear 

Reactor Theory , p. 265, D. Van Nostrand Company, Inc., 
Princeton, Toronto, London, and New York, 1952. 

oi / 

— G. C. Pomraning, "On the Validity of the Constant Source 
Assumption for the Cell Problem," Nucl. Sci . Eng. 18 , 
400-403 (1964) . 

22 / 

— D. C. Leslie, see Footnote 19. 
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The boundary conditions applied are 



Vcp m (r 1 ,E) = 0 , 



(2.14) 



and 




(2.15) 



D m (E) = l/32™ r (E) is the moderator diffusion coefficient , 



and d(E) is the energy dependent linear extrapolation length. 
The linear extrapolation length may be determined from d Q (E), 
the linear extrapolation length into a black cylindrical 



probability, P e (E) . This has the added advantage of 
eliminating the need to state anything explicit about the 
shape of the flux within the fuel element . The fuel element 
spatial flux is implied in any calculation of the escape 
probability, but by dealing strictly with P fi (E) one can 
extract and utilize in a simple manner the accuracy of 

r\ / j 

transport theory calculations . — ' 

Considering the escape probability of the fuel element 

25 / 2 

as its albedo, — and including the term in V cp which can 

— ^ A. M. Weinberg and E. P. Wigner, The Physical Theory of 
Neutron Cha in Reactors , p. 765, The University of 
Chicago Press, Chicago, 1958. 

— ! K. B. Cady and M. Clark, Jr., "Neutron Transport in 

Cylindrical Rods," Nucl. Sci. Eng., lj$, 491-507 (1964). 

?S / 

— ' S. Glasstone and M. C. Edlund, p. 129, see Footnote 20. 



^(E) is the moderator macroscopic absorption cross section 

a 
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be obtained in the derivation of Fick's Law 



26/ 



one can write 



cp m (r o ,E)-2D m (E)Vcp m (r o ,E) + G(E) V 2 cp m (r o , E) 
cp m (r o ,E)+2D m (E)Vcp m (r o ,E) + G(E) V 2 <p m (r o , E) 



(2.16) 



For a black cylindrical hole, P g (E) is by definition zero so 
that G(E) is given by 



Vcp m (r E) 

G(E) = - y - -S 

V 2 cp m (r E) 



2D (E) -d (E) 
m x ' o 7 



(2.17) 



In the actual cell, where P e (E) is known from transport theory 
calculations, substitution of Eq . (2.17) into Eq . (2.16) 
yields with rearrangement 



<p m (r Q ,E) 
Vcp m (r o ,E) 



d(E) = 



4-D (E)P (E) 



m 



l-P e (E) 



+ dfE) . (2.18) 



An approximate expression for d Q (E) is given in Appendix C 
while the calculation of the fuel element escape probability 
is discussed in Chap. 4. 

The relative fuel absorption, f (E)/[ 1-f (E) ], may now be 
found from the solution of Eq . (2.13) utilizing the boundary 
conditions in Eqs . (2.14), (2.15), and (2.18). For a 
cylindrical cell, Eq . (2.12) may be expressed as 



l-f(E) 

HTE r 



N(E)r (y 2 -l) 

2 i , 

2D m (E)Vcp m (r 0 ,E) 



(2.19) 



where y = r^/r Q . Denoting £^(E)/D m (E) by /c 2 , it follows 
that 



26 / 



A. M. Weinberg and E. P. Wigner, pp . 192-194, see 
Footnote 23 . 
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Vcp m (r o ,E) 



N(E)/Z®(E) 



x(E)/« + d 0 (E) 



4 D m( E ) p e ( E ) 

l-P e (E) 



where 

X(E) 



K 1 (/cr 1 )I 0 (/cr 0 ) + I 1 (/cr 1 )K 0 (/cr o ) 
K 1 (/cr 1 )I 0 (/cr 0 ) - I 1 (/cr 1 )K 1 (/cr 0 ) 



( 2 . 20 ) 



( 2 . 21 ) 



V Ip K 0 ,, and K-^ are the modified Bessel functions of the 
first and second kind,, respectively. Substitution into 
Eq . (2.19) yields 

l-f (E) 2r 0 (y 2 -l)^(E)P e (E> r „ , 

f(E) TT71) + -y(y -D^xCE) 



r (y 2 -l)Z™(E)d (E) 

+ 2d itr 9 1 - (2 - 22) 

m 



Insertion of Eq . (2.22) into Eq . (2.9) thus attempts to 
account for the heterogeneous effects in the lattice through 
use of diffusion theory in the moderator and transport theory 
in the fuel element. The formulation is not restricted to 
the use of diffusion theory in the moderator region but it 
is utilized because of its calculational ease and good approxi- 
mation to one velocity transport theory. 
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2.3 Formulation In Terms of Escape Probabilities 

The expression, Eq . (2.22), which approximates the 

relative fuel absorption, and hence, the heterogeneous effects 

in the cell may also be arrived at by the formulation of the 

energy dependent absorption rates in terms of moderator and 

fuel element escape probabilities . This approach is an 

extension of Amouyal, Benoist, and Horowitz's one velocity 

27 / 

thermal utilization theory — to the polyenergetic case. 

One begins by defining the following energy dependent 
escape probabilities, all of which are calculated assuming 
the adjacent region is black: 

P e (E) is the escape probability of neutrons incident 
on the fuel element in the actual lattice directional dis- 
tribution, 

P (E) is the escape probability of neutrons in the 
moderator at energy E due to collisions at other energies, 
and, 

P m (E) is the escape probability of neutrons incident 
on the moderator from the fuel element . 

These escape probabilities are schematically depicted 
in Fig. 2.1 for a Wigner-Seitz cylindrical cell. If one 

27 / 

— A. Amouyal, P. Benoist, and J. Horowitz, "Nouvelle 
Methode de Determination du Facteur d ' Utilization 
Thermique d'une Cellule," J. Nucl. Energy 79 
(1957) . 
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Figure 2.1 The Fuel Element and Moderator Escape Prob- 
abilities 
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considers a spatial source of neutrons from scattering which 
has been normalized to one: 



N(r,E)dV = 1 , 



(2.23) 



then reference to Fig. 2.1 shows that the current into the 
fuel element, J_ (E) , may be expressed as 



or. 



AJ_ (E) = P fi (E) + P^(E)P m (E)AJ_(E) , 



AJ_ (E) = 



V E > 



1 -P e (E)P m (E) 



(2.24) 



(2.25) 



where A is the surface area of the fuel element . The net 
current into the fuel element is given by 



J n (E) = J_(E) - P e (E)J_(E) , 



> o that 



AJ n (E) = 



[l-P e (E)]P g (E) 

l-P e (E)P m (E) 



(2.26) 



(2.27) 



Since the source has been normalized to one neutron, Eq . (2.27) 
is the ratio of the absorption rate in the fuel element to the 
absorption rate in the cell. Therefore, 

[l-P e (E)]P s (E) 
l-P e (E)P m (E) * 



f (E) = 



(2.28) 
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Rearranging Eq . (2.28) or substituting Eq . (2.27) into Eq . 
(2.12) yields: 



l-f(E) 

"ner * 



P e (E) 

T^“(e7 



i-P^E) 

T^nr 



i.-p s ( e ) 

V E > 



(2.29) 



Writing the heterogeneous effects in terms of these 
escape probabilities has several advantages . In the first 
place a certain amount of physical transparency is gained by 
separating the contributions of the fuel and moderator into 
easily defined quantities. Secondly, P g (E) is a functional 
of the fuel element properties and its calculation depends 
only slightly on the moderator through the directional dis- 
tribution of the current entering the fuel elements. P C (E) 
and P m (E) are functionals of the moderator and are very 
nearly independent of the fuel element properties . In 
addition, there are no restrictions on this simple model 
and the escape probabilities may be calculated with as much 
sophistication as is desirable. Calculational techniques 
for determining the moderator escape probabilities, P (E) 
and P m (E), are discussed in Chap. 3, while the fuel element 
escape probability, P e (E), is discussed in Chap. 4. 



CHAPTER 3 



MODERATOR ESCAPE PROBABILITIES 
3.1 P s (E) 

In this chapter the moderator escape probabilities, 

P„(E) and P (E), as defined in Sec. 2.3 are discussed, 
s m 

It will be shown that by making appropriate assumptions , 

the formulation of the heterogeneous effects in terms of 

escape probabilities , Eq . (2.29), is identical to the 

expression derived utilizing diffusion theory in the 

moderator and charac terizing the fuel element with its 

escape probability, P g (E), as given in Eq . (2.22). 

The simplest model for P (E) in cylindrical geometry 

s 

is obtained by assuming, as in Sec. 2.2, a spatially flat 

source of neutrons in the moderator. Again as an expediency 

it is assumed that the flux shape in the moderator is 

given by one velocity diffusion theory applied at each 

energy. One finds P (E) by calculating the net current 

into a black cylindrical hole of radius r . For a spatial 

source normalized to one neutron, P (E) is given by 

s 

p s (E) = AD m (E) v cp m (r o , E ) • (3.1) 

Solving the diffusion equation, Eq . (2.13), with the 
boundary conditions given by 

Vcp m (r 1 ,E) = 0 , (3.2) 



- 21 - 
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and. 



yields 



cp m (r Q ,E) 

Vcp m (r o ,E) 



d o< E > ' 



l/[V m Z™(E)] 

Vcp m (r E) - . ( 3 . 3 ) 

X(E)/* + d 0 (E) 

V m is the moderator volume, = Z™(E) /D m (E) , and y(E) is 

given by Eq . (2.21). In accordance with Eq . (3.1), P (E) 

s 

is then given by 

p s (E) - § • , (3.4) 

r o (y 2 -l)/c X (E) + <d Q (E) 

or rearranging,, 

1-P S (E) r 2 o r o 2 

- P /£)' = (y "l) ,c d (E) 4- -J- (y “1) ^x(E) - 1 . 
s v ' 

(3.5) 

As before, y = r p/ r Q * Figure 3.1 shows P g (E) as calculated 
for the Cornell University Zero Power Reactor cores . The 
effect of various water to uranium dioxide volume ratios 
is shown as well as the energy dependence. The parameters 
for the cores utilized in these calculations are given in 
Appendix A, while the cross section data incorporated in 
these computations are discussed in Chap. 7. 



- 









Sec . 3.1 



23 




Figure 3.1 The Moderator Escape Probability, 
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3.2 P (E) 
m ' ' 

The moderator escape probability for neutrons incident 
on the moderator from the fuel element, P m (E), may be cal- 
culated from P_(E) through a reciprocity relationship under 
the following assumptions. If the flux in the moderator is 
isotropic, and if P (E) is calculated using a spatially 

o 

uniform source, as it has been, then P (E) and P (E) are 

s m 

related by ---- 2 - 9 / 



i-VE) 

P s< E > 



ttIH 

’ MjW X . 



where for cylindrical cells 



(3.6) 




= ^(y 2 -l) . 



(3.7) 



P m (E) is shown in Fig. 3.2 as calculated from Eq . (3.6) 

utilizing P (E) as defined by Eq . (3.6) and depicted in 
s 

Fig . 3.1. 

Insertion of Eqs . (3.5), (3.6), (3.7) into Eq. (2.29) 



yields an expression for the relative fuel absorption: 

1-f.m 2r o (y 2 -l)Z>)P e (E) ^ 

f (E) 1-P_ (E) + 2 



+ 



r o(y 2 -l)^( E ) d 0 ( E ) 



2D (E) 
m 



- 1 



(3.8) 



28 / 

- — 7 K. M. Case, F. de Hoffmann, and G. Placzek, Introduction 
to the Theory of Neutron Diffusion , Vol. I, United 
States Government Printing Office, 1953 . 

— ^ K. B. Cady, "Neutron Transport in Cylindrical Rods," 
Massachusetts Institute of Technology Thesis, 1962. 
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Figure 3.2 The Moderator Escape Probability 
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which is identical with Eq . (2.22) as derived in Sec. 2.2. 

This is the expression which will be incorporated in the 
calculations discussed in Chap. 7 in an attempt to account 
for the heterogeneous effects of the lattice. Due to the 
physical transparency and separability of effects achieved 
through the consideration of the various escape probabilities , 
all future discussion will be in terms of them. 



CHAPTER 4 



THE FUEL ELEMENT ESCAPE PROBABILITY 

4.1 The Fuel Rod Escape Probability 

As was noted in Sec. 2.2, it is in many respects desir- 
able to characterize the fuel element in terms of its 
escape probability, P e (E), thereby eliminating the need 
to explicitly define the flux in the various regions of the 
element. In addition, computation of P g (E) allows one to 
complete the calculation of the relative fuel absorption as 
defined in Eq . (3.8). 

In the most general case, the fuel element is composed 
of a cladding region, a void region, and a fuel rod or fuel 
meat region. In the determination of P e (E), the fuel rod 
may be characterized by its escape probability, P r (E) . The 
effect of the void may be accounted for by strictly geometric 
factors, while the cladding effects may be expressed in 
terms of several transmission factors. The results derived 
from the calculation of any of the above factors are of 
course dependent upon the assumption of an angular distri- 
bution of the current entering each respective region. 

The fuel rod escape probability depends on the fuel rod 
radius, a = Z^(E)r^, measured in mean free paths (mfp's) and 
the ratio, c = S^(E)/2^(E) of the fuel scattering cross 
section to total cross section. The simplest expression 
for P r (E) is derived from diffusion theory and it is just 
the albedo of the rod: 



- 27 - 
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l-2D f (E)p(E) 

P r (E) = , 

l+2D f (E)p(E) 

where 



3(E) = /c 



I 1 (/c f r f ) 

f To&pfJ 



(4.1) 



(4.2) 



I and I-, are the modified Bessel functions of the first kind, 
while = 2^(E)/D^(E)j D^(E) being the diffusion coefficient 
of the fuel. The quantity iCf may be obtained either from P-1 
theory—^ 

/c f /s£(E) =V3(l-c) , (4.3) 

or from Wigner's "elementary theory/' 



c tanh~ 1 (/c E /sj:(E)) =/c f /S^(E) . 



(4.4) 



Transport theory expressions for the fuel rod escape 
probability have also been developed from Boltzmann's one 
velocity transport equation under the assumption of 
isotropic scattering in the laboratory system. This results 



in the integral transport or Peierl's equation: 

-2. I r-r ' I 



<p(r) = cp u (?) + cZ t 



dr' e 



cp(r ' ) 



47r|r-r' | 



. i 2 



(4.5) 



30 / 

— A. M. Weinberg and E. P. Wigner, The Physical Theory of 
Neutron Chain Reactors , The University of Chicago 
Pr e s s , Chicago" 1958 . 
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where cp u (r) is the uncollided flux due to incident surface 
sources . Escape probabilities in terms of the solution of 

Eq . (4.5) have been considered by several authors . — 2 ’■ ' 

This work utilizes the treatment by Cady for an isotropic 
flux incident upon the fuel rod. The calculational techni- 
ques utilized in the computation of P r (E) are quoted below.— ^ 
For thin rods of less than 0.5 mfp's a very accurate 
expression results from the uniform flux approximation. In 
terms of Placzek's collision probability, — ' M qo , this yields 
the following expression for P r (E): 

P r (E) = l-2a(l-c) 

is calculated in terms of the modified Bessel functions 
oo 

of the first and second kind with an argument of a: 

M oo - 1 - I ( I l K 1 + l- 2 a [I 1 K 0 + l-a( I 1 K l +I 0 K 0 )> ) • 

(4.7) 



1-M 



oo 



1-cM 



oo 



(4.6) 



For rods which have radii, a, greater than 0.5 mfp's, 
the parabolic flux approximation yields more accurate results . 
In this case ( 



P (E) = l-2a(l-c) { l-(l-c) 



M - c (M M 0 9 - 3M 9 ) 
oo v oo 22 oz 



1-c(4Moo+4M 2 2- 6 Mo 2 )+c 2 (4MooM22- 2 Mo 2 2 )J 



\ 






\ 



(4.8) 



where the following quantities are calculated in order: 



— ! K. B. Cady, "Neutron Transport in Cylindrical Rods," 
Massachusetts Institute of Technology Thesis, 1962. 

— ^ K. B. Cady and M. Clark, Jr., "Neutron Transport in Cylin- 
drical Rods," Nucl. Sci. Eng. lj$, 491-507 (1964). 
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and. 



A = a[I 1 (a)K 1 (a) + I 0 (a)K Q (a)] , 
B = 23(1^ + 1 - A) , 



C = 



= fd^ + 1 - B) , 



D = I (I 2 K o - 1 + B + C) 
4 2I 2 K 1 

E = |(— + C - D) , 

(l -i 4l 9 K 9 

F = y(— j ^ h E ) •> 

a a 



M = 1 - C , 
oo 



M ~ = 1 - 2C + E , 
02 



M 22 = 1 - 3C + 3E - F . 



(4.9) 

(4.10) 

(4.11) 

(4.12) 

(4.13) 

(4.14) 

(4.15) 

(4.16) 



Consideration of the asymptotic expansions for the moments, 
Mq 0 > M q2 , an< * ^22* yields the following approximations for the 
fuel rod escape probability. For the uniform flux approximation. 



P r (E) * = 1 , (4.17) 

r 1 + 2a(l-c) 1 + 22^(E)r f 



and for the parabolic flux approximation, 

P (E) % - 2a . 

l+4a(l-c) + a 2 (l-c) 2 + a 3 (l-c) 4 

(4.18) 

— ^ G. W. Stuart and R. W. Woodruff, "Method of Successive 
Generations," Nucl. Sci. Eng. JLO, 388 (1961). 
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Figure 4.1 shows the fuel rod escape probability as cal- 
culated utilizing the expressions stated above. The effect 



used in the computations are those of the low enrichment., 
high density uranium dioxide fuel of the Cornell University 
Zero Power Reactor. These parameters are listed in Appendix 
A, while cross section data are discussed in Chap. 7. The 
numerical methods needed for calculation of the modified 
Bessel functions are outlined in Appendix C. 



Any void space which is present between the fuel rod 
and the cladding is accounted for by considering the fraction 
of neutrons entering the void space from the cladding which 
are heading toward the fuel rod. The escape probability for 
the fuel rod and the void, P (E), may then be expressed as 



where £ is the fraction of the neutrons headed toward the 
fuel rod. For an assumed isotropic flux emergent from the 



235 

of the 0.27 U resonance is plainly evident. The parameters 



4.2 Void Considerations 



P^E) - (1- O + £P r (E), 



(4.19) 





(4.20) 



— ^ K. B. Cady, see Footnote 31. 



35 / 

— K. M. Case, F. de Hoffman, and G. Placzek, Introduction 
to the Theory of Neutron Diffusion , Vol. I, United 
States Government Printing Office, 1953. 



— ^ K. B. Cady, p. 161, see Footnote 31. 
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Energy , E(eV) 

Figure 4.1 The Fuel Rod Escape Probability, P r ,(E) 
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where is the fuel rod radius and r c is the inner radius 
of the cladding. Examination of Eqs . (4.19) and (4.20) shows 
that the correct limits exist in that P (E) goes to P r (E) 
as r c goes to r and P (E) goes to one as r^ goes to zero. 

4.3 Cladding Effects 

The fuel element cladding is treated in accordance 

37 / 

with the cylindrical thin region theory developed by Cady. — 
The neutron reflection and transmission characteristics of 
the cladding can be divided into several fractions . For 
neutrons incident on the outer surface of the cladding, 

1. g is the fraction of neutrons, which return 
uncollided to the moderator, 

2. T is the fraction transmitted uncollided through 
the cladding, 

3. r. is the fraction scattered in and transmitted 

1X1 

through the cladding, 

4. r . is the fraction scattered in and reflected 

out 

out of the cladding back to the moderator, and 

5. l-(g+T+r^ n +r out ) t * ie rema i n i- n g fraction which 
is absorbed in the cladding. 

Similarly, for neutrons incident on the cladding from the 
fuel rod region. 



37 / 



K. B. Cady, p. 63, see Footnote 31. 
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1. x is the fraction transmitted uncollided through 
the cladding into the moderator, 

2. p^ n is the fraction scattered in the cladding and 
reflected back into the fuel rod and void region, 

3. P out is the fraction scattered in and transmitted 
through the cladding, and, 

4. l"( T+ Pin + Pout^ t ^ e rema i n i n 8 fraction which is 
absorbed in the cladding. 

In terms of these energy dependent fractions, the 
escape probability for the entire fuel element may be 
expressed as 

P e (E) = g+r out + (T+p 0uC )J 0 (E) . (4.21) 



J q (E) is the net current entering the inner surface of the 
cladding from the fuel rod and is a function of P rv (E), 
the escape probability of the void and fuel rod region: 



J 0 (E) 



( T+r in> P rv< E > 

i'PinPrv^ 



(4.22) 



The computational receipes used for the calculation 
of the various fractions mentioned above are stated below. 
These schemes assume isotropic fluxes on both sides of 
the cladding and they are performed at each energy where 
applicable. The outer radius of the cladding is denoted 
by r , the inner cladding radius is r . Macroscopic cross 
sections refer to the cladding material. 
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« ■ r c /r o ’ 


(4.23) 


C ■ r o- r c ’ 


(4.24) 


e - i - i 2 , 


(4.25) 


F - \ sin' 1 % , 

O 7T ^ 


(4.26) 


F l - ? - 


(4.27) 


f 2 - F o + 1 & • 


(4.28) 


f 3 - | % (i-e 2 /3) , 


(4.29) 


n _ tt 6 2 

D i - zr 5 • 


(4.30) 


VV * ln(4/v/e) , 


(4.31) 


°2 - f f 1 ' & ' &A- |)e 2 ' &A- 


y 




(4.32) 


E 1 ’ < F 2 ' 1 D l )/F l ' 


(4.33) 


L 2 = (| F 3 - 20F 1 - 4D 2 )/F 1 , 


(4.34) 


t = f^Ce^Fjl]; + ^(e)^)^ 2 , 


(4.35) 


c i = (jt) ( f F o + ^ 3 


(4.36) 


C 2 = d 2 /| 2 , 


(4.37) 


^1 = 1 C 1 " ^ ’ 

x\ = 29 - | £ - 4C C 2 , 


(4.38) 

(4.39) 



(4.39) 
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t = l-zl(E)r 0 *\ +^(z£(E)r Q ) 2 , 

7 = 0.577215 Euler's Constant, 

x = 2^(E)t , 

7 2 

p u = \ [7 + ln(x) - -j] + 1 - -g- 

X 3 X 4 x 5 

+ 7$ ‘ JZU + JEW ’ 

and, 

P unif 



1 - 



u 

Z S< E > 

Z^(E) 



d-P u ) 



(4.40) 

(4.41) 

(4.42) 

(4.43) 

(4.44) 



For x £ .05 , 



1/-C, 



g = (1-F X ) - 2Z^(E)r 0 (l-F 2 ) + ^(Sj;(E)r o ) 



while for .05 <x 0.2 , 

g = -.0212 - .0266 ln(x)/| 



2.2 



2g( E ) 

r in + r out = 77777 (1 -S- T)P unif 



^(E) 



2 s (E) 

Pin + Pout _ 77777 ^ P unif 



^(E) 



7T (1 ' f 3 ) j 
(4.45) 



(4.45) 

(4.46) 

(4.47) 



f. = — sin -1 

m 7r 



2| 

,nt, 



(4.48) 
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and. 



r in " 


^in^ r in + r out^ J 


(4.49) 


r out 


= < r in + r out> - r in ' 


(4.50) 


Pin = 


^in ^Pin + Pout^ ' 


(4.51) 


Pout 


" ^Pin + Pout^ " Pin ' 


(4.52) 



Calculations have been made using the above techniques 

for the fuel elements of the Cornell University Zero Power 

Reactor, the specifications of which are given in Appendix A. 

The results of these computations in the form of the fuel 

element escape probability are shown in Fig. 4.2. In these 

calculations, the parabolic flux approximation, Eq . (4.8), 

was used for the calculation of P r (E) . Again the effect of 
235 

the U resonances can be seen. 

In terms of previously defined quantities, the ratio 
of the number of neutrons absorbed in the fuel rod to the 
number absorbed in the fuel element is given by 



rod abs . rate 
element abs . rate 



J 0 (E)[l-F rv (E)] 

P rv (E)[l-F e (E)] 



(4.53) 



Similarly the ratio of the number of neutrons absorbed in 
the cladding to the number absorbed in the fuel element 
can be defined by 



Cornell ZPR 
Fuel Element 



38 
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Figure A. 2 The Fuel Element Escape Probability, 
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clad abs . rate 
element abs . rate 



J 0 (E)[l-F rv (E)] 

P rv (E)[l-P e (E)] 



(4.54) 



S a< E > 

2=(E)+E=(E)P U 



[(l-g-T)+J 0 (E)(l-x)]. 



(4.55) 



The results of performing the calculation expressed by 

Eq . (4.53) are shown in Fig. 4.3 for the Cornell fuel 

elements. Inspection of Fig. 4.3 shows that only about 1% 

of the neutrons absorbed by the fuel element are absorbed 

in the cladding. This is because of the low absorption 

cross section of the aluminum cladding. On the other hand, 

38 / 

in the Brookhaven uranium dioxide fuel elements, — (see 
Chap. 8), where stainless steel is used as the cladding 
material, about 12% of the neutrons absorbed in the fuel 
element are lost to the cladding. 



— ^ H. C. Honeck, "The Calculation of Thermal Utilization 
and Disadvantage Factor in Uranium Water Lattices," 
Nucl . Sci. Eng. 18, 49-68 (1964). 
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Figure 4.3 The Ratio of the Fuel Rod Absorption Rate to the Fuel Element Absorption Rate 



CHAPTER 5 



MODELS FOR THE SCATTERING KERNEL 

5.1 The Monatomic Gas Model 

By formulating the lattice heterogeneous effects in 

terms of contrived absorption rates as in Secs. 2.1 and 2.3, 

one is left free to incorporate as much sophistication as 

is desired into the calculation of the energy transfer or 

scattering kernel to be used in Eq . (2.9). With this type 

of an approach one can determine directly the effect of 

different mathematical models on the neutron spectra and 

lattice integral properties . One can then decide on the 

basis of computer time used and the accuracy of the results 

whether or not use of the more complicated models for the 

kernel is desirable for a particular computation. 

Numerous mathematical models for the scattering 

39 40 41 42 43 44/ 

kernel have been proposed . — 1 1 1 1 — ■ In an attempt 

to bracket the degrees of computational difficulty and the 
accuracy of the results predicted with these many models, 
we have limited ourselves to calculations with only two 
of them; the Wigner-Wilkins monatomic gas model— ^ and the 
Nelkin water kernel ^ 

— ! E. P. Wigner and J. E. Wilkins, Jr., "Effect of the 

Temperature of the Moderator on the Velocity Distri- 
bution of Neutrons with Numerical Calculations for H 
as the Moderator," AECD-2275 (1944). 

— ' ^ J. E. Wilkins, Jr., "Effect of Temperature of the 

Moderator on the Velocity Distribution of Neutrons 
for a Heavy Moderator," CP-2481 (1944). 
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The choice of the monatomic gas model is mainly 
dictated by its ease of calculation. This model treats 
the moderator as a free gas in thermal equilibrium at 
the temperature of the moderator. Although crystalline 
or molecular binding effects in the thermalization prop- 
erties of the moderator are neglected, infinite homogeneous 
media calculations— ^ with this model have yielded spectra 
which are essentially consistent with those measured 
experimentally . 

The cross section for scattering from an energy E' 

to an energy E through an angle with a cosine of p may be 

48 / 

derived from the short collision time approximation — 
applied to the single oscillational mode of the trans- 



lational motion of the scattering particle. For the 

:er 

1/2 



49 / 

situation where E' is greater than E, this results in — 



a. 



cr(E' ^ E, p) = -^j- 



EA 



. E'ttTX 



exp 



ZPRi (E ‘ E ' + X> 



(5.1) 



— ^ M. S. Nelkin, "Scattering of Slow Neutrons by Water," 
Phys . Rev. 119 , 741-746 (1960). 

— ^ A. Radkowsky, "Temperature Dependence of the Thermal 

Transport Mean Free Path," Physics Quarterly Report, 
(April, May, June, 1950), ANL-4476, p. 89, Argonne 
National Laboratory (July, 1950) . 

— t H. D. Brown and D. S. St. John, "Neutron Energy 
Spectrum in D 2 O, " DP-33 (1954). 

— I The Scattering Law , Proceedings of the Brookhaven Con- 
ference on Neutron Thermalization, BNL 719, Vol. I, 
Brookhaven National Laboratory (1962) . 
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where o(E' — ^-E,^) is in barns /(eV-unit cosine). The bound 
atom cross section in bams is a^, T is the moderator 
temperature in electron volts , and A is the ratio of the scat- 
tering nucleus mass to the neutron mass. The quantity X 
is related to the momentum transfer in that 

X = E+E' - 2[m/Ie T . (5.2) 



A similar result can be obtained for upscattering in energy 
by employing the principle of detailed balance, Eq . (2.5) 
as is discussed below. 

The scattering kernel is found by the angular inte- 
gration of Eq. (5.1), 

1 



a(E'— *»E) 



a(E '-3-E, p)d|_i . (5.3) 



-1 



This operation can be carried out analytically and it results 
in the formulation of the Wilkins ' heavy gas model— ^ for 
the scattering kernel. For down scattering in energy. 



E. P. Wigner and J. E. Wilkins, Jr., see Footnotes 39 
and 40 . 

— ! M. S. Nelkin, see Footnote 41. 

— ^ L. G. de Sobrino and M. Clark, Jr., "A Study of the 
Wilkins Equation," Nucl. Sci. Eng. 10 , 388 (1961); 
"Comparison of the Wilkins Equation with Experiments 
on Water Systems," Nucl. Sci. Eng. 10, 377 (1961). 

— ^ M. S. Nelkin, see Footnote 41. 

— ! F. D. Federighi and D. T. Goldman, "KERNEL and PAM - 
Programs for Use in the Calculation of the Thermal 
Scattering Matrix for Chemically Bound Systems," 
KAPL-2225, Knolls Atomic Power Laboratory (1962). 

— ^ J. E. Wilkins, Jr., see Footnote 40. 
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E' > E, the monatomic gas scattering kernel is expressed in 



barns /eV as 



51/ 



a(E 1 — *-E) = 



where 



and. 



a f 7] 

2E ' 



2 / 



ERF(r)Z-pz') + ERF(rjz+pz') 



+ exp(z ' 2 -z 2 ) [ERF(r]Z ' -pz) -ERF(pz+rjz ' ) ] 



z Z = E/T , 
z' 2 = E'/T , 

V = -j(A+l)A ly/2 , 

p - -|(A-1)A" 1/2 . 



(5.4) 



(5.5) 

(5.6) 

(5.7) 

(5.8) 



ERF(z) is the error function defined by 

z 



ERF(z) 



h'l 







O 

while Of is the free atom scattering cross section 
moderating medium. 



(5.9) 
for the 



a f 



% 




(5.10) 



The scattering kernel for upscattering, E' < E, is cal- 
culated utilizing detailed balance as expressed in Eq . (2.5). 

In the case where one is considering moderation by 
hydrogen, Eq . (5.4) reduces to the Wigner-Wilkins model for 
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52 / 

a gas of mass one. — Then for E' 



a(E'— »-E) = 



> E, 

ERF(z) . 



(5.11) 



The use of the monatomic gas model has another advantage 



in that the scattering cross section defined by 

‘00 



a s ( E ) " 



c(E-^E')dE' , 



(5.12) 



is an analytic function given by 

1 



a g (E) — 



(1 + -4^) ERF(z) + - 

2y J^ry 



(5.13) 



where 




(5.14) 



This analytic function for the scattering cross section is 
useful because it allows one to check on the accuracy of 
the numerical integration scheme employed. This aspect is 
more fully discussed in Sec. 5.2. 

The calculations performed utilizing the monatomic 
gas model for the scattering kernel are discussed in 
Chap. 7. Calculation of the error function necessary in 
the computations of this kernel is detailed in Appendix C. 



— M. J. Poole., M. S. Nelkin, and R. S. Stone, "The Measure- 
ment and Theory of Reactor Spectra, " Progress in 
Nuclear Energy , Series I, Vol. 2, pp. 91-164, Pergamon 
Press, New York, London, Paris, and Los Angeles, 1958. 

52/ 

— ' E. P. Wigner and J. E. Wilkins, Jr., see Footnote 39. 
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5.2 The Nelkin Water Kernel 

The second mathematical model for the scattering kernel 

53 / 

utilized in these calculations is that proposed by Nelkin — 
for protons bound in water. Although the computer cal- 
culation of the Nelkin water kernel takes a relatively 
long time, its use is justified in that infinite homogeneous 
media thermal spectra calculated with this model are in good 
agreement with those measured experimentally In addition, 

diffusion coefficients calculated with the Nelkin model— ^ 
agree well with measured values ' 

Physically, the mathematical approximation proposed by 
Nelkin for water treats the scattering molecule translation 
with the free gas model utilizing a translator mass (M) of 
18. The effects of the atomic motions within the molecule 
are approximated by simple harmonic oscillations, the 
phonon frequencies of which determine the relative strength 
of the interaction with a particular oscillator. The 
oscillators treated are as follows: A hindered rotational 

oscillator of mass (M ) 2.32 and frequency (O^) corresponding 
to an energy (oi = *fiO ) of 0.06 eV, where'll is Planck's 
Constant divided by 2 tt ; a vibrational oscillator of 
mass (My) 5.84 and an energy (a>^) of 0.205 eV ; and a 

— / M. S. Nelkin, see Footnote 41. 

— ^ J. R. Beyster, et al., "Measurement of Neutron Spectra 
in Water, Polyethylene, and Zirconium Hydride," Nucl. 

Sci. Eng. £, 168 (1961). 
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second vibrational oscillator of mass (M^ = M v /2) 2.92 and 
energy (ac^) of 0.481 eV. Different mathematical approxi- 
mations may be made to treat the translational, vibrational 
and rotational modes, dependent upon the incident neutron 
energy and the energy or momentum transfer. However, a 
single expression may be used to characterize the differential 
scattering kernel for scattering from an energy E' to an 
energy E (E 1 > E) through an angle with a cosine of p,. 
Parameters within the Nelkin water kernel then change 
dependent upon the specific approximation made. 




m=0 



ZZ 








nao 

2t 

e I n ( z ) exp 




(E-E' 



+ 



X 

M 



+ mx>2 



+ £a>i 



+ nao) 




(5.15) 



— ^ G. P. Calame, "Diffusion Parameters of Water for Various 
Scattering Kernels," Nucl. Sci. Eng. _19, 189 (1964). 

— ^ E. Starr and J. Koppel, "Determination of Diffusion 

Hardening in Water," Nucl. Sci. Eng. 14, 224 (1962). 
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Equation (5.15) gives the differential scattering kernel in 

barns /(eV-unit cosine) with as the bound atom cross 

section in barns . The I (z) are the modified Bessel 

n 

functions of the first kind and order n. As before, the 
quantity X is given by 

X = E+E'-2|r \f~EE' . (5.16) 

One may allow for all of the transitions for the 
translational and rotational modes, while utilizing a 
phonon expansion for the first (co^ = .205 eV) vibrational 
mode and the elastic (unexcited) limit for the second 
(o>2 = .481 eV) vibrational mode by inserting the following 
values into Eq . (5.15): 

I = I + -J_ + _L_ 

A co r M r tanh(o) r /2T) M v“l *V°2 



£ = T/M , 



z 



X 

Mcosinh(oo/2T) 



1 

H 




9 



mm = 0 , 

m = 2 , 

CO = 0) r , 



(5.17) 
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and 



M = M. 



r 



where T is the moderator temperature in eV. The phonon 
expansion for the second vibrational mode is included if 
mm is set equal to 2 in the above set of parameters . 

To allow for all transitions of the translational 
and first vibrational modes , while treating the rotational 
mode with the short collision time approximation and the 
second vibrational mode with a phonon expansion, the 
following parameters are inserted into Eq . (5.15): 



1 . 1+2 



A V>2 



? M ’ M r e to r /T_ 1 ’ 2 




z 



X 



Mds inh (cd / 2T ) 



(5.18) 



1 1 
M " M 



= 1 + A 

M M r 



mm = 2 



ii = 0 



CD as oo 



1 



and 



M = M v . 
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The oxygen present in the water is accounted for by adding 
the monatomic gas model with mass 16., Eq. (5.1), to the 
Nelkin water kernel, Eq . (5.15). 

Angular integration of the differential scattering 
kernel, a(E' — *»E,|j.), was performed numerically to obtain 
the scattering kernel in the form utilized in lattice cal- 
culations. For the non diagonal terms in the scattering 
matrix (E'=j=E), the Gauss quadrature formula was used for 



the angular integration: 

1 



a (E ' — -E) = 



kk 



a(E ' — »-E, |a)dp. 



-1 



k=l 



w. 



a(E' — ^E,(i k ) + a(E' — *-E,-|i k ) 



(5.19) 



where w k is the weighting coefficient corresponding to p, k , 
a zero of the Legendre polynomials For terms where 

the energy transfer is large (E « E'), an eight point 
(kk = 4) integration is used, while for small energy 
transfers where the fine structure of the differential 
kernel is a rapidly varying function of angle, a sixteen 
point (kk = 8) quadrature integration is performed. 

Calculation of the scattering kernel diagonal terms 
(E = E'), where there is no energy transfer, requires some 

— ^ A. N. Lowan, N. Davis, and A. Levinson, "Tables of the 
Zeros of the Legendre Polynomials of Order 1-16 and 
the Weight Coefficients for Gauss's Mechanical Quad- 
rature Formula," Bull. Am. Math. Soc . 48, 739 (1942). 
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modification of the above procedure. Inspection of Eq . 
(5.15) shows that for E' = E, X = 2E(1 -p), and there is 
a singularity in the differential scattering kernel at 

H = 1. 

The use of the Gaussian quadrature angular integration 

scheme introduces numerical errors in this situation. 

Therefore the following calculational technique is utilized 

58 / 

as suggested by Federighi and Goldman. — ' Equation (5.15) 
is considered in the limit as X goes to zero (p goes to 1) . 
The singularity is subtracted from the differential kernel 
and integrated analytically, the remainder being integrated 
with the sixteen point Gaussian quadrature rule. The 
analytically integrated singularity is then added to the 
result. This process is indicated below: 



and X = 2E(l-p). The diagonal terms are then approximated 
by 



1 



1 




(5.20) 



where. 



f(n) = c(E— E,n) - “7= , 

4/71x 



(5.21) 



. 58/ 



F. D. Federighi and D. T. Goldman, see Footnote 49. 
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a(E — s-E) 



8 




f (^ k )+f (-|a k ) 



k=l 



+ 



g b 

2ME 



(5.22) 

In the calculations of the Nelkin water kernel made for 
incorporation into the lattice studies Eqs . (5.17) were 
utilized for incident energies , E', less than 0.33 eV, 
while Eqs. (5.18) were used for E' greater than 0.33 eV. 

This is essentially the same approximation made by Honeck 

59 / 

for the water kernel calculations incorporated in THERMOS. — 
The phonon expansion for the second vibrational mode along 
with Eqs. (5.17) is not used because only a small increase 
in accuracy at higher incident energies is accomplished 
with a large (by a factor of 3) increase in computer time 
through its utilization. 

Figure 5.1 shows the results of a calculation which 
utilized both Eqs. (5.17) and (5.18) in Eq . (5.15) for an 
incident energy of 0.366 eV. The results are expressed in 
barn-s/unit lethargy, Ea(E ' — *-E), for down scattering, 

E' > E. However, the basic difference between the two 
approximations is plainly evident in Fig. 5.1. By con- 
sidering all transitions for the hindered rotator, Eqs. 
(5.17), a scattering kernel is predicted with a complicated 
fine energy structure having oscillations at energy 

— I H. C. Honeck, "THERMOS - A Thermalization Transport 
Theory Code for Reactor Lattice Calculations," 

BNL 5826, Brookhaven National Laboratory (1961) . 



Ea(E' — (barns /unit lethargy) 
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’igure 5.1 Approximations Used in the Calculation of the 

Nelkin Water Kernel 
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intervals (E f — E) equal to the frequency of the rotator, 

<*> r = 0.06 eV. In going to the short collision time approxi- 
mation for the hindered rotator, Eqs . (5.18), this fine 
structure is considerably smoothed, however the cross 
section derived by the energy integration of the kernel, 

Eq . (5.12), is essentially the same for both models. The 
choice of 0.33 eV as the energy above which Eqs. (5.18) will 
be used and below which Eqs. (5.17) will be used is made 
with this fact in mind, but primarily on the basis of 
computer economy. In order to predict accurate cross 
sections utilizing all transitions for the hindered 
rotator, it is necessary to use a fine energy mesh in the 
calculations (AE < .06eV), but if one is interested in 
energies in the range up to 1.0 eV this means that an 
excessively large number of energy points must be used 
in the calculation. This problem is further discussed 
below and in Chap . 7 . 

The first moment of the scattering kernel, a-^(E' — »-E) 
is calculated for down scattering utilizing the same 
techniques as above, where. 



As before, a sixteen point Gaussian quadrature angular 
integration is made for the off diagonal terms when the 
energy transfer is small, while eight points is considered 



1 




(5.23) 
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sufficient for large energy transfers. The diagonal terms * 
E' = E* for the first moment are calculated using the 
function defined in Eq . (5.21): 



a ] _(E 



— E) 




(if(tx)d^ + —— 

4 fwi 



1 




V 2E(l-(i) 



Wi 



n k f (^ k )-w k f ^ a k> 



k=l 






(5.23) 

The upscattering terms for both the scattering kernel 
and its first moment are calculated using detailed balance* 

Eq . (2.5). It is to be noted that a slight* but 
acceptable inconsistency is introduced in this manner* 
since the use of detailed balance means that in the up- 
scattering portion of the scattering matrices* the 
approximation used* Eqs . (5.17) or Eqs . (5 . 18)* actually 
depends on the neutron final energy rather than its initial 
energy . 

Energy integration was performed on the scattering kernel 
and its first moment to determine the average cosine of 
the scattering angle* ]I(E) : 



h( E) 




a 1 (E ^E')dE' 



o( E ^E)dE' 



(5.24) 



o 
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A trapezoidal rule numerical integration scheme. Appendix C, 
is used for this calculation and the results are shown in 
Fig. 5.2. Figure 5.2 also shows the transport cross section 
of water as calculated using the expression 

a tr (E) = a a (E) + a s (E)[l-H(E)] . (5.25) 

In the energy integration of the scattering kernel, 

Eq. (5.12), to obtain the scattering cross section, a (E), 

s 

certain numerical inaccuracies were encountered. Specifically, 
numerical integration of the Nelkin water utilizing the 
trapezoidal integration scheme discussed in Appendix C con- 
sistently yields scattering cross sections which were higher 
than those measured experimentally— ^ for energies greater 
than 0.1 eV. Essentially the same percentage error is seen 
when the numerically integrated monatomic gas kernel is com- 
pared with the analytic expression for the cross section 
given by Eq. (5.13). These numerical errors are attri- 
butable to the fact that at higher incident energies, the 
diagonal terms of the scattering matrix are large compared 
with the first off diagonal terms for the energy mesh used 
in the calculation as discussed in Chap. 7. Use of the 
trapezoidal integration rule in this region therefore yields 
scattering cross sections for both kernels which are higher 

— ^ D. Hughes and B. Schwartz, "Neutron Cross Sections," 

BNL 325, 2nd ed., Brookhaven National Laboratory 
(1958) . 
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.Figure 5,2 The Average Cosine of the Scattering Angle and the Transport Cross 

Section for Water 
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than experimentally measured or analytically calculated 
values . Auxiliary calculations showed that good agreement 
between cr g (E) calculated with the Nelkin water kernal and 
measured experimentally or between a (E) calculated by the 
monatomic gas model and a (E) given by Eq . (5.13) could 

U 

be obtained by either utilizing a much finer energy mesh 
or by using more sophisticated integration schemes. However, 
both of these approaches create large increases in the com- 
putational time since the kernel calculation time is pro- 
portional to n where n is the number of energy points con- 
sidered, and the more complicated integration schemes which 
involve the fitting of exponentials or higher order 
polynomials to the calculated kernel points require 
additional computer time due to the uneven energy spacing 
used in the calculations . It is therefore desirable to 
retain the simplicity of the trapezoidal rule and correct 
for its inability to handle accurately the large diagonal 
terms encountered at high energies . This is done by 
recognizing that the actual value of the diagonal terms, 

E' = E, has no effect on the spectrum calculated with 
Eq . (2.1) or Eq . (2.9). This is seen by substituting the 
definition of the scattering cross section, Eq . (2.2) 
into Eq. (2.1) or (2.9) and letting E go to E* . Further- 
more, integration of both sides of Eq . (2.1) or Eq . (2.9) 
over the energy variable E indicates the scattering cross 
section used in calculations must be the one derived from 
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the scattering kernel by the numerical integration scheme. 

Since some of the integral properties that are to be 
calculated are strongly dependent upon the use of an 
accurate scattering cross section, we propose to modify 
the calculated diagonal terms in such a way that the 
scattering cross section calculated with the trapezoidal 
numerical integration scheme will equal the measured values 
for the Nelkin water kernel and Eq . (5.13) for the monatomic 
gas model. Thus, the diagonal terms, a*(E— ^-E), utilized 
in the calculations will be given by: 

a*(E— *“E) = C(E)o(E — ^E) , (5.26) 

where a(E — »-E) is calculated using Eq . (5.22). The correction 
factors, C(E), are easily found for the trapezoidal inte- 
gration rule from the difference between the measured or 
analytic scattering cross section and the scattering cross 
section using the diagonals, a(E — ^E). These correction 
factors are assumed to be temperature independent and those 
calculated for room temperature may therefore be applied to 
the diagonals, a(E— *-E) calculated at any temperature. It 
is to be noted that this technique is not applied to the 
calculation of ]I(E), Eq . (5.24), since the most important 
effect is the relative magnitude of the correctly calculated 
values of c^(E' — *-E) and a(E' — *"E) . 

Appendix B gives a listing of the scattering matrix 
as calculated using the Nelkin water model. The physical 
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temperatures of 20°C and 40°C are considered. 

The source of thermal neutrons, S(E), due to down 
scattering from energies above the thermal cutoff energy, E c , 
can be calculated in accordance with Eq. (2.3) by taking the 
short collision time approximation on all of the modes of 
oscillation in the Nelkin model. Assuming a 1/E flux above 
E c and neglecting upscattering, Eq. (2.3) yields an un- 
normalized source of: 

S(E) = ERF/I7T e , (5.27) 



where ERF(z) is defined in Eq . (5.9). 

The effective temperature, T , is given by: 

T T , E r E 1 2E 2 
i e “M + M r + M v + M v 



(5.28) 



where the E's are given for each mode of oscillation as: 



E = 



/ 





(5.29) 



Because of the relative magnitudes of and o^, T £ may be 
approximated by: 



T 



e 






(5.30) 



At a physical temperature of 20 °C the effective temperature 
has a value of 0.1169 eV . 



CHAPTER 6 



THERMAL LATTICE SPECTRA AND INTEGRAL PROPERTIES 

Having selected appropriate mathematical models for the 
scattering kernel, Eq. (5.4) or Eqs . (5.15) and (5.19), and 
for the relative fuel absorption, f (E) / [ 1-f (E) ] , Eqs. (3.8) 
and (4.21), one proceeds to look for solutions to Eq . (2.9), 
which is rearranged below. 



S (E) +/ S(E'^-E)cp m (E' )dE' 



?“(E) = 



2“(E) + 2“ l (E) 



-,m. 



1 + 



f(E) 



1 - f(E) 



( 6 . 1 ) 



where S(E) is given in its unnormalized form by Eq . (5.27) and. 



Z(E'— *-E) = Na(E'— E) , (6.2) 

N being the number density of the moderator. Solution of 
Eq . (6.1) for the spatially averaged, energy dependent thermal 
flux in the moderator is accomplished with Program COUTH 
(Cornell University Thermalization) , a Fortran 63 code written 
for use with the Control Data Corporation 1604 digital com- 
puter. A Fortran statement of COUTH appears in Appendix E, 
while the numerical methods employed are discussed in 
Appendix C. This chapter describes the lattice variables and 
integral parameters calculated in COUTH. The results of a 
set of actual lattice calculations are described in Chap. 7. 
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Determination of the spatially averaged flux in the 
moderator through the solution of Eq . (6„1) leads to the 
calculation of several other useful quantities. The 
spatially averaged thermal flux in the fuel, cp^(E), and in 
the cladding, cp^E), may be determined from the moderator 
spectrum by equating absorption rates in the various 
regions. Utilizing the relationships in Eqs . (4.53) and 
(4.54), 



-f , 



V m 


' f(E) 




V f sf(E) 


_l-f(E) 





J 

P 



rv (E)(l-P e (E)) 




(6.3) 

m f 

where V and V are the volumes of the moderator and fuel 
respectively, S^(E) is the fuel absorption cross section, 
5“f?V(l-f) is the absorption rate in the fuel element, 

a 

and J (1-P ) /(1-P )P is the ratio of absorptions in the 

fuel meat to absorptions in the fuel element. Similarly 
for the cladding: 



?°(E) 



v m 


f(E) 




V C 2®(E) 


l-f(E) 





J o( E )< 1 - p ry( E >> 

P rv (E)(l-P e (E)) 







(6.4) 

Q C 

where V is the cladding volume and 2 (E) is the cladding 

a 

absorption cross section. 

Average neutron velocities in each region of the cell 
may now be defined utilizing the relationship between 
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neutron density and flux, 

N(E) = ?(E)/v(E) , (6.5) 

where v(E) is the neutron velocity. The average velocity in 
the moderator, v m , is then expressed as 



v(E)N m (E)dE 



9 m (E)dE 



v m E 



( 6 . 6 ) 



N m (E)dE 



9 m (E)/v(E)dE 



Similar expressions define the average velocities in the fuel, 
and in the cladding, v £ . 

Definitions of disadvantage factors follow directly by 
expressing the total spatially averaged neutron density, N, 
as shown below for the moderator: 



if 1 



N m (E)dE = / 9™(E) /v(E)dE 






(6.7) 



Using analogous expressions for the total densities in the 
fuel, N r , and the cladding, N c , one can calculate neutron 
density disadvantage factors for the moderator, 6™; 

8* = lf7N f , (6.8) 

and for the cladding 5^; 

6^ = N° /N f . (6.9) 
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The disadvantage factors are related to the familiar flux 
disadvantage factors through the relationships: 

6® =^/T f = v m 5™/v f , (6.10) 

and 

6 9 = * ^ c 6 n^f J (6-11) 



where m, c, and f again refer to moderator, cladding, and 
fuel, and 



-m 

9 



V 11 (E)dE 



( 6 . 12 ) 



In terms of quantities already defined the thermal 
utilization, f, and the more basic quantity, qf, may be 
determined, where 



and. 



thermal absorption rate in the fuel rod 

f = 

thermal absorption rate in the cell 



(6.13) 



thermal fission neutron production rate in the cell 

T)f = • 

thermal absorption rate in the cell 

(6.14) 

These definitions lead to the following expressions: 
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f 




zf(E) cp f (E)dE 



Z”(E) 



1 + 



m 



m 



cp m (E)dE 



(6.15) 



and. 



v V £ l Z £ (E) cp f (E) / [ l+a(E) ]dE 



Vf = 



(6.16) 



V 



.m 



c 

s ~ 


-1 


/ E^(E) 


1 + f .(E) 


1 + l-f(E) 



cp m (E)dE 



where v is the number of fission neutrons produced per 
thermal fission, and a(E) is the energy dependent capture 
to fission ratio for the fuel. 

A thermal neutron lifetime, may be defined as 

the ratio of the total number of thermal neutrons in the 
cell to the thermal absorption rate in the cell, 



tfV™ + + N f V f 



£ t ~ 



(6.17) 



V 1 



m 



s i ;‘(E) 






W) 



9 m (E)dE 



where N is defined by Eq . (6.7) for each region- 
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Since one of the prime concerns in reactor calculations 
is the accurate determination of various reaction rates,, it 
is desirable to define average thermal cross sections which 
are suitable for use in calculations utilizing the thermal 
flux or the 2200 meter/second flux. In terms of the 
spatially averaged flux or neutron density in a given region, 
a reaction rate, RR, per unit volume of the region is 
expressed by: 



E 



RR = 



2(E)cp(E)dE = 



2(E)N(E)v(E)dE 



(6.18) 



As before, superscripts on the variables will denote the 
region considered (m = moderator, c = cladding, f = fuel), 
while subscripts will denote the interaction considered 
(a = absorption, s = scattering, t = total, tr = transport, 
f = fission) . Thus, the thermal absorption rate in the 
moderator will be given by 



RR m = 
a 



2®(E)9 m (E)dE = 



z“(E)Tf"(E)v(E)dE . (6.19) 

a 



In the notation of Wescott,— ^ reaction rates may now be 
formulated in terms of an effective cross section, 2, and 



— I C. H. Wescott, "The Specification of Neutron Flux and 
Effective Cross Sections," AECL 352 (1956). 
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a mean cross section, Z, for use in a calculation with the 
spatially averaged 2200 meter per second flux, 9 , and the 
spatially averaged thermal flux, cp T , respectively. A 
given reaction rate is now defined by 



RR = Zcp Q = 2cp T , 



( 6 . 20 ) 



where, 



<p^ = Nv = v 
o 00 



N(E)dE = v Q / <p(E)/v(E)dE , 



( 6 . 21 ) 



and, 



cp T = N v = / N(E)v (E)dE = / cp(E)dE . (6.22) 



The velocity v Q denotes 2200 m/sec, while v is the mean 
velocity defined by Eq . ( 6 . 6 ). This leads to the definition 
of the effective cross section as: 



Z(E)cp(E)dE 



Z = 



vj N(E)dE 



(6.23) 



while the mean cross section is defined as 
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2 = 



2(E)cp(E)dE 




(6.24) 



Consideration of Eqs . (6.23) and (6.24) shows that Zv Q = Zv. 
Applying this nomenclature to absorption rate per unit 
volume of the moderator as an example yields: 




— m 

= 2 cp 

a 



— m 

S a 



( 6 . 25 ) 



where 



/ 2®(E)^ m (E)dE 

o a 

v 0 J 9 m (E)/v(E)dE 





E c 

/ 5®(E)V“(E)dE 

o 

E c 

/ cp m (E)dE 
o 

^ n (E)/v(E)dE , 




(6.26) 



(6.27) 



(6.28) 



o 
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and. 



Z™ v = Z m v 
a o am 



( 6 . 29 ) 



( 6 . 30 ) 



This approach is extended to define effective and mean 

-h 

cross sections for the homogenized cell, Z n and Z , sub- 
scripts denoting the interaction considered. In this 
manner 



dE 



h -'o 



V m 2 m (E)cp m (E) + V C Z C (E)cp C (E) + V f Z f (E)cp f (E) 



Z“ = 



dE 

o j vTEJ 

o 



v”V"(E) + vV^E) + V £ <p f (E) 



( 6 . 31 ) 



and. 



dE 



s h 



V^CE^CE) + V c X c (E)cp c (E) + V £ 2 £ (E)<p £ (E) 



dE 



V I V n (E) + V c cp C (E) + V f cp f (E) 



( 6 . 32 ) 
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The mean velocity of neutrons in the homogenized cell is de- 
fined as 



dE 



vV© + vV^E) + vV(E) 



,c — c 



f-f , 



V, = 



V, 



dE 

vTEj 



VV'(E) + V^CE) + V f cp f (E) 



(6.33) 



so that. 




(6.34) 



To complete the calculation of the thermal parameters * 

2 

one would like to define a thermal diffusion area, L , for 

the homogenized cell in a meaningful way. In the one 

2 

velocity, one region problem, L may be defined unambiguously 
as 

L 2 = D/2 , (6.35) 

a 

where D = 1/32^^.. However, in the polyenergetic, homogenized 

2 

cell problem, one must decide whether L in Eq. (6.35) should 
be calculated utilizing mean or effective cross sections, and 
whether D should be calculated by averaging the transport 
cross section or the transport mean free path. 
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6 9/ 

Larrimore— 'has indicated the desirability of utilizing 

U 

the mean homogenized cross section, 2 :, in the calculation 

a 

0 AQ A/ £ c / 

of L . He also notes that several authors — 1 1 ' argue in 

favor of averaging the transport mean free path, while 
Sjbstrand— ^ has shown a physical basis for averaging the 
transport cross section in the determination of D. The paradox 
between these two approaches has recently been clarified by 
Pomraning— ^ and on the basis of that discussion this work 
averages the reciprocal of a homogenized transport cross 
section, 2 “ r (E), over a homogenized spectrum in the deter- 
mination of the diffusion coefficient, D. The homogenized 
transport cross section is defined as 



/t , n _ 2® r (E)7 n (E)V m + E^ r (E)q) C (E)V C + 2^ r _(E)cp f (E)V f 
^tr^) ~ -in /TlX „m , — c /t ,\ttC , — f/T^\T,f 



tr' 



9 U (E)V + 9 (E)V^ + 9 (E)V 



(6.36) 



while the homogenized flux, 9 (E), is defined by 



iV) = 



— Ill / ■p \ t rill 

9 (E)V 



+ 



— c 

9 



(E)V ( 



v m + V C + V 



9 f (E)V f 



(6.37) 



— J. A. Larrimore, "Temperature Coefficients of Reactivity 
in Homogenized Thermal Nuclear Reactors," Massachusetts 
Institute of Technology Thesis (1962) . 

— ^ A. M. Weinberg and E. P. Wigner, The Physical Theory of 
Neutron Chain Reactors , p. 514, The University of 
Chicago Press, Chicago, 1958. 

— t C. D. Petrie, M. L. Storm, and P. F. Zweifel, "Calculation 
of Thermal Group Constants for Mixtures Containing 
Hydrogen," Nucl. Sci. Eng. 2, 728-744 (1957). 
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The diffusion coefficient is then expressed as 

? h (E)/s£ r (E)dE 

D= 2_ . (6 .3 8) 

? h (E)dE 

The thermal diffusion area is calculated from 

L 2 = D/fO* , (6.39) 

a 

where 2^ is defined by Eq. (6.32). 

The calculation of all of the parameters and variables 
discussed in this chapter is carried out by program COUTH. 
The results derived from these calculations for an actual 
lattice system are discussed in Chap. 7. 





— K. S. Singwi and L. S. Kothari, "Transport Cross Section 
of Thermal Neutrons in Solid Moderators/' Second 
Geneva Conference Proceedings, 16^ 325 (1958). 

— ! N. G. SjOstrand, "Definition of the Diffusion Constant 
in One-Group Theory," J. Nucl. Energy, Part A, 

Reactor Science, 151 (1960). 

— ^ G. C. Pomraning, "On the Energy Averaging of the Diffusion 
Coefficient," Nucl. Sci. Eng. 1£, 250 (1964). 



CHAPTER 7 



APPLICATION OF THE METHOD TO ACTUAL LATTICES 

7.1 Room Temperature Calculations 
In this chapter we discuss the numerical results obtained 
with program COUTH when the cell theory formulated in Secs . 

2.1 and 2.3 is applied to an actual lattice system. The aim 
of this work is to calculate at various temperatures , the 
spatially averaged thermal spectra, and the integral parameters 
defined in Chap. 6 for all of the cores incorporated in the 
Cornell University Zero Power Reactor facility The 
Cornell ZPR cores are light water moderated and utilize high 
density, low enrichment uranium dioxide fuel, clad in aluminum. 
Several water to uranium oxide volume ratios (nominally 1:1, 
1.5:1, 2:1, 3:1 and 4:1) are provided with sets of grid plates 
which vary the spacing (lattice pitch) between the fuel 
elements arranged in a hexagonal array. More specifically, 
the fuel rods consist of 0.60 inch diameter uranium dioxide 
fuel of density 10.38 gm/cc, enriched to 2.096 atom per cent 
of U . The fuel rods, formed by 48 one inch long fuel 
pellets, are clad in 0.028 inch thick aluminum to form a 
fuel element with an outside diameter of 0.666 inches. A 
list of the parameters of the various lattices is given in 
Appendix A. 

— t D. D. Clark, "The Cornell University Nuclear Reactor 
Laboratory," CURL-1, Cornell University (1961). 
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Initial calculations have been carried out utilizing 
as a scattering kernel both the Nelkin water kernel, 

Eq . (5.15), and the monatomic gas model, Eq . (5.11), for 
a gas of mass one. The two physical temperatures of 20°C 
and 40°C have been considered. The results of an elementary 
study utilizing the monatomic gas model have previously been 
reported ^ 

In the solution of the neutron balance equation, 

Eq . (6.1), to determine the spatially averaged moderator 
flux, ^(E), the heterogeneous effects of the lattice are 
accounted for by the relative fuel absorption, f (E) / [ 1-f (E) ] , 
as given in terms of escape probabilities by Eq . (2.29). 

The moderator escape probabilities, P g (E) and P m (E), are 
given by Eqs . (3.4) and (3.6), respectively. This yields 
Eq . (3.8) for the relative fuel absorption. The fuel element 
escape probability to be used in Eq . (3.8) is calculated by 
means of Eq . (4.21), which in turn requires the calculations 
indicated by Eqs. (4.22), (4.19), and (4.8). Equation (4.8) 
calculates the fuel rod escape probability utilizing the 
parabolic flux approximation. 

— ^ S. S. Berg, "Initial Experiments on the Cornell University 
Zero Power Reactor Cores," Cornell University Thesis 
(1964). 

— t K. B. Cady and C. R. Mac Vean, "A Simplified Formulation 
of Space-Energy Cell Theory," Trans. Am. Nucl. Soc . 6_, 
234-235 (1963). 
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The source of neutrons due to down scattering from 
energies above E c is that given by Eq . (5.27), normalized 
to one neutron per second per unit volume of the moderator. 
That is, 



The calculations were performed utilizing an energy 
mesh containing forty- five energy points ranging from 
0.001 eV to a thermal cutoff energy, E c , of 1.609 eV. In 
the stipulation of a useful energy mesh a compromise must 
be reached between several conflicting requirements . In 
the first place, a reasonable range of energies must be 
covered. Secondly, the total number of energy mesh points, 
n, may not be too large, since the computer time necessary 
for the calculation of the scattering kernel increases as 
n squared. Furthermore, it is desirable to have a finely 
spaced energy mesh in the region above the thermal energy, 
E Q = kT, where the 1/E tail joins that part of the spectrum 
resembling a Maxwellian. Finally, the use of the Nelkin 
water kernel places certain restrictions on the size of 
the energy spacing, AE, between the energy points in the 
mesh. With the use of the approximation allowing for all 
transitions of the hindered rotator, Eqs . (5.17) (utilized 
for incident neutron energies less than 0.33 eV) it is 



erf/eTt 



e 



S(E) - -j. 




(7.1) 



e 



o 
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necessary to use a AE which is less than the frequency 
(oo^ = 0.06 eV) of the hindered rotator. This requirement 
insures that the fine structure of the water kernel is 
preserved in the calculations and that energy integration 
of the scattering kernel will yield a fairly accurate 
scattering cross section. In order to satisfy this require- 
ment, the energy spacings utilized in the calculations are 
less than one quarter of cd^ or .015 eV for all energies 
below 0.33 eV. In light of these restrictions, the energy 
mesh used in the calculations consists of 35 energy points 
below E = .33 eV with 34 energy intervals corresponding 
to equal intervals in velocity. Each of these low energy 
velocity intervals corresponds to 0.1 v Q or 220 m/sec. 

For energies above E = 0.33 eV the ten energy intervals 
correspond to velocity intervals which exponentially 
increase from 220 m/sec to 2200 m/sec as energy increases. 

A complete listing of the energy mesh utilized is given 

in Appendix B . 

23 5 

The U absorption and fission cross sections used 
are those compiled by Suich.— ^ The absorption cross 

no o 

sections for U , the moderator and the cladding are all 
assumed to be 1/v and the 2200 m/s values of these cross 
sections are given in Appendix A. The scattering cross 

— f J. Suich, "Temperature Coefficients in Heterogeneous 
Reactor Lattices," Massachusetts Institute of 
Technology Thesis (1963) . 
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sections for all of the materials other than the moderator 

are calculated utilizing the analytic expression, Eq . (5.13) 

derived from the monatomic gas model. Scattering cross 

sections for the moderator are computed by integration of 

the scattering kernel in accordance with Eq . (5.12). The 

use of the monatomic gas model yields a water scattering 

cross section given by Eq . (5.13), while the use of the 

Nelkin water kernel results in the scattering cross section 

72/ 

given by BNL 325. — The free atom scattering cross sections 
used are given in Appendix A. 

The transport cross sections are calculated using 

E tr (E) = 2 a (E) + V E) [1 ' ^ (E)] * (7 ‘ 2) 

The average cosine of the scattering angle, p(E), is cal- 
culated from the Nelkin water kernel, Eq. (5.24) and Fig. 
(5.2), for the water moderator. The water transport cross 
section derived from the Nelkin model is also used in the 
calculation of the moderator escape probabilities with use 
of the monatomic gas kernel. For all other materials ]I(E) 
is given by 

U(E) - H = 2/3 A , (7.3) 

A being the mass number of the scattering material. 

Calculation of the relative fuel absorption f(E)/ 

[ 1-f (E) ] , Eq. (2.29), using the P g (E) shown in Fig. (3.1), 

72 / 

— 'D. Hughes and B. Schwartz, "Neutron Cross Sections," 

BNL 325, 2nd Ed., Brookhaven National Laboratory (1958) 



Sec. 7.1 



78 



the P m (E) shown in Fig. (3.2), and the P e (E) shown in Fig. 
(4.2), yields the apparent lattice absorption cross section. 



Z™(E) 



1 + f(E)/[l-f(E)] 



This quantity is shown in 



Fig. 7.1 for the five Cornell ZPR cores. Comparison with 
the absorption cross section of pure water, Z m (E), shows 

a 

the non 1/v behavior of this apparent cross section as well 

235 

as the effect of the 0.273 eV resonances in U .It is 
this non 1/v behavior which yields a temperature coefficient 
of the thermal utilization, f. 

The apparent lattice absorption cross section shown 
in Fig. (7.1) is substituted into the neutron balance 
equation, Eq . (6.1), along with either the monatomic gas 
or Nelkin water scattering cross section. The computer 
solution of Eq . (6.1) is discussed in Appendix C, the result 
of the solution being the spatially averaged moderator 
spectrum, cp^E), evaluated at the discrete energy points 
in the energy mesh used. Figure 7.2 shows the results of 
this calculation utilizing the Nelkin water kernel for the 
five water to uranium oxide volume ratios of the Cornell 
ZPR. Since all of the results are normalized to one neutron 
per second slowing down per unit volume of the moderator. 
Fig. 7.2 shows the flux suppression as well as the spectral 
hardening experienced in the tighter (lower water to oxide 
volume ratios) cores. The difference between moderator 
spectra predicted with the monatomic gas model and the 
Nelkin water kernel is shown in Fig. 7.3 where the y“(E) 
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Figure 71 Apparent Lattice Absorption Cross Section 
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Figure 7.2 Spatially Averaged Moderator Thermal Spectra 
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Figure 7.3 Comparison of the Moderator Spectra Predicted with the 
Nelkin Water and Monatomic Gas Models 
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predicted by these two models are plotted for the extreme 

lattices of the ZPR; that is the 1:1 and 4:1 water to oxide 

volume ratio cores. Figure 7.3 indicates that harder 

spectra are predicted with the use of the Nelkin water 

kernel than with the use of the monatomic gas kernel. 

Using the spatially averaged moderator flux, the 

spatially averaged fluxes in the fuel, cp (E), and the 

cladding, cp c (E), may be calculated using Eqs . (6.3) and 

(6.4) respectively. Results typical of this calculation 

using the Nelkin water kernel are shown in Fig. 7.4 where 

cp^E), qj c (E), and q>^(E) are plotted for the ZPR 3:1 core. 

Spectrum hardening as well as depression, giving rise to 

disadvantage factors, is experienced as one goes from the 

moderator to the cladding to the fuel. 

The average neutron velocities in the moderator, 

cladding, and fuel as defined in Eq . (6.6) are shown in 

Fig. 7.5 for all of the ZPR cores, these calculations 

being made with the Nelkin water kernel. Again the 

spectral hardening produced in the tighter cores is 

evident. Figure 7.5 also shows that for the tighter 

cores v c is well represented by the relationship which 
73 / 

Honeck — has verified for the Brookhaven uranium metal 
cores, 

— ^H. C. Honeck, "The Calculation of the Thermal Utiliza- 
tion and Disadvantage Factor in Uranium/Water 
Lattices," Nucl. Sci. Eng. JL8, 49-68 (1964). 
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cp(E), n/cm -sec-eV 
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Figure 7.4 Spatially Averaged Cell Spectra 
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Figure 7.5 Average Neutron Velocities 
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v c “ 2^ v f + V m^ ‘ ^ 7 ‘ 4 ^ 

The error in the result predicted by Eq . (7.4) increases 
with the water to oxide volume ratio to a maximum of 1.8% 
for the 4:1 core. 

The average neutron velocities , v^, for the homogen- 
ized cell, Eq. (6.33), are shown for calculations with 
both the monatomic gas model and the Nelkin water kernel 
in Fig. 7.6. Again the harger spectra predicted by the 
Nelkin water kernel is evident as the average neutron 
velocities for the homogenized cell calculated with the 
monatomic gas model are lower than those given by the 
Nelkin model by 5.4% for the 1:1 core and 2.4% for the 
4 : 1 core . 

The neutron density disadvantage factors for the 
moderator, 5™, and for the cladding, 5^, as defined in 
Eqs . (6.8) and (6.9) respectively, are shown in Fig. 7.7 
as calculated with both the monatomic gas kernel and the 
Nelkin water model . In this case, the discrepancies 
between the results computed with the two models is 
entirely attributed to the softer spectra produced with 
the monatomic gas model since the transport cross section 
for water used in both calculations was that calculated 
with the Nelkin water kernel. The 5™ predicted with the 
monatomic gas kernel are 1.4% higher for the 1:1 core and 
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figure 7.6 Average Neutron Velocity for the Homogenized 

Cell 



Neutron Density Disadvantage Factors 
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Figure 7 Neutron Density Disadvantage Factors 
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0.8% higher for the 4:1 core than those calculated with the 
Nelkin model. Similar ly, the differences in the cal- 
culated by the two models is 0.8% for the 1:1 core and 
0.3%, for the 4:1 core. The fact that the disadvantage 
factors monotonically increase as a function of water to 
oxide volume ratio is expected— ^ since a non-reflecting 
cylindrical cell boundary condition is incorporated in this 
formulation . 

The results of the calculation of the thermal utili- 
zation, f, and rjf defined by Eqs . (6.15) and (6.16) 
respectively for the Nelkin water model are shown in 
Figs. 7.8 and 7.9. The results predicted with the mon- 
atomic gas model closely approximate those computed with 
the Nelkin water kernel, the differences being less than 
0.2% for all cases of f, r), and rjf. This indicates that 
these integral properties are essentially insensitive to 
the scattering model used in the averaging process . 

Figure 7.10 shows the variation of thermal lifetimes, 
Eq. (6.17), as a function of water to oxide volume ratio 
calculated with the Nelkin water kernel and the monatomic 
gas model. Again relatively small variations, always less 
than 0.4%, between the Nelkin water and monatomic gas kernel 
calculations are present. The lifetimes are seen to decrease 
in the tighter cores because of the increased apparent 
absorption. Fig. 7.1, present in these cores. 



74 / 



H. C. Honeck, p. 56, see Footnote 73. 
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Figure 7.8 Thermal Utilization 
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Figure 7.9 Integral Properties, r\ and 
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Figure 7.10 Thermal Lifetime 
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Figure 7.11 shows the effective absorption cross 
section, 2 , Eq . (6.31), and the mean absorption cross 

a 

™h 

section, 2", Eq . (6.32), for the homogenized cell as 

a 

calculated using the Nelkin water and monatomic gas 
models. Variations encountered in the calculation of 
these results when the monatomic gas model was utilized 
were less than 0.4% for the effective absorption cross 
sections and 2.0% for the mean absorption cross sections. 

The thermal diffusion coefficient, D, and the 
thermal diffusion length, L, as calculated by averaging 
a homogenized mean free path in accordance with Eq . (6.38) 
are shown in Fig. 7.12 using both the Nelkin water and 
monatomic gas kernels. From Fig. 7.12 it is seen that 
the monatomic gas model predicts a D which is about 3.0% 
less, and an L which is about 4.0%, less than those cal- 
culated with the Nelkin water kernel . The differences 
between these results are strictly due to the relative 
hardness of the spectra predicted by the two models, since 
the transport cross sections utilized were the same in 
both cases . It is still not clear how these calculated 
diffusion lengths would compare with those measured in a 
buckling experiment; however, it is noteworthy that the 

minimum in L for the 2:1 core is the same phenomenon 

75 / 

measured for the Brookhaven uranium metal lattices . — 

— ^ R. L. Hellens and H. C. Honeck, "Summary and Preliminary 
Analysis of the BNL Slightly Enriched Uranium, Water 
Moderated Lattice Measurements," Light Water Lattices, 
p. 54, International Atomic Energy Agency, Vienna (1962). 
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Figure 7.11 Effective and Mean Absorption Cross Sections for 

Homogenized Cell 
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Figure 7.12 Thermal Diffusion Coefficient and Length for the 

Homogenized Cell 
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This also indicates that the customary receipe for cal- 

culating a homogenized diffusion area, L = (1-f) L — 

lacks validity in that it predicts a monotonically 
2 

increasing L with increasing water to oxide volume 
ratios . 

A tabular listing of the parameters presented in this 
section appear in Appendix B as calculated with the Nelkin 
water kernel at 20°C. A discussion of the accuracy of 
these results is given in Chap. 8. 



7.2 Temperature Coefficients 

In order to estimate the effect of temperature on the 
various integral parameters defined in Chap. 6 and pre- 
sented for 20°C in Sec. 7.1, all of the calculations were 
remade at a physical temperature of 40°C. In these cal- 
culations any changes in the fuel and cladding densities 
and dimensions have been neglected. In addition, the 
expansion of grid plate supporting the fuel elements leading 
to a change in the equivalent cell radius has also been 
neglected. All of the results are expressed in terms of 
temperature coefficients defined for a given variable, x, 
as 



1 Ax 1 

x AT " x (20 u ) 



x (40° ) - x(20° ) 0 „-l 

— — c 



(7.5) 



— ^ A. M. Weinberg and E. P. Wigner, The Physical Theory of 
Neutron Chain Reactor s , p. 724, The University of 
Chicago Press, Chicago, 1958. 
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The temperature coefficient of the average neutron 
velocity in the homogenized cell, v^, Eq . (6.33), is shown 
in Fig. 7.13 for both models of the scattering kernel. 

This indicates that although the Nelkin model predicts a 
harder homogenized spectrum at room temperature, the 
monatomic gas spectra hardens at a faster rate with 
increasing temperature. It is therefore possible to 
postulate that at some elevated temperature the integral 
parameters dependent primarily on the degree of spectral 
hardening will be the same for both models of the scattering 
kernel. The temperature coefficient of v^ predicted with 
the monatomic gas model is 23.8% higher for the 1:1 core 
and 7.8% higher for the 4:1 core than those predicted by 
the Nelkin water kernel. 

Figure 7.14 shows the temperature coefficients of 
the moderator and cladding disadvantage factors, 6™ and 
5 , these parameters having been defined in Eqs . (6.8) and 
(6.9) respectively. As expected, they are negative; the 
effects of harder spectra being seen in Fig. 7.7. The 
monatomic gas model predicts larger negative temperature 
coefficients because of the greater relative spectral 
hardening encountered with that model . The maximum differ- 
ence between the scattering models occurs in the 1:1 core 
where the temperature coefficient of 6™ predicted by the 
monatomic gas model is 30% higher than the prediction by 
the Nelkin water kernel. In the case of 6^ the results 
differ by 20%. 
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Figure 7.13 Temperature Coefficient of the Average Neutron 
Velocity for the Homogenized Cell 
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Figure 7.14 Temperature Coefficients of the Neutron Density 

Disadvantage Factors 
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The variation of the temperature coefficients of the 
thermal utilization, f, Eq . (6.15), and r)f, Eq . (6.16) is 
shown in Fig. 7.15, the difference between the two curves 
for each scattering model being equal to the temperature 
coefficient of r\. Again the results predicted are more 
positive with the use of the monatomic gas model because of 
the greater relative spectral hardening given by that model . 
It is to be noted that the temperature coefficient of tj f 
is slightly negative for the 1:1 core. 

The final temperature coefficients worthy of note are 
those of the diffusion coefficient, D, Eq . (6.38), and the 
thermal diffusion length, L, Eq. (6.39). These are shown 
in Fig. 7.16; the difference between the monatomic gas and 
the Nelkin water models being quite evident . Again, the 
temperature coefficients predicted by the monatomic gas 
model are higher than those produced by the Nelkin water 
kernel. The deviation is the most pronounced for the 1:1 
core where the monatomic gas results are higher than the 
Nelkin water kernal results by the following percentages : 
temperature coefficient of D, 68%; and temperature co- 
efficient of L, 34%. 

The differences between the monatomic gas and the 
Nelkin water kernel results are more pronounced for the 
D and L temperature coefficient calculations than for any 
of the other parameters previously discussed. This may 
be attributed to the fact that these calculations are 
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Figure 7.15 Temperature Coefficients of the Integral Prop 

erties, f and rjf 
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Figura 7,16 Temperature Coefficients of the Thermal Diffusion 
Coefficient and Length for the Homogenized Cell 
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doubly sensitive to the degree of spectral hardening pre- 
dicted by the two models * since the homogenized cell trans- 

h 

port cross section, 2“ r (E), Eq . (6.36), used in these cal- 
culations is derived by flux and volume weighting the 
transport cross sections of the various regions in the 
cell. The reciprocal of this homogenized transport cross 
section is then averaged over a homogenized spectrum, 

Eq. (6.37), making the resultant D highly sensitive to 
the degree of hardness in the spectra utilized in these 
averaging processes . 

An estimate of the accuracy of the temperature 
coefficient calculations is discussed in Chap. 8, while 
a tabular listing of all of the temperature coefficients 
calculated with the Nelkin water kernel is given in 
Appendix B . 



CHAPTER 8 



ACCURACY OF THE RESULTS 

In an attempt to estimate the accuracy of the results 

predicted with program COUTH, as presented in Chap. 7 and 

listed in Appendix B, two lines of approach have been 

utilized. The first of these involved the calculation by 

COUTH of several integral properties of the Brookhaven 

uranium dioxide lattices, a comparison then being made between 

77 / 

these results and those calculated by Honeck — with the 
THERMOS transport theory code. The second of these was an 
investigation of the sensitivity of the Cornell ZPR results 
to changes in the inputs to the COUTH code such as the den- 
sities, cross sections, and dimensions of various regions 
within the lattice cell. 

78 / 

The Brookhaven cores—' studied are generically related 
to the Cornell ZPR cores in that they are light water mod- 
erated and utilize low enrichment uranium dioxide fuel. 

Five water to uranium oxide volume ratios are provided, 
these being 1.291:1, 1.600:1, 2.050:1, 2.803:1, and 4.000:1. 
The fuel rods are 0 .444 inches in diameter and clad with 
0.028 inch thick stainless steel giving the fuel element 
an outside diameter of 0.500 inches. The uranium dioxide 
fuel is enriched to 3.0 atom per cent, its density being 
9 .3 gm/cc . 

H. C. Honeck, "The Calculation of the Thermal Utilization 
and Disadvantage Factor in Uranium/Water Lattices," 

Nucl. Sci. Eng. 18^, 49-68 (1964). 
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The integral thermal parameters compared are the thermal 
utilization, f, the ratio of the number of fission neutrons 
produced to the number of thermal neutrons absorbed in the 
fuel, 7], and the moderator neutron density disadvantage 
factor, 8™. The THERMOS calculation of these parameters 
utilizes a thirty point energy mesh ranging from 0.00025 to 
0.632 eV, while COUTH uses a 41 point energy mesh ranging 
from 0.001 to 0.645 eV, as discussed in Sec. 7.1. Both 
programs utilize the same cross section input data and both 
perform their calculations utilizing a Nelkin water model 
for the scattering kernel at 20°C. 

A comparison between the results predicted by THERMOS, 
using integral transport theory, and COUTH, using the 
simplified cell theory, are shown in Fig. 8.1 for the thermal 
utilization, f. In Fig. 8.1 it is seen that the f's cal- 
culated with the two codes compare very well, the maximum 
deviation in the results occurring in the 4.0:1 core where 
the THERMOS f is 0.39% higher than the f calculated with 
COUTH . 

The results of the calculation of 7) by COUTH and 
THERMOS are seen in Fig. 8.2 to be even closer than those 
for f. In this cas'e the maximum difference between the 
results predicted by the two techniques is found in the 
1:1 core, the THERMOS result here being 0.19% higher than 
the COUTH estimate. 



78 / 



H. C. Honeck, see Footnote 77. 
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Water/Oxide Volume Ratio 
Figure 8.1 Thermal Utilization 
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In the case of the moderator disadvantage factor* o™* 
the agreement between the two codes is not quite as close 
as in the calculations of f and q. Figure 8.3 shows the 
results derived* indicating that THERMOS predicts higher 
disadvantage factors for the tighter cores than COUTH* while 
this role is reversed with the larger wafer to oxide volume 
ratios. The THERMOS 6™ for the 1.29:1 core is 0.62% higher 
than the COUTH number* while for the 4.0:1 core the COUTH 
6™ is 0.97% higher than that given by THERMOS. Honeck— ^ 
assigns a combined uncertainty to the THERMOS calculation 
of the disadvantage factor of 1.5%* meaning the COUTH cal- 
culation of this parameter has an uncertainty less than 
2.5%. The deviations between COUTH and THERMOS in this 
calculation are believed to be attributable to COUTH's use 
of simple diffusion theory in the moderator region* com- 
bined with an overestimation of the fuel element escape 
probability. Resolution of these differences requires 
further study and is discussed in Chap. 9. 

The parametric study made on the Cornell ZPR for 20°C 
using the Nelkin water kernel investigated the sensitivity 
of the calculated parameters to changes in the input data . 
Changes where made in the fuel density* the fuel enrichment* 
the fuel rod radius* the cladding radius* the cladding 
scattering and absorption cross sections* the moderator 
density* and lattice spacing or pitch. In an attempt to 
79 / 

— H. C. Honeck* p. 60* see Footnote 77. 
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Water/Oxide Volume Ratio 

Figure 8,3 Neutron Density Moderator Disadvantage Factor 
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determine the effect of using various numerical integration 
schemes, the results utilizing an integration method which 
fit exponentials, i.e., straight lines on log paper, were 
also studied. The effect of these changes were noted for 
several parameters, these being: The thermal utilization, 

f; the ratio of the number of fission neutrons produced to 
the number of neutrons absorbed thermally in the cell, qf; 
the moderator neutron density disadvantage factor, 6™; the 
thermal diffusion length, L; the thermal lifetime the 

mean neutron velocity for the homogenized cell, v^; and the 
mean total cross section for the homogenized cell, The 

results of the parameter variation study are shown in 
Table 8.1 where the percentage changes in the input para- 
meters required to produce a positive 0.1% change in the 
resultant integral properties are listed. The results 
presented in Table 8.1 represent the most sensitive cases, 
the error analysis having been carried out on the 1:1 and 
4 : 1 ZPR cores . Thus the percentage changes shown which 
involve changes in the fuel element parameters come from 
the 1:1 core analysis, since a relatively large portion of 
the unit cell is composed of fuel for this core. The 1:1 
core is also most sensitive to changes in lattice pitch. 

On the other hand results from the 4:1 core are given for 
the sensitivity of parameters to changes in the water 
density, water playing a relatively more important role 
in this core. 



PERCENTAGE CHANGES IN INPUT PARAMETERS REQUIRED TO CHANGE INTEGRAL PROPERTIES 

BY +0.1 PERCENT 
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The use of the numerical integration scheme involving 
the fitting of exponentials rather than straight lines as 
in the trapezoidal integration rule produces the following 
maximum percentage changes in the integral parameters: 

f, 0.004%; 7jf, .009%; 6™, 0.12%; L, 1.1%; .069%; 

— 

v^, 0.127 o ; and Z^;, 2.5%. The large percentage change in 
““h 

z“ may attribute to the fact that the numerical integration 
scheme using exponentials predicts scattering cross sections 
for water by integrating the scattering kernel which are as 
much as 3.0% higher than the experimentally measured 
values . 

Table 8.1 shows that the numerical results are 

essentially insensitive to the cladding cross section 

data,, while they are most sensitive to the lattice pitch. 

80 / 

The actual tolerances given—' for the Cornell ZPR lattice 
pitches are less than 0.10% so that the output parameters 
will vary no more than about 0.2% because of variations 
in this dimension. The fuel element parameter, density, 
enrichment, fuel rod radius, and fuel element radius are 
known to within 0.5% and reference to Table 8.1 indicates 
that in general variations of the thermal integral prop- 
erties will be less than 0.5%. 

In light of the results of comparing the predictions 
of the THERMOS and COUTH codes and the results of the 
input parametric studies, the following probable uncer- 
tainties may be assigned to the 20°C calculations performed 
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by COUTH utilizing the Nelkin water kernel as a scattering 
model . 



1 . 


f, 7), and 7)f; 0.5%. 






2. 


Disadvantage factors for 


the moderator 


and 




cladding, 2.0%. 






3. 


Thermal lifetimes; 3.0%. 






4. 


Mean velocities; 3.0%. 






5. 


Averaged cross sections. 


including the 


calculation 



of D and L; 57 0 . 

With regard to the room temperature calculations made 



utilizing the monatomic gas model, the deviations from the 
Nelkin water kernel results were quoted in Chap. 7. Using 
these percentages, one can then assign the following prob- 
able uncertainties to the 20°C monatomic gas calculations . 

1. f, r), and 7)f; 0.8%. 

2. Moderator and cladding disadvantage factors; 4.0%. 

3. Thermal lifetimes; 4.0%. 

4. Mean velocities; 8.0%. 

5. Averaged cross sections, including the diffusion 

coefficient and length; 10.0%. 

On the basis of these two sets of uncertainties, one 
must reconcile the accuracy obtained with the computer time 
required for a given calculation. A monatomic gas scattering 
kernel, in the form of a 45 by 45 matrix takes only two and 
one half minutes of computer calculation, while the Nelkin 
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water scattering kernel takes on the order of fifteen 
minutes . Thus if one is undertaking a detailed tempera- 
ture study wherein multiple scattering kernels must be 
calculated, use of the monatomic gas model is advocated 
for one primarily interested in f and rjf. However, if 
emphasis is to be placed on the other integral parameters, 
and in particular the averaged cross sections, then the 
more complicated Nelkin model should be used. 

The assignment of uncertainties to the temperature 
coefficients presented in Chap. 7 presents a slightly 
different problem. In this case one must estimate the 
uncertainties in a calculation at a given temperature and 
then estimate the uncertainty in the degree of spectral 
hardening predicted by a certain model . The fact that 
the THERMOS and COUTH results at room temperature are 
in agreement for all values of water to fuel volume ratios 
indicates that COUTH can adequately handle spectral harden- 
ing effects if the scattering kernel model can predict the 
appropriate temperature changes in the energy transfer 
process . We therefore, apriori, assign uncertainties to 
the temperature coefficient calculations as being twice 
those stated above for the room temperature calculations 
with both the Nelkin water model and the monatomic gas 
kernel . 

80 / 

— S. S. Berg, "Initial Experiments on the Cornell University 
Zero Power Reactor Cores," Cornell University Th 
(1964) . 



CHAPTER 9 



SUMMARY AND CONCLUSIONS 

A simplified polyenergetic cell theory has been formu- 
lated for use in the determination of the spatially averaged 
energy dependent thermal fluxes in the various regions of a 
lattice cell. The derived spectra are then utilized in the 
calculations of the thermal integral parameters and average 
cross sections required for reactor computations . 

The cell theory, as formulated, postulates an infinite 
moderator region with the absorption cross section of this 
region appropriately modified to account for the neutron 
leakage into and absorption by the fuel element in the 
actual lattice. The modifications to the moderator absorption 
cross section, have been formulated both in general terms 
an in terms of moderator and fuel element escape probabilities, 
the latter approach offering physical transparency and ease 
of calculation. 

Analytic expressions for the escape probabilities have 
been presented, integral transport theory having been applied 
in the fuel element region, while diffusion theory was utili- 
zed in the moderator region. Using these analytic expressions 
the theory was applied to actual lattices in the form of the 
Cornell University Zero Power Reactor cores to determine room 
temperature parameters and their temperature coefficients . 
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In an attempt to estimate the accuracy of these results, 
the method was also applied to the Brookhaven uranium 
dioxide cores and the results then compared with those 
predicted by multigroup transport methods . The remarkably 
good comparison between these two sets of calculations in 
conjunction with a parametric study which yielded a relative 
low sensitivity of the calculated parameters to changes in 
the input data has led to the assignment of rather small 
uncertainties in the results calculated with the simplified 
cell theory. 

Although the estimated uncertainties tend to instill 
confidence in the calculational results, one must insure 
that these events are not just fortuitous and that the 
simplified cell theory is in fact predicting behaviors 
which closely approximate the physical processes taking 
place . 

In this regard, there are three areas in which further 
work is necessary in the simplified cell theory as pre- 
sently formulated. The first of these concerns the validity 
of the assumptions made in the calculations of the modera- 
tor escape probabilities, or alternatively, in the calcu- 
lation of gradient of the moderator flux at the fuel 
element surface. In these calculations it was assumed that 
the shape of the flux in the moderator was given by mono- 
energetic diffusion theory applied at each energy point. 
Although multigroup transport calculations have given some 
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support to this assumption,, its accuracy can certainly be 
questioned particularly in the tighter cores where fuel 
element spacing becomes comparable with a mean free path. 
Thus there is a need to accurately determine the flux 
gradient, or more precisely, the net current as a function 
of energy at the fuel element moderator interface. Any 
study such as this should include the effects of a non 
uniform spatial source of slowing down neutrons. The use 
of multigroup transport theory specif icially oriented 
toward this problem is suggested for this future study. 

The second area of concern involves the calculation 
of the fuel rod escape probability. In the expression 
utilized in these calculations an isotropic flux was 
assumed incident upon the fuel rod . This in turn leads 
to escape probabilities which are too large; the flux in 
an actual lattice is more peaked toward the rod. The use 
of smaller escape probabilities for the fuel rod would 
lead to larger values of f and qf, along with larger 
values for the calculated disadvantage factors . Thus 
there is need for at least a parametric study utilizing 
all of the analytic approximations for the fuel rod 
escape probability presented. One could then turn his 
attention to variations caused with changes in angular 
distribution of the neutrons incident on the fuel 



element . 
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The final area where further study is needed concerns 
the calculation of temperature coefficients . The immediate 
need is to compare temperature coefficients calculated with 
integral transport theory with those predicted by the 
simplified cell theory as was done with the room temperature 
case. However, if both calculations utilize the same 
scattering model, it is believed the comparison of the 
temperature coefficient computations will be comparable 
with those for room temperature, since the simplified cell 
theory has demonstrated its ability to handle spectral 
hardening in five different cores at 20°C. The accuracy 
of the temperature coefficient results thus reduces to a 
question of how well the scattering model predicts spectral 
hardening as a function of temperature. This question is 
perhaps best answered by the comparison of pure moderator 
spectra predicted by the model at various temperatures 
with those measured experimentally. 



APPENDIX A 



LATTICE PARAMETERS FOR THE CORNELL UNIVERSITY ZERO POWER 

REACTOR 

This section details the input parameters which were 
utilized by program COUTH in the calculations described 
in Chap. 7 for the Cornell University Zero Power Reactor 
cores . Appendix D deals with the Fortran FORMATS used in 
COUTH for reading data into the program. 

Table A.l lists the fuel element parameters for the 
Cornell ZPR while Table A. 2 gives the dimensional param- 
eters of the five hexagonal lattices of the ZPR. Material 
concentrations and cross sections are given in Table A. 3. 
The absorption cross sections, cr , listed are evaluated 
at E 0 = .0253 eV while the scattering cross sections given 

are the free atom values, . The number of neutrons 

235 

produced per thermal fission for U , v, is 2.43 while 

235 

the value of the U fission cross section at E Q is 
81 / 

580.9 barns. — 7 Appendix B gives a complete listing as 

a function of energy for the scattering and transport 

cross sections of water, the absorption, fission, and 

235 

scattering cross sections of U , and the scattering 

238 

cross sections of U and oxygen. 



— - H. C. Honeck, "The Calculation of Thermal Utilization 
and Disadvantage Factor in Uranium/Water Lattices," 
Nucl . Sci. Eng. 18 , 66 (1964). 
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TABLE A . 1 

CORNELL UNIVERSITY ZERO POWER REACTOR FUEL ELEMENT 



PARAMETERS 




Fuel Material 


Uranium Dioxide 


Fuel Rod Diameter, in: 


0 .600±0 .003 


Fuel Density, gm/cc: 


10 .381±0 .52% 


Fuel Enrichment, Weight Percent U 


2 .070±0 . 10 


oot; 

Fuel Enrichment, Atom Percent U : 


2 .096±0 .90 


Oxygen to Uranium Number Ratio: 


2.01:1 


Fuel Molecular Weight: 


270 . 26 (physical) 


Cladding Material: 


6061-76 aluminum 


Cladding Density, gm/cc: 


2.70 


Cladding Inner Diameter, in: 


0 .610±0 .002 


Cladding Thickness, in: 


0 .028±0 .003 


Fuel Element Diameter, in: 


0 . 666± . 006 



a/ 

— S. S. Berg, "Initial Experiments on the Cornell University 
Zero Power Reactor Cores," Cornell University Thesis 
(1964). 
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D. Clark, "The Cornell University Nuclear Reactor Laboratory," CURL- 
Cornell University (1961) . 
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TABLE A. 3 

CORNELL UNIVERSITY ZERO POWER REACTOR MATERIAL CON- 
CENTRATIONS AND CROSS SECTIONS 



Material 


Concentration 
(atoms /barn-cm) 


(barns ) 


(barns 


u 235 


4.8498xl0" 4 


679 .0 


10.0 


y238 


2.2654xl0" 2 


2.71 


8.3 


Oxygen in UO 2 


4.6510xl0" 2 


0.0 


3.76 


Aluminum 


6.0205xl0“ 2 


0.230 


1.40 


Hydrogen in ^0 




0.328 


20.4 


20°C 


6.6746x10 2 






40°C 


6.6346xl0~ 2 






Oxygen in ^0 




0.0 


3.76 


20°C 


3 .3373xl0“ 2 






40°C 


3.3173xl0" 2 







a/ 

— H. C. Honeck, "The Calculation of the Thermal Utilization 
and Disadvantage Factor in Uranium/Water Lattices/' 
Nucl. Sci. Eng. 18 , 66 (1964). 



APPENDIX B 



NUMERICAL RESULTS FOR THE CORNELL UNIVERSITY ZERO 

POWER REACTOR 

Chapter 7 discussed some of the results of the cal- 
culations made on the Cornell University Zero Power Reactor 
using the simplified cell theory as incorporated in COUTH. 
This section presents all of the results in a tabular form 
derived from the calculations utilizing the Nelkin water 
model for the scattering kernel. 

Table B.l lists the energy , velocity and lethargy 
mesh used in the calculations by COUTH. The derivation 
of this mesh is discussed in Sec. 7.1, and all of the cal- 
culations are carried out using the lethargy mesh. 

Table B.2 gives the scattering kernel in matrix form 
as derived from the Nelkin water model, Eq. (5.15), for 
a physical moderator temperature of 20°C. The scattering 
matrix is expressed in barns /unit lethargy and it is 
related to the scattering matrix in barns /eV by, 

tf(uj— ^u i ) = E ± a (E j — 9 “E ;L ) . (B.l) 

The subscripted indices i and j correspond to the lethargy 
and energy points listed in Table B.l. The diagonal terms 
of the scattering matrix listed in Table B.2 are those cal- 
culated with Eq . (5.22) and they are not modified in 
accordance with Eq . (5.26). 
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TABLE B.l 

ENERGY, VELOCITY, AND LETHARGY MESH UTILIZED BY COUTH 



I 


FNEPfiY 


DFLT A F 


VFLOrlTY 


DELTA V 


1 FTHaRGY 


DELTA u 
















1 


1 • 6o9E 000 


3.782E-001 


17547 . 0 


2200 . 0 


. 0 000 


.2679 


2 


1. 23l£. QQQ 


2 . 644E-001 


15347.0 


1747.5 


.2679 


.2416 


3 


9 . 6*7E-Q01 


1 . 873E-0 01 


13599.5 


1388.1 


.5097 


.2153 


4 


7 * 7 o 4 F - o n 1 


1 . 344F-0 nl 


1 221 1 . 4 


1102.6 


.7250. 


.1893 


5 


6.4«inE-noi 


9.770E-002 


11108.7 


875.8 


.9143 


.1642 


6 


5 • 4 7 3 E- q 0 1 


7.189F-002 


10232.9 


605.7 


1.0785 


.1408 


7 


4 • 7«i4E"001 


5.350E-002 


9537.2 


552-6 


1.2194 


• 1194 


8 


4 . 24- 9E~ H 0 1 


4 . 0 ??P * n 02 


89.84^6- 


439 . Q 


1 .3387 


• 1 002 


9 


3 . 81 7E-001 


3. 051E-002 


8545.6 


348.7 


1.4389 


. 0833 


10 


3.5i?F-f)0l 


2.333F-002 


8197.0 


277.0 


1 .5292 


. 0687 


11 


3 • 279E-n 0 1 


1 . 796E-0 02 


7920.0 


220-0 


1.5910 


* 0563 


12 


3 • 0 09E* n 0 1 


1 . 746E-0 02 


7700 . 0 


220 . 0 


1.6473 


. 0580 


13 


2.9P4E-001 


1.695E-002 


7480.0 


220.0 


1.7053 


. 0597 


14_ 


2 • 755E~o01 


1 .644E-002 


7260*0 


220 * 0 


1 .7650 


.0615 


15 


2*5ooE-OOl 


1 .594E-002 


7040.0 


220.0 


1.82*5 


. 0635 


1 6 


2 , 431 E-nni 


1 . 543F-0n2 


6820 . 0 


220 . 0 


1 .89o0 


. 0656 


17 


2.277E-001 


1 .493E-002 


6600 . 0 


220.0 


1.9556 


. 0678 


18 


9.1?flE-onl 


1 . 442F-0n2 


6380 . 0 


290 . 0 


9.0234 


. 0702 


19 


1 . 9B3E-001 


1.391E-002 


6160 . 0 


220.0 


2.0936 


. 0727 


2Q 


1 . Ba4F-nnl 


1 . 341E-0 0? 


5940 . 0 


220 • 0 


2-16*3 


. 0755 


21 


1 • 7i nE-noi 


1.290E-OO2 


5720.0 


220.0 


2.24i8 


. 0784 


22 


1 . 5ai F-nQl 


1 . 240F-0n2 


*5 ‘S 0 0 , 0 


2j>0 - 0 


2 T 3 2 0 3 


. 081 * 


23 


1 . 4«?7E-n0l 


1.189E-002 


5280 • 0 


220.0 


2.4019 


. 0851 


24 


1 . 3 3 8 E - 001 


1 , 138F-0Q? 


5 n * 0 - 0 


220 . 0 


2 .4870 


. 0889 


25 


1 •2?4E-n0l 


1.088E-002 


4840 . 0 


220 . 0 


2.5759 


. 0930 


26 


1 . li (SE- o-0 1 


1 . 037E-0 02 


462p ..0 


-22Qj-Q 


2 . 660 Q 


jlXI9_76 


27 


1 • Oi ?E-noi 


9 . 866E- 0 0 3 


4400.0 


220.0 


2.76*6 


.1026 


?8 


9.13?E-o02 


9.36nE-0n3 


418 0 . 0 


290 . 0 


2.8691 


.1081 


29 


8.1QRE-002 


8.854E-003 


396o . 0 


220 . 0 


2.9773 


.1143 


30 


Z* 3ilE~o02 


fl.348E-0f)3 


374Q.0- 


220--4- 


_ 3 ^0 9 1 6 


.1219 


31 


6.476E- n o2 


7.842E-003 


3520 . 0 


220 *0 


3.2128 


• 1291 


32_ 


5 . 609 E - 0 Q 2 


7 . 336F- 0 p 3 


33 0 Q . Q 


290 . 0 


3.34-1 9 


# 1. 3 8 0 


33 


4 . 9*flE-o02 


6.830E-003 


3080.0 


220 . 0 


3.4709 


.1482 


3 4 


4 . ?7=;F-no? 


A , 394F-Qo3 


28*0-0 


?20 ■ 0 


3,6281 


*1601 


35 


3 • 64 3E- n 0 2 


5.818E-003 


2640 . 0 


220 . 0 


3.7882 


.1740 


36 


3 . 0 M F-r) 0 ? 


5 , 3l3F-no3 


? 4 2 Q . Q 


220*0 


3.9622 


.194)6 


37 


2.530E-002 


4 . 8o7F-On3 


2200.0 


220 *0 


4.l5p9 


.2107 


38 


2* n49E-QQP 


4 .3(11 F- 0 p3 


1 98 0 . Q 


220_. 0 


4.3636 


.2356 


39 


1 . 61 9F- n 02 


3.795E-003 


176o . 0 


220-0 


4.5991 


.2671 


4f| 


1.24nE-n02 


3.289E-0n3 


1540.0 


220.0 


4.86*2 


.3083 


41 


9* 107E-0O3 


2.783E-003 


1320.0 


220 . 0 


5.1745 


.3646 


42 


6 • 3 p4 F - Q 0 3 


2.277E-0O3 


1100 . 0 


220.0 


5.5391 


.4463 


43 


4.048E~o03 


1 .771E-003 


880 . 0 


220 • 0 


5.9854 


.5754 


44 


2.277F-003 


1 . 265F-003 


660 . 0 


220.0 


6 . 56 n 8 


.8109 


45 


l.Ol 2F-003 


. 0 0 OE 0f)0 


44Q.0 


. 0 


7.3717 


. 00 0 0 
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Table B.3 gives the Nelkin water kernel scattering 

matrix for a moderator physical temperature of 40°C. Its 

format is the same as that described for Table B.2. 

Table B.4 lists as a function of energy the scattering 

cross section, Eq . (2.2), and the transport cross section, 

Eq. (5.25), for water. The energy dependent scattering 

23 5 238 

cross sections for U , U , and oxygen are also listed 
in Table B.4 as calculated using Eq . (5.13) with the param- 
eters listed in Table A. 3. 

Table B.5 gives as a function of energy the absorption 
and fission cross sections— ^ for as well as the 

uranium dioxide fuel fission and absorption cross sections. 
One plus the capture to fission ratio, [l+a(E)], for the 
fuel is also listed in Table B.5, along with the normalized 
slowing down source of neutrons, S(E), Eq. (7.1 ), used in 
the calculations . 

Table B.6 through B.10 list the energy dependent 
thermal spectra calculated by COUTH for a moderator 
physical temperature of 20°C. The moderator, Eq. (6.1), 
fuel, Eq. (6.3), and cladding, Eq . (6.4), spectra are 
given for the five Cornell ZPR cores . 

Table B.ll gives the integral properties defined in 
Chap. 6 for the five ZPR cores at a temperature of 20°C, 

82 / 

— ' J. Suich, "Temperature Coefficients in Heterogeneous 
Reactor Lattices," Massachusetts Institute of 
Technology Thesis (1963) . 
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TABLE B .6 

SPATIALLY AVERAGED THERMAL SPECTRA FOR THE CORNELL 

ZPR 1:1 CORE 
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moderator 


CLADD i NG 


FUEL 






spectrum 


spectrum 


spectrum 
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j .6093E 00 


5.4517E-01 
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5.3873E-01 
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1.2311E 00 


6 . 1583E-01 


6.1323F-01 


6 . 0196E-01 
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9.6687E-01 


8 . 0549E-01 
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1 . 7 8 1 4 E 0 0 
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10 
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12 
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2.5621E 00 


2 . 490 OF 00 


2.3562E 00 


13 


P.9244E-01 


2.6780E 00 


2.5986E 00 


2.4525E 00 


14 


2.7549E-01 


2 . 8 72 0 E 00 


2.7858E 00 


2.6273E 00 


15 


2.59n5E-0l 


3.1313E 00 


3 . 0 4 0 0 E 00 


2.8708E 00 


16 


2.4311E-01 


3.3868E 00 


3.2923F 00 


3.H52E 00 


17 


2.2768E-01 


3 . 68 0 5E 00 


3.5815F 00 


3.3940E 00 


18 


1275E-01 


4.0585E 00 


3 . 950 IF 00 


3.7445E 00 


19 


1 . 9833E-01 


4.5281E 00 


4 . 4 0 53F 00 


4.1731E 00 


20 


1 .8442E-01 


5.0873E 00 


4.9465E 00 


4.6815E 00 


21 


1 . 7lnlE-0l 


5.7628E 00 


5.5982F 00 


5.2910E 00 


22 


1.5811E-01 


6 . 6l62E 00 


6.4207F 00 


6. 0586E 00 


23 


1 .4571E-01 


7.7214E 00 


7.4837F 00 


7 . o 48 0 E 00 


24 


1 .3382E-01 


9.1329E 00 


8.R386E 00 


8 . 3053E 00 


25 


1 . 2244E- 0 1 


1.0860E 01 


1.0490E 01 


9.8281E 00 


26 


1 .1156E-01 


1.2890E 01 


1.2422F 01 


1 . 1 5 9 9 E 01 


27 


1 . 0119E-01 


1.5239E 01 


1.4646E 01 


1.3619E 01 


28 


9.1325E-02 


1.7990E 01 


1.7239F 01 


1.5961E 01 


29 


8 . 1965E-02 


2.1205E 01 


2.0247F 01 


1.8650E 01 


30 


7.3H0E-02 


2.4876E ol 


2.3666F 01 


2.1683E 01 


31 


6.4762E-02 


2.8784E 01 


2.7259F 01 


2.4807E 01 


32 


5.6920E-02 


3.2558E ol 


3 . o 6 76F 01 


2.7716E 01 


33 


4 . 95P3E-02 


3.5759E 01 


3.3489F 01 


3 . nOOlE 01 


34 


4.2753F-02 


3 . 8 i 32E 01 


3 , 5455E 01 


3 . i 375E 0 i 


35 


3.6429E-02 


3.9529E ol 


3.64U5F 01 


3.1844E 01 


36 


3 - 0610E-02 


3.9848E ol 


3.6256F 01 


“3.127/E 01 


37 


2 . 5298E-02 


3.8994E ol 


3.4941F 01 


2 . 9651E 01 


38 


2. 0491E-02 


“3 . 6909E 01 


5.2445E 01 


2.6991E 01 


39 


1 .61O1E-02 


3 . 3619E 01 


2.8878F 01 


2.3436E 01 


40 


1 . 2396E-02 


2 . 9266E ol 


2 .4426F 01 


1.9 20”3E 0i 


41 


9 . 1072F-03 


2.4io9E 01 


1 . 9 4 4 9F 01 


1.4659E 01 


42 


f .3244F-03 


1.8522E 01 


1.4278E 01 


1 • 0146E 01 


43 


4 . 0476F-03 


1.2954E 01 


9.3837E 00 


6.1077E 00 


44 


2.2768F-03 


7.804 0 E 0 0 


5.2262F 00 


2 .952 2 E 00 


45 


1 . 0119E-03 


3.7721E 00 


2 . i 714p 00 


9 . 4648E-01 
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TABLE B . 7 

SPATIALLY AVERAGED THERMAL SPECTRA FOR THE CORNELL 

ZPR 1.5:1 CORE 



I 


Energy 


M0D6RA 1 Oh 


ClAI)D i \G 




F U6L 








SPECTRUM 


S?FC 1 9 1 J M 


SPECTRUM 


1 


1 . 6o 93E oo 


5 . 4770E- 


01 


5 . 47 /9f- 


01 


5.41016- 


01 


2 


j . 2 3 1 1 E 0 0 


6.2156E- 


01 


0.1 843F- 


01 


6 . p 7 0 6 E - 


01 


3 


0.6667E-G1 


8 . 1724E- 


01 


8.P963F- 


01 


7 .91086- 


01 


4 


7.79416-01 


l . 0 2 1 7 E 


00 


1 . Ql36F 


U0 


9.9151E- 


01 


5 


6.45filF-0l 


1 .2606E 


00 


1 . 2 48 3f 


0 0 


1.21936 


00 


6 


8 .47316-01 


1 . 4943E 


00 


1.47730 


00 


1 . 43876 


00 


7 


& . 7542E-01 


1.7191E 


00 


1 . ft944F 


00 


1 . 6438E 


00 


8 


-a . 21P2E-U1 


1 .9332E 


oo 


1 . P983F 


0 0 


1 . 83386 


uo 


9 


7 . 8170E-01 


2.1327E 


00 


2 . f)869F 


00 


2 • 00556 


00 


10 


3 . 5119E-01 


2.3190E 


00 


2.7592F 


00 


2-15826 


00 


11 


3.27R6E-01 


2 . 5420E 


00 


2 . 4 6 7 3 F 


00 


2 . 3 4 5 4 1 


00 


12 


3 . 09906-01 


2.6954E 


00 


2 . 6 0 7 5F 


00 


2.45736 


00 


13 


2 . 9244E-01 


2.8287E 


00 


2 . 7 3 1 2 F 


0 0 


2 . 5776E 


00 


14 


?. 7549E-01 


3- . 0 4 6 5 E 


00 


2.9399F 


0 0 


2.7726E 


00 


15 


2 , 59f)5E-0l 


3 .3369E 


00 


3 . 2230F 


00 


3.04366 


00 


16 


?. 4311F-01 


3 . 6 3 o 4 E 


00 


3.5112F 


00 


3.32236 


00 


17 


2,27686-01 


3.9778E 


00 


3 . 8511F 


00 


3 . 6499E 


00 


18 


2 . 1275E-01 


4 . 4 3 6 6 E 


00 


4.29570 


0 0 


4 . o 722E 


uo 


19 


i , 9833E-01 


5 . U 2 4 7 E 


00 


4 . P62PP 


0 0 


4 . fto5 76 


oo 


20 


1 . 8442E-01 


5 . 7 5 2 3 E 


00 


5 . 5615f 


0 0 


5 . 2ft3oE 


0 0 


21 


i . ''lQl E-0l 


6.6668E 


oo 


6 . 43 /9f 


0 0 


6.n845E 


00 


22 


1 . 5 8 1 1 E - 0 1 


7 . 8 5 7 7 E 


oo 


7 . 5 7 7 3 F 


0 0 


7 . i 50 U 6 


00 


23 


1 . 4571E-01 


9 . 4 351 E 


00 


9 , f j 8 2 9 F 


00 


8.85416 


00 


24 


1 .33826-01 


1.14936 


01 


1 . 1 04?f 


01 


1 . o 3 7 6 E 


01 


25 


1 .22446-01 


1 . 4 g78E 


01 


1 . 34V1F 


01 


1.96416 


01 


26 


i .11566-01 


1.7210E 


01 


1 . ft444p 


01 


1 .53556 


01 


27 


i . OU9E-01 


2.0941E 


ol 


1 . 9937F 


01 


1.85386 


01 


28 


o .13956-09 


2. 5381 E 


01 


2 . 4 0 7 2 <= 


01 


2.22876 


01 


29 


p . 1965F-0? 


3 . 0 6 1 3 E 


ol 


2 . 8 8 V 6 F 


01 


2.661/6 


01 


30 


7.3H.0E-U2 


3.6609E 


01 


3 . 43e9F 


0 1 


3.1 5 0 o 6 


01 


31 


ft. 47626-09 


4 . 3 o 7 3 E 


01 


4 . n 2 2 2 f 


01 


3.66036 


01 


32 


5 , 6 9 2 0 F - l) 2 


4 , 9486E 


01 


4 . R898F 


01 


4.i 4 6 9 E 


01 


33 


4.95«3F-U2 


5 . 521 9E 


01 


5 . n 8 i.l 1 F 


01 


4 .55096 


01 


34 


4.27536-02 


5 .90216 


01 


5 . 4 6 .j 3 F 


01 


4.82316 


01 


35 


,7,64296-02 


6.29636 


01 


5 . ft 6 1 3 F 


01 


4.95246 


01 


36 


3 . 0 6 1. f- F - U 2 


6 . 4 3 6 1 E 


01 


5 . ft 9 2 1 F 


01 


4.91036 


01 


37 


7.529PF-U2 


6.37626 


01 


5 . 822 7F 


01 


4 .ft6666 


01 


38 


2 , u 4 Q 1 6-02 


6 . 0 9 9 7 E 


ol 


5.1 4 H 0 p 


01 


4 . 90276 


01 


39 


i . 6toiF-U2 


5 . 6 o 58E 


01 


4 . 8 8 7 4 F 


01 


3 . 723UE 


01 


40 


1 . 23966-02 


4 .91546 


01 


3 . 8 it 3 F 


01 


3 . o 4 4 4 E 


01 


41 


Q. 10726-03 


4.U712E 


01 


3 . f> 6 m 5 F 


01 


2 .71286 


01 


42 


ft. 3244 E- 03 


3 . 1386E 


01 


4 . 2 3 1 2 F 


01 


1.53596 


01 


43 


4 . 0476F-0 3 


2.1971E 


01 


1 . 4 4 -J 7 F 


01 


9 . 79686 


0 0 


44 


?.27ft8E-03 


1.33376 


01 


7 . 8 5 3 2 F 


uo 


4,43636 


0 0 


45 


1 . 0119E-03 


6.3281E 


00 


3 . i 4 v5o 


00 


1 - 37286 


0 0 
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TABLE B .8 

SPATIALLY AVERAGED THERMAL SPECTRA FOR THE CORNELL 

ZPR 2:1 CORE 



I 


ENERGY 


moderator 

spectrum 


CLADDING 
SPEC f RUM 


fuel 

spectrum 


1 


1 . b (J o 3 fc o n 


5 .4898E- 


01 


5 . 48H7E- 


01 


5.4207E- 


01 


2 


1 . 2 3 1 1 F 0 0 


6 . 2 4 4 8E - 


01 


6 , p 0 8 7 e “ 


01 


6. 0946E- 


01 


3 


Q e 66^7fr-0l 


8.2325E- 


01 


8.1 492E- 


■01 


7 . 9 6 0 6 E - 


01 


4 


7.7941E-U1 


1 . Ci 3 0 1 E 


00 


1 . n2Q9F 


00 


9 . 9 8 6 5 E - 


01 


5 


* . 4 5 p 1 F - 1) 1 


1 .2730E 


00 


l.?596E 


00 


1.2298E 


00 


6 


5 . 4 7llE- ui 


1 .5h8E 


00 


1.4924E 


00 


i . 4534E 


00 


7 


4 , 754 2F-U1 


1 . 7 4 2 9 E 


00 


1 .714BF 


00 


1 . 6635E 


00 


8 


4.21O2E-01 


1 ,9644E 


00 


1 . 9253E 


00 


1. 8593E 


00 


9 


7 ■ 8 1 7 [! E - 0 1 


2 . 1725E 


00 


2 . i 2 U 2F 


00 


2 . o 375E 


00 


10 


3 . 5li 9E-01 


2.3689E 


00 


2.3003E 


00 


2.1975E 


00 


11 


7 . 27O6F-01 


2 . 6 o 38E 


00 


2.5160F 


00 


2.3936E 


00 


12 


3 , ii 9o Of-- 01 


2.7681E 


00 


2 . 6668E 


00 


2.5234E 


00 


13 


2 . 9244F-01 


2.9124E 


00 


2 . 7996F 


00 


2.6422E 


00 


14 


p . 7549F-U1 


3 . 1 4 6 U E 


00 


3.n220F 


00 


2 . 850 IE 


00 


15 


?.59o5F-Ol 


3 .4581E 


00 


3 . 3248F 


00 


3.1398E 


00 


16 


? . 4 3 1 1 E - 0 1 


3 . 7806E 


00 


3.6399F 


00 


3.4441E 


00 


17 


?,27f>8E-Cl 


4 . 1711E 


00 


4 . n2U0F^. 


00 


3 . 80 96E 


00 


18 


? .1275E-01 


4 .6967E 


00 


4 . 5266F 


00 


4 . 2 91 0 E 


00 


19 


1 . 9833E-H1 


5.3856E 


00 


5 . 1 862F 


00 


4.9128E 


00 


2 0 


1 . 8442F-G1 


6.2613E 


00 


6. 0233F 


00 


5 . 7 0 0 7E 


00 


21 


i . 7lni E-Ol 


7 , 391 4 E 


00 


7 . f)999F 


00 


6 . 7102E 


00 


22 


1 , 5813 E- 01 


8 . 8 9 2 4 E 


00 


8 . 5 2 7 1 F 


00 


8 . o 462E 


00 


23 


i .4571F-01 


1 . 0 9 0 9 E 


01 


1 . 0438F 


01 


9.R308E 


00 


24 


i . 33R2F- 01 


1.3573E 


01 


1 .2953F 


01 


1.2176E 


01 


25 


i . 2244E-U1 


1 .6971E 


01 


1 . 6151F 


01 


1.5133E 


01 


26 


i . 1156E-01 


2 . 1 1 6 3 E 


01 


2 . n 0 6 8 f 


01 


1.8739E 


01 


27 


1 . 0119F-01 


2 . 6233E 


01 


2.4787F 


01 


2 • 3 0 3 0 E 


01 


28 


Q.1325E-02 


3.2322E 


01 


3 . p 3 7 6F 


01 


2.8124E 


01 


29 


P . 19iS5F-02 


3.9526E 


01 


3 . P932F 


01 


3 . 4019E 


01 


30 


7.3HOF-02 


4.7797E 


01 


4 . 440 OF 


01 


4 . p 68 6 E 


01 


31 


47R2F-0? 


5.6769E 


01 


5 . ?359f 


01 


4.7649E 


01 


32 


5 .6 9 ? 0 E - 0 2 


6.5794E 


01 


6 . 0 1 6 5f 


01 


5.4377E 


01 


33 


4 . 95P3F-02 


7 . 4 q 7 1 6 


01 


6 . 7 0 9 0 F 


01 


6 . o 1 0 2E 


01 


34 


4 . 2753F-02 


8 . G958F 


01 


7 . 2462F 


01 


6.4123E 


01 


35 


3 . 6429F-02 


8.5931E 


01 


7 . S601F 


01 


6.6199E 


01 


36 


3 . 0 6 1 0 F - 0 2 


8 . 8511E 


01 


/ . P 3 / 9 F 


01 


6. 5887 E 


"01 


37 


2 , 52Q8F- u? 


8 . 6270E 


01 


7 . 4250F 


01 


6.3009E 


01 


38 


? . 04O1F-U2 


8 .4 9l8E 


ol 


6.9214F 


01 


5.7580E 


01 


39 


1 . 61Q1E-02 


7 .8403E 


01 


6 . j 5 7 4 F 


01 


4 . 9 97 IE 


01 


40 


i , 23Q6E-02 


6.8997E 


01 


5.1782F 


01 


4 . p 71 0E 


01 


41 


0 . 1072F-03 


5.7299E 


01 


4 . Q 8 1 4f 


01 


3.p762E 


01 


42 


P . 3244F-U3 


4 . 4 2 4 i E 


01 


2 . 9431F 


01 


2. p9l3E 


01 


43 


4 . U476F-03 


3 . 0975E 


01 


1 . R818F 


01 


1 . ?249E 


01 


44 


P.27P8E-03 


1 .8772E 


01 


1 • f) 0 7 4 F 


01 


5.6906E 


00 


45 


1 .0119F-U3 


8.8666E 


00 


3 , q564F 


00 


1 . 7245E 


00 
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TABLE B .9 

SPATIALLY AVERAGED THERMAL SPECTRA FOR THE CORNELL 

ZPR 3:1 CORE 



I 


energy 


MODERATOR 


cladding 




fuel 








SPECTRUM 


SPEC! RUM 


spectrum 


1 


j . 6093E 00 


5 .5o25E- 


01 


5.4979E" 


01 


5.4298E- 


01 


2 


1 .2311E 00 


6.2740E- 


01 


6 . 2299?" 


01 


6.H54E- 


01 


3 


Q.6667E-0-L 


8.2930E- 


01 


8.1943F- 


01 


8 • 0045E- 


01 


4 


7.7941. E-Oj 


1 . 0386E 


00 


1.0275E 


00 


1.0051E 


00 


5 


A . 45nlE-0l 


1.2855E 


00 


1.2694E 


00 


1.2394E 


00 


6 


5.4731E-01 


1 . 5295E 


00 


1 .5062F 


00 


I.4668E 


00 


7 


4 . 7542E-01 


1.7671E 


00 


1.7333E 


00 


1.6814E 


00 


8 


4 . 2192E-U1 


1.9963E 


00 


1 • 9 4 9 1 E 


00 


1.8824E 


00 


9 


3 . 8170E-01 


2.2135E 


00 


2.1503E 


00 


2.Q664E 


00 


10 


3.5H9E-01 


2 . 4 2 0 8E 


00 


2.3379F 


00 


2.2334E 


00 


11 


3.2786E-Q1 


2.6693E 


00 


2.5651E 


00 


2.4384E 


00 


12 


3 . 0990E-01 


2.8464E 


00 


2.7231E 


0 0 


2.5767E 


00 


13 


2 . 9244E-01 


3 . 0 0 52E 


00 


2.8671E 


00 


2.7059E 


00 


14 


?. 7549E-01 


3.2604E 


00 


3 • 1 0 7 7 E 


00 


2.9308E 


00 


15 


?.59o5E-0l 


3 . 6o43E 


00 


3.4384E 


00 


3.2471E 


00 


16 


2.43HE-01 


3.9724E 


00 


3.7952E 


00 


3.5910E 


00 


17 


P.2768E-01 


4.4338E 


00 


4.2404E 


00 


4.0184E 


00 


18 


2.1275E-01 


5.0712E 


00 


4.8493E 


00 


4.5969E 


00 


19 


1 .9833E-01 


5.9324E 


00 


5.6662E 


00 


5.3675E 


00 


20 


1.8442E-01 


' 7 . 0 6 7 3 E 


00 


6.7408E 


00 


6.3797E 


00 


21 


1 . 7lnlE-0l 


8.50O9E 


00 


8.1684F 


00 


7.7201E 


00 


22 


1 .5811E-01 


1. 0640E 


01 


l.OlO^E 


01 


9.5357E 


00 


23 


1 .4571E-01 


1.3451E 


01 


1.2740F 


01 


1.1999E 


01 


24 


1 . 3382E-01 


1 . 7220E 


01 


1 .6262E 


01 


1.5281E 


01 


25 


1 . 2244E-01 


2.2H0E 


01 


2. 0 796F 


01 


1.9484E 


01 


26 


1 .1156E-01 


2.8257E 


01 


2.6453E 


01 


2.4701E 


01 


27 


1. 0119E-01 


3.5016E 


01 


3.3342E 


01 


3.1004E 


01 


28 


9.1325E-02 


4.4981E 


01 


.4 . 1 6 2 6 E 


01 


3.8540E 


01 


29 


R.1965E-02 


5.5868E 


01 


5.1313F 


01 


4.7266E 


01 


30 


7.3H0E-02 


6 . 8387E 


01 


6.2343E 


01 


5.7119E 


01 


31 


6.4762E-02 


8 . 2 0 50 E 


01 


7 .4117E 


01 


6.7450E 


01 


32 


5 . 692QE-U2 


9.5986E 


01 


8.5787E 


01 


7 . 7 5 0 9 E 


01 


33 


4.9583E-02 


1 . 0909E 


02 


9.6262E 


01 


8.6236E 


01 


34 


4 . 2753E-Q2 


1.2036E 


02 


1 . o458e 


02 


9. 2 5 "4 6 E 


01 


35 


3.6429E-02 


1.2889E 


02 


1 . 0966F 


02 


9.5918E 


01 


36 


3. 0610E-02 


1.3381E 


02 


1 ,1087e 


02 ' 


9.5639E 


01 


37 


2.5298E-02 


1.3434E 


02 


1. 0?7lE 


02 


9.1404E 


01 


38 


?. 0491E-02 


1.2997E 


02 


1.0010E 


02 


8.3273E 


01 


39 


1 .6191E-02 


1.2054E 


02 


8.8597E 


01 


7.1902E 


01 


*0 


1 .2396E-02 


1.0644E 


02 


7.3949F 


01 


5 . 8138E 


01 


41 


9 . 1072E-03 


8.8606E 


01 


5.7754E 


01 


4.3530E 


01 


42 


f .3244E-U3 


6T8496E 


01 


4'. j 12 7 e 


01 


2.922^41“ 


01 


43 


4 . 0476E-03 


4.7951E 


01 


2 .5877E 


01 


1.6843E 


01 


44 


2.2768E-03 


2 . 9 0 0 2E 


01 


1 . 3586E 


01 


7.6744E 


00 


45 


1 .0119E-03 


1 .3635E 


01 


5.227QF 


00 


2.27836 


00 
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TABLE B.10 

SPATIALLY AVERAGED THERMAL SPECTRA FOR THE CORNELL 

ZPR 4:1 CORE 



I 


energy 


moderator 


CLADL) 1 




FUfc L 








spectrum 


spectrum 


SPECTRUM 


1 


1 . 6 o o 3 E 0 1) 


5 . 5o90E- 


01 


5 . 5 0 1 5 F - 


01 


5 . 4 333E- 


01 


2 


1 , 23t IE Uq 


6 .2890E- 


01 


6 ,?3b0F- 


•oi 


6.1 2336- 


•01 


3 


9 , 86676-01 


8 .32416- 


01 


8 .21216- 


01 


8 . n?2uE- 


01 


4 


7 . 7941 F- Ji 


1 . 0430E 


00 


1 . f)3u 3F 


0 0 


1 . 0078E 


00 


5 


6 , 45n] E- 0 1 


1 .2920E 


00 


1.2736F 


00 


1 . 9435E 


00 


6 


5.47J1F-U1 


1.53866 


00 


1 .5120F 


00 


1.4725E 


00 


7 


4 . 7542E-U1 


1.77956 


00 


1 .74096 


00 


1 .6989E 


00 


8 


4 . 21O2F-01 


2.01296 


00 


1 . 05b ?E 


00 


1.8918E 


00 


9 


3 . 8 1 7 0 E - 0 1 


2 .23506 


00 


2.16276 


00 


2. f)783E 


00 


10 


3 . 5H9E-01 


2.44846 


00 


2.35356 


00 


2.9483E 


00 


11 


3.2786F-01 


2 . 7 o 496 


00 


2.58546 


00 


2.4576E 


00 


12 


,3.099 0 6- 01 


2.89026 


00 


2.74856 


uo 


2.60076 


OU 


13 


9 , 9244F-01 


3 .05926 


00 


2.90006 


00 


2. 7369E 


00 


14 


2.7549E-01 


3.33046 


00 


3.1 5336 


00 


2.9739E 


00 


15 


2 . 59q5E-01 


3.6990E 


00 


3.50546 


00 


3 . 3 1 0 3E 


OU 


16 


? . 4 3 1 1 F - l) l 


4.10496 


00 


3 . 89b06 


0 0 


3.6364E 


00 


17 


2.2768E-01 


4 . 6265E 


00 


4 .39576 


00 


4 . i 6 5 6 E 


00 


18 


? ,12756-01 


5 . 3 6 G 3 E 


00 


5 . 0914F 


00 


4 .82646 


00 


19 


±. 98336-01 


6 . 3 7 2 IE 


00 


6 , Q 43 56 


00 


5.72526 


00 


20 


1 ,84426-01. 


7.73656 


00 


7 . 32546 


00 


6.9331E 


00 


21 


1 . 7lnlE-0l 


9 .5931E 


00 


9 . Q6216 


00 


8 . 5648E 


00 


22 


1 .58116-01 


1.21546 


01 


1 . j 4516 


01 


1 . n 3 0 5 E 


01 


23 


1 .48716-01 


1.56856 


01 


1.47236 


01 


1 . 3B70E 


01 


24 


i .3382F-U1 


2.04566 


01 


1 .91416 


01 


1 . 79886 


01 


25 


1 .2244E-01 


2.67066 


01 


2.48676 


01 


2*32996 


01 


26 


1 .11566-01 


3 .4641E 


01 


3 . ?0 7 76 


01 


2.99536 


01 


27 


1 . 011.9F-01 


4.4483E 


01 


4.09196 


01 


3.80496 


01 


28 


9 . 13956-02 


5.64756 


01 


5 . 15856 


01 


4.77616 


01 


29 


R. 19656-02 


7 . 07496 


01 


6.40476 


01 


5.8995E 


01 


30 


7.3H0E-02 


8.71786 


01 


7 .82216 


01 


7.16676 


01 


31 


6 .47626-02 


1.05166 


02 


9 . 3344p 


U1 


8.49476 


01 


32 


5 6920 E- 02 


1 .23646' 


02 


1 . 0 83 66 


02 


9.79076 


01 


33 


4.95R3F.02 


1.41226 


02 


1.21926 


02 


1. 09226 


02 


34 


4 .27536-02 


1 .56586 


02 


1 . 3273F 


02 


1-17486 


02 


35 


3 . 6429E-0? 


1.68446 


02 


1 .39266 


02 


1 . 2 181E 


02 


36 


3 . 06i r.F-02 


1.75576 


02 


1.40806 


02 


1.21366 


02 


37 


? . 52086-02 


1.76876 


02 


1.36vJ6p 


02 


1.1571E 


02 


38 


2 .04916-02 


1.71576 


02 


1 . ?623r 


02 


1 . 0501E 


02 


39 


1 .61O1F-02 


1.59476 


02 


1-11176 


02 


9.02186 


01 


40 


1 . 2 3 o 6 F - U 2 


1.41046 


02 


9 . 21966 


U1 


7.94836 


01 


41 


O.1Q72E-03 


1.17536 


02 


7 . 1499F 


01 


5.3989E 


01 


42 


6.32446-03 


9 . (IR916 


01 


5.04986 


U1 


3.59626 


01 


43 


4 , 04766-03 


6 . 3611E 


01 


3.14326 


01 


2.04596 


01 


44 


2 . 2 7 6 8 F - 0 3 


3.84306 


01 


1.63206 


U1 


9.9237E 


00 


45 


1.011 9R-03 


1 .80256 


01 


6 . 93836 


0 0 


2.71836 


00 
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TABLE B.ll 

CORNELL UNIVERSITY ZERO POWER REACTOR INTEGRAL PROPERTIES 



Integral 

Properties 




Water/Oxide Volume 


Ratio 






1:1 


1.5:1 


2:1 


3:1 


4:1 


f 


.92185 


.88785 


.85480 


.79331 


.73687 


V 


1.72300 


1.72867 


1.73165 


1.73461 


1.73605 


T)f 


1.58835 


1.53480 


1.48021 


1.37607 


1.27925 


c m 

0 

n 


1.25149 


1.30644 


1.35120 


1.42173 


1.47787 


5 C 

n 


1.14087 


1.14917 


1.15325 


1.15677 


1.15807 


D, cm. 


.3699 


.3110 


.2770 


.2402 


.2206 


Lj cm. 


2.0086 


1.9356 


1.9253 


1.9659 


2.0298 


|j.sec . 


28.69 


35.17 


41.53 


53.43 


64.37 


v . m/sec. 
m J 


3667. 


3318. 


3134. 


2945. 


2849. 


v , m/sec. 


3802. 


3478. 


3310. 


3146. 


3067. 


V£, m/sec. 


3971. 


3626. 


3446. 


3265. 


3177. 


v h , m/sec. 


3802. 


3426. 


3223. 


3012. 


2902. 
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while Table B.12 gives their temperature coefficients as 
defined by Eq . (7.5). 

The effective cross sections, Eqs . (6.23) and (6.31), 
are given at 20°C for the five cores in Table B.13, the 
temperature coefficients being given in Table B.14. 

Table B.15 lists the 20 °C values of the mean cross 
sections, Eqs. (6.24) and (6.32), while Table B.16 lists 
the temperature coefficients of these cross sections . 
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TABLE B . 12 

CORNELL UNIVERSITY ZERO POWER REACTOR INTEGRAL PROPERTY 

TEMPERATURE COEFFICIENTS 



Temperature 

Coefficients 




Water/Oxide Volume Ratio 




(10" 4o C _1 ) 


1:1 


1.5:1 


2:1 


3:1 


4:1 


f/f 


+ . 266 


+ .546 


+ .813 


+1.37 


+1.95 


nh 


- .310 


-.386 


-.411 


- .429 


-.440 


uf/uf 


- .044 


+ . 160 


+ .402 


+ .941 


+1.51 


6 m /6 m 
n' n 


-2.76 


-3.71 


-4.46 


-5.56 


-6.35 


S C /6 C 

n n 


-1.53 


-1.78 


-1.94 


-2.08 


-2.16 


D/D 


+3.23 


+4.49 


+5.76 


+7.27 


+8 . 38 


L/L 


+5.51 


+6.89 


+7.82 


+9.10 


+9.90 


i t /£ t 


+ .048 


-.411 


-.807 


-1.34 


-1.63 


v /v 
m' m 


+8.32 


+10.1 


+11.2 


+12.6 


+13.3 


v c /v c 


+7.76 


+9.34 


+10.3 


+11.4 


+12.1 


Vf/^f 


+6.92 


+8.55 


+9.58 


+10.7 


+11.5 


. +! 
, X! 

”l> 


+7.76 


+9.63 


+10.9 


+12.3 


+13.0 
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TABLE B . 13 



CORNELL UNIVERSITY ZERO POWER REACTOR EFFECTIVE CROSS SECTIONS 



Effective 
Cross Sections 




Water /Oxide Volume 


Ratio 




(cm -1 ) 


1:1 


1.5:1 


2:1 


3:1 


4:1 


a 


.0219 


.0219 


.0219 


.0219 


.0219 


2 C 

a 


.0139 


.0139 


.0139 


.0139 


.0139 


a 


.3641 


.3680 


.3700 


.3719 


.3728 


z m 

s 


4.1967 


4.0195 


3.9254 


3.8295 


3.7806 


2 C 

s 


.1470 


.1346 


.1283 


.1220 


.1190 


S 


.6685 


,6112 


.5811 


.5510 


.5364 


tr 


2.9376 


2.8927 


2.8688 


2.8445 


2.8322 


AC 

27 

tr 


.1572 


.1451 


.1389 


.1328 


.1298 


zr 

tr 


1.0182 


.9660 


.9385 


.9110 


.8976 


•fttn 

Z t 


4.2186 


4.0414 


3.9473 


3.8513 


3.8025 


2 C 

Z t 


.1609 


.1485 


.1421 


.1358 


.1328 


^f 

2 t 


1.0327 


.9792 


.9511 


.9230 


.9092 


Z f 


.2852 


.2618 


.2637 


.2655 


.2664 


2 n 

a 


.1584 


.1292 


.1095 


.0851 


.0706 


2 h 

s 


2.4138 


2.6804 


2.8600 


3.0800 


3.2103 


2V 

tr 


1.9152 


2.0965 


2 2227 


2.3807 


2.4755 




2.5722 


2,8097 


2.9694 


3.1651 


3.2809 


yll 

Z f 


. 1036 


.0816 


.0667 


.0482 


.0 372 
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TABLE B . 14 

CORNELL UNIVERSITY ZERO POWER REACTOR EFFECTIVE CROSS SECTION 

TEMPERATURE COEFFICIENTS 



Temperature 

Coefficients 




Water/Oxide Volume 


Ratio 




(10" 4o C -1 ) 


1:1 


1.5:1 


2:1 


3:1 


4 : 1 


i m /t m 

a a 


-2.97 


-2.97 


-2.97 


-2.97 


-2.97 


i c /z c 

a 1 a 


0.0 


0.0 


0.0 


0.0 


0.0 


-If 

z 1 7z 

a a 


0.0 


0.0 


0.0 


0.0 


0.0 


i m /z m 
s s 


+2.09 


+2.88 


+3.34 


+3.84 


+4.11 


t c /2 c 

s s 


+8.92 


+9.43 


+10.4 


+11.6 


+12.2 


Z f /Z f 
s s 


+7.06 


+8.60 


+9 .62 


+10.8 


+11.5 


-Atn /«m 
Z tr'"tr 


+1.29 


+1.60 


+1.78 


+1.97 


+2.07 


v c /vC 

^tr'^tr 


+7.19 


+8.54 


+9.39 


+10.4 


+10.9 


i f /Z f 
Z tr />J tr 


+4.07 


+4.76 


+5.19 


+5.69 


+5.97 


• 

yin /'yin 

V^t 


+2.07 


+2.85 


+3.31 


+3.80 


+4.06 


iy±i 


+7.21 


+8.59 


+9.43 


+10.4 


+10.9 


i{,tl 


+4.12 


+4.81 


+5.26 


+5.76 


+6.04 


i\n\ 


-1.63 


-1.85 


-1.99 


-2.13 


-2.21 


ih /g h 


-.06 


+ .43 


+ .78 


+1.35 


+1.63 


th/S 11 

s s 


+1.77 


+2.28 


+2.63 


+3.10 


+3.40 


/z h 


+1.30 


+1.34 


+1.38 


+1.48 


+1.58 
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Table B . 14 (continued) 



Temperature 

Coefficients 




Water/Oxide Volume 


Ratio 




(10' 4o C _1 ) 


1:1 


1.5:1 


2:1 


3:1 


4:1 


i h /s h 

z t /z t 


+1 .66 


+2.19 


+2.56 


+3.05 


+3.37 




-.10 


+ .55 


+1.20 


+2.28 


+3.23 
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TABLE B . 15 

CORNELL UNIVERSITY ZERO POWER REACTOR MEAN CROSS SECTIONS 



Mean 

Cross Sections 




Water/Oxide Volume Ratio 




(cm -1 ) 


1:1 


1.5:1 


2:1 


3:1 


4:1 


z m 

a 


.0131 


.0145 


.0154 


.0164 


.0169 


z c 

a 


.0080 


.0088 


.0092 


.0097 


.0099 


z f 

a 


.2017 


.2233 


.2363 


.2506 


.2582 


z m 

s 


2.5177 


2.6648 


2.7560 


2.8605 


2.9190 


2 C 

s 


.0851 


.0852 


.0852 


.0853 


.0853 


S f 

s 


.3 704 


.3708 


.3710 


.3713 


.3714 


~m 

Z tr 


1.7624 


1.9178 


2.0141 


2.1248 


2.1867 


Ml 
Pt O 


.0910 


.0918 


.0923 


.0929 


.0931 


tr 


.5641 


.5860 


.5992 


.6138 


.6215 


sin 

^t 


2.5309 


2.6794 


2.7713 


2.8768 


2.9359 


K 


.0931 


.0939 


.0944 


.0950 


.0953 


4 


.5721 


.5940 


.6073 


.6219 


.6296 


4 


.1430 


.1588 


.1684 


.1789 


.1844 


z h 

a 


.0917 


.0830 


.0747 


.0622 


.0535 


Ml 

CO 


1.3970 


1.7214 


1.9524 


2.2499 


2.4337 


2 tr 


1.1084 


1.3463 


1.5174 


1.7391 


1.8766 


4 


1.4887 


1 . 8044 


2.0271 


2.3120 


2.4872 


S f 


.0599 


.0524 


.0455 


.0352 


.0282 
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TABLE B.16 

CORNELL UNIVERSITY ZERO POWER REACTOR MEAN CROSS SECTION 

TEMPERATURE COEFFICIENTS 



Temperature 

Coefficients 




Water/Oxide Volume 


Ratio 




( io - 4 o c -1 ) 


1:1 


1 . 5:1 


2:1 


3:1 


4:1 


^m/gm 
a a 


- 11.0 


- 12.7 


- 14.0 


- 15.0 


- 15.7 


Z C /Z C 
a 1 a 


-7 .49 


- 9.13 


- 9.78 


- 10.8 


- 11.6 


z f /z f 
a' a 


- 8.13 


- 9.78 


- 10.9 


- 12.2 


- 12.9 


Z m /Z m 
s ' s 


- 6.16 


- 7.11 


- 7.75 


- 8.51 


- 8.94 


Z c /Z c 
s' s 


+ .16 


+ .14 


+ .12 


+ .11 


+ .10 


Z f /Z f 
s s 


+ .11 


+ .11 


+ . 10 


+ .09 


+ . 08 


Z™ / Z ? 
tr' tr 


- 6.96 


- 8.36 


- 9.28 


- 10.3 


- 10.9 


Z c /Z c 
‘ J tr /Z tr 


- .60 


- .76 


-.92 


- 1.08 


- 1.18 


rf /z£ 

tr' tr 


- 2.83 


- 3.67 


- 4.23 


- 4.92 


- 5.31 




- 6.19 


- 7.15 


- 7.78 


- 8.55 


- 8.98 




- .54 


- .75 


-.90 


- 1.05 


- 1.15 


H r 4 


- 2.79 


- 3.61 


- 4.18 


- 4.85 


- 5.24 




- 8 .46 


- 10.2 


- 11.3 


- 12.6 


- 13.3 


ih /g h 


- 7.69 


- 8.97 


- 9.83 


- 10.6 


- 11.1 


z h /z h 

s s 


- 5.91 


- 7.19 


- 8.01 


- 8.93 


- 9.40 


i h /z h 

IT 


- 6.36 


- 8.11 


- 9.23 


- 10.5 


- 11.2 
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Table B.16 (continued) 



Temperature 

Coefficients 




Water/Oxide Volume 


Ratio 




(10 -4o C _1 ) 


1:1 


1.5:1 


2:1 


3:1 


4:1 




-6.01 


-7.28 


-8.08 


-8.97 


-9.44 


i h c /2 h c 


-7.76 


-8.87 


-9.44 


-9 . 66 


-9.70 



APPENDIX C 



NUMERICAL METHODS IN COUTH 

This appendix describes the computational methods 
employed by the COUTH code in solving Eq . (6.1) for the 
spatially averaged moderator flux ^(E), as well as the 
subsequent calculation of all of the parameters defined 
in Chap. 6. A Fortran listing of COUTH appears in 
Appendix E utilizing the Nelkin water model for the 
scattering kernel. The input to COUTH is discussed in 
Appendix D. 

COUTH, as presently instrumented begins its cal- 
culations with the computation of the energy, velocity, 
and lethargy mesh as described in Sec. 7.1 and listed 
in Appendix B. Inspection of these meshes shows that 
the lethargy mesh has the finest relative spacing in 
the region of maximum interest, .025 eV <£ E < .6 eV, 
while both the velocity and energy mesh tend to have the 
finest spacing at low energy end of the mesh. For this 
reason all calculations in COUTH are done in lethargy, 
u, space utilizing the relationships 



cp(u) = Ecp(E) , 


(C.l) 


a(u'— -u) = Ea(E ' — E) , 


(0.2) 


S(u) = 2(E) . 


(0.3) 
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The heterogeneous effects of the lattice. Sec. 2.3, are 
incorporated into COUTH in terms of the moderator and fuel 
element escape probabilities . The escape probability of 
source neutrons, P (E), is calculated in the subprogram FUNCTION 

b 

PSLOF utilizing diffusion theory as expressed by Eq . (3.4). 
P m (E), the moderator escape probability of neutrons incident 
from the fuel element is calculated utilizing the reciprocity 
relationship, Eq . (3.6), by the subprogram FUNCTION PMODF. 

The linear extrapolation length into a black hole, d , 
necessary for these calculations is approximated by — 7 

d 0 (E)2^ r (E) - X - °' 6 - 229 + 0.7104 , (C.4) 

l+2.482™ r (E)r 0 

where r Q is radius of the black hole. This expression equals 
4/3 at r Q = 0, Zaretsky's— ^ value of 0.889 at r Q = 1, 
and 0.7104 for r Q approaching infinity. Leslie's 
SPECTROX method—^ corresponds to the special case of A = 2/3. 

The fuel rod escape probability, P r (E), is calculated 
from the parabolic flux approximation, Eq . (4.8). This com- 
putation is carried out by the subprogram FUNCTION PR0DF. 

The subprogram SUBROUTINE PESCF calculates the fuel element 
escape probability, P e (E), utilizing the thin region transport 
theory embodied in Eq . (4.21). 

In the explicit calculation of escape probabilities, the 
modified Bessel functions of the first and second kind and of 
order 0, 1, and 2 are required. Calculation of the functions 
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I Q (z), I-^(z), ^(z), K Q (z), K-^(z), and K^(z) is accomplished 

by the subprogram SUBROUTINE BESS . 

The modified Bessel functions of the first kind, I (z), 

86 / 

are computed in the following manner. — 7 Noting that for 

any given z, I n (z) approaches zero as n gets very large, we 

set I n ( z ) equal to zero and I n _^(z) equal to 1x10 ^ where 

for any z, n is determined as follows. For z < 45, n = 

[2z+10], where [x] is the largest integer contained in x. 

For z > 45, n is taken equal to 100. Lower order I's are 

87 / 

then calculated utilizing the recurrence relationship — 7 

Vl( z > ■ T + ■ < c - 5 > 



Recognizing that the I's calculated in this manner bear the 
correct relationship to each other, but lack normalization, 
the results are normalized through the expression 

oo n 



e = 



I (z) 
m v 



I D ( Z ) + 2 



m=-oo 



Y. 

m=l 



i_(z) 

m 



(C . 6) 



This same technique is utilized in the calculation of the 
I's utilized in the computation of the Nelkin water kernel, 

Eq. (5.15). 

The zero order modified Bessel function of the second 

88 / 

kind, K Q (z), is computed using the following approximations. — 7 



— ^ A. M. Weinberg and E. P. Wigner, The Physical Theory of 

Neutron Chain Reactors , p. 765, The University of cKTeago 
Press, Chicago" 1958 . 

— ' 1 D. F. Zaretsky, PICG 5, 529 (1955). 
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For z > 5 : 



IT 


1/2 -z 


l 2 

1 - -I 


l 2 -3 2 l 2 • 3 2 -5 2 , 1 


1^1 


c 


1 TTSx + 


2 I ( 8x ) 2 3 ! (8x) 3 



K o (z) % 



(C.7) 

the number of terms used in this expansion being set equal to 
2[z] . For z < 5, 

50 

ViW * < c - 8 > 



K q (z) 



m=0 



where 



8 = -(7 + log z/2). 



8 n = 2 , 



(C .9) 



7 = Euler's Constant = 0.577215, . .., (C.10) 



(C.ll) 



and , 



o _ m-1 R 

0 “ ZZ 0 -l 

m m in- 1 



m > 2 (C . 12) 



The first and second order modified Bessel functions of 



the second kind are calculated using the Wronskian relation- 
Ship 



K (z)I 
m v ' 



m+1 



(z) + K 



m+1 



< Z >V Z > 




(C . 13) 



— ' D. C. Leslie, "The SPECTROX Method for Thermal Spectra in 
Lattice Cells," AEEW-M211, United Kingdom Atomic Energy 
Authority (1962) . 

— ^ F. D. Federighi and D. T. Goldman, "KERNEL and PAM-Programs 
for Use in the Calculation of the Thermal Scattering 
Matrix for Chemically Bound Systems," KAPL-2225, Knolls 
Atomic Power Laboratory (1962) . 
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The calculation of the slowing down source, S(E), as 
given in Eq . (7.1), and of the monatomic gas model for the 
scattering kernel, Eq . (5.9), requires the evaluation of 
the error function for various arguments . The error 
function, ERF(z), is defined by 



ERF(z) 





(C . 14) 



In the COUTH code, ERF(z) is calculated by FUNCTION ERF 

90 / 

in the following manner. — ' For |z|>4. 



ERF(z) 




(C . 15) 



while for 2 |z| <4, 



ERF(z) 






1.3 

77 



1.3.5 



+ 



1.3. 5. 7 
2.2 4 z 8 




(C . 16) 



For |z| <2, 



ERF(z) % — 




z 2(n-l) z 40 

"(2n-l) + 2TIDTJT4TT 



(C.17) 



87 / 

- — ' H. B. Dwight, Tables of Integrals and Other Mathematical 
Data , p. 189" The Macmillan Company, New York, 1961 . 

88 / 

— M. Goldstein and R. M. Thaler, "Recurrence Techniques for 
Calculation of Bessel Functions," Mathematical Tables 
and Other Aids to Computation, 13, 102 (1959). 
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Numerical integration in the lethargy space used in 
COUTH is carried out by the subprogram FUNCTION TRAP. The 
trapezoidal rule is utilized to simplify the calculations 
with the uneven lethargy spacing encountered. Integration 
of an arbitrary function of lethargy, f(u), which has been 
evaluated at the discrete points f(u^), i = 1,2,3, ...n, 
is approximated by 



u 



n 



f (u)du 



'X. 



U 



n-1 



i=l 



f (u i )+f (u i+1 ) 



Au. 



(C - 17) 



The lethargy intervals, Au^, are defined by 



Au^ = u^-^-u^, i = 1,2, 3,... n-1 . 



(C . 18) 



Calculation of the scattering kernel using the monatomic 
gas model, Eq . (5.4), is performed by the subprogram SUB- 
ROUTINE FKERN. Only the down scattering in energy (or up- 
scattering in lethargy) portion of the kernel is calculated 
in FKERN, the upscattering portion being calculated by 
detailed balance, Eq. (2.5), in the main program of COUTH. 
Calculation of the monatomic gas kernel in the form of a 
45 by 45 matrix as needed with the lethargy mesh used takes 
about minutes of computer time. If the Nelkin water 
kernel is used, then the scattering kernel, Eq. (5.15), is 
calculated with the subprogram SUBROUTINE WATER. Again, 
only the downscattering portion of the matrix is computed, 

89 / 

■ — M. Goldstein and R. M. Thaler, see Footnote 88. 
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the upscattering portion being calculated by detailed 
balance in the main body of COUTH. The actual methods 
used in the calculation of the Nelkin water kernel are 
discussed in Sec. 5.2. Since the calculation of 45 by 45 
matrix representing the Nelkin water kernel takes approxi- 
mately 15 minutes at a given temperature, the down- 
scattering portion of the matrix for a physical tempera- 
ture of 20 °C has been reproduced on cards so that numerous 
lattice studies may be made at this temperature without 
necessitating the recomputation of the water kernel. Con- 
trol cards read into COUTH (see Appendix D), determine 
whether the downscattering portion of scattering kernel 
shall be calculated using the Nelkin model, or whether the 
Nelkin results for 20°C shall be read directly in. 

For calculations at temperatures above 20 °C, the 

average cosine of the scattering angle for water, p(E), 

-4 

is decreased by 1x10 for every centigrade degree increase 

in temperature above 20°C at all energies. The slowing 

down source of neutrons, S(E), is calculated according to 

Eq . (7.1), while the moderator scattering cross section, 

Z m (E), is calculated using Eq . (2.2) and as discussed in 

Sec. 7.1. The neutron balance equation, Eq . (6.1), 

rewritten below in terms of the lethargy mesh is then solved 

91/ 

using an iterative process . — 

— ! H. B. Dwight, p. 136, see Footnote 87. 

— ^ M. J. Poole, M. S. Nelkin, and R. S. Stone, "The Measure- 
ment and Theory of Reactor Spectra," Progress in Nuclear 
Energy , Series I, Vol. 2, p. 99, Pergamon Press, New York, 
London, Paris, and Los Angeles, 1958. 
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u 



/ 



c 




Z(Uj— •-u.j.) cp (Uj)dUj 




o 



(C.19) 




The cutoff lethargy, u^, is 7.372 corresponding to a low 
energy of .00101 eV with the thermal cutoff energy, E c , of 
1.609 eV being the zero lethargy point. The relative fuel 
absorption f / [ 1- f ] is calculated from the escape probabilities 
in accordance with Eq . (2.29). 

An initial guess is made that the moderator spectrum 
is a constant in lethargy: 



This guess is then normalized to satisfy neutron conservation 
as expressed in Eq . (2.4). 



Wi) - 1 , 



i = 1, 2, 3, . . . , n = 45 . 



(C .20) 




o 



(C . 21) 



The normalized flux, '^(u^) is then substituted into the right 
hand side of Eq . (C.19) and the next generation approximation 
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to the moderator flux, cp^(u^), is calculated. This process 
is repeated utilizing the general expressions 




u 



S( Ui ) + / Z(uj— «-u 1 )'5”(Uj)d(Uj) 



2>i> + *S<«i> 



m, 



1 + 



f(u ± ) 

r-f(u.) 



and. 



i = 1,2.. .n. 



(C . 22) 







L 



i=l, 2, . . .n . 
(C . 23) 



The equation is considered solved when the successive 
iterations satisfy the following convergence criterion at 
every other lethargy point: 




i 1,3, 5, 7, . . . , n * 
(C .24) 



The convergence criterion, e, utilized in the COUTH cal- 
culations is 1x10"^. In general solution of the integral 
equation, Eq . (C.19), has taken about 13 iterations with use 
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of the monatomic gas model, while the use of the Nelkin water 
kernel with its more complicated fine structure usually 
requires about 27 iterations for solution. 

The fuel and cladding spectra, the mean velocities, 
the disadvantage factors, the lifetime, the thermal utilization, 
7)f , L, D, and the averaged cross sections are calculated in a 
straightforward manner as defined in Chap. 6. 

An entire lattice calculation, including the computation 
of the escape probabilities, the spatially averaged thermal 
spectra, and all of the parameters defined in Chap. 6 takes 
about 40 seconds of computer time with use of the Nelkin water 
kernel, and about 30 seconds with the monatomic gas model. 
Calculational times for the kernel are not included in these 
figures . Compiling time for the Fortran source program of 
COUTH runs about two and one half minutes . 



APPENDIX D 



INPUT-OUTPUT FORMAT FOR COUTH 

A Fortran statement of the COUTH code for use with the 
Nelkin water kernel as a scattering model is listed in 
Appendix E. This section describes the input required 
and the output received for a typical lattice calculation. 

The input data is listed below, the first column being the 
item or data card number, the second column is the Fortran 
FORMAT as it appears in COUTH, the third column is the 
Fortran list of the input variable, and the fourth column 
is a description of each variable. 

1. (14) IND Identification 

number. If IND =1, 
the previously cal- 
culated down 
scattering portion 
of the Nelkin water 
kernel at 20°C will 
be used in the cal- 
culation. If IND > 2, 
the Nelkin water 
kernel will be cal- 
culated at the tempera- 
ture specified in 
Item 2 . 



- 160 - 
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2. (2F10.6,I5,E12.4) T 



RH(J> 

N 



EPSI 



3. (13F6.1) 



XA(I) 



4. (13F6.1) 



XF(I) 



The moderator tempera- 
ture in degrees abso- 
lute . 

The moderator density 
at T in gm/cc. 

The number of lethargy 
intervals . N = 44 gives 
mesh in Table B.l. N 
must be < 54. 

Convergence criterion, 
e, Eq. (C . 24) . 

235 

U absorption cross 

section in barns . 

235 

U fission cross 

section in barns. 



5. (7E11.3) 



SK(J,I) Downscattering portion 

of Nelkin model 
scattering matrix for 
20°C . 



AMUBAR(I) Average cosine of the 
scattering angle in 
water, p(E^) . 



6. (20F4.3) 
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7. (11F7.4) CF(I) 

8. (8F8.4.E16 .4) DEN 

ENR 

RO 

RCI 

RF 

AMC 

ssc 

SAC 

FNC 



Correction factors to 
the SK(JjI) diagonals 
see Sec . 7.1. 

Fuel rod density in 
gm/ cc . 

Fuel enrichment in 
weight per cent . 

Fuel element diameter 
in inches . 

Fuel cladding inner 
diameter in inches . 

Fuel rod diameter in 
inches . 

Cladding scatterer 
mass in amu . 

Cladding free atom 
cross section in barns . 
Cladding 2200 m/sec 
absorption cross 
section in barns . 
Cladding number density 
in atoms /barn -cm . 

80 character title. 



9. ( 10A8) 



TITLE 
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10. (2F10.5) W(j) The nominal water to 

oxide volume ratio 
for the lattice under 
consideration . 

AP The lattice pitch or 

rod to rod spacing (in.) 
in a hexagonal lattice. 

COUTH is instrumented to reread a series of items 10 until 
a negative value of W(j) is encountered. If this occurs, 
items 8 and 9 will be reread with the subsequent read in of a 
series of items 10. Calculation will stop when the computer 
encounters a negative value of W((), item 10, followed directly 
by a negative value for DEN, item 8. This allows one to cal- 
culate the scattering kernel at a given temperature, then con- 
sider a fuel element of certain specifications and calculate 
parameters for all of the lattices utilizing this fuel element. 
New fuel element specifications are then read in and the 
lattices utilizing this fuel element are analyzed. There is 
no limit on the number of fuel elements or the number of 
lattices for each fuel element type which may be considered. 

The output for a typical 1 lattice calculation using a 
45 point energy mesh is outlined below. The first column 
represents the page number after the execution of the source 
program, the second column gives the Fortran list of the output 



co lo vo 
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variable, and the third column gives a description of the 
output variable. 



1. I 
E(I) 
DE(I) 
V(I) 

DV (I) 
U(I) 

DU (I) 

2 . BLANK 
SK(J, I) 



8. E (I) 

SX(I) 

ST(I) 

SS25(I) 

SS28(I) 

SS16 (I) 

AMUBAR(I) 



Index identifying lethargy mesh points . 
Energy mesh point, E^. 

Energy interval, AE^ - E i+1 -E i . 

Velocity mesh point, v^ 

Velocity interval, Av^ . 

Lethargy mesh point, u^ . 

Lethargy interval, Au^. 



Scattering matrix, a(U(J) ^-U(I)), 

J running horizontally and I running 
vertically . 

Energy mesh point, E^. 

Moderator scattering cross section in barns . 

Moderator transport cross section in barns . 
235 

U scattering cross section in barns. 

238 

U scattering cross section in barns . 

0 16 scattering cross section in barns . 
Average cosine of the scattering angle, 

]I(E i ) , for water. 
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9. E(I) 


Energy mesh point, E^ . 


XA(I) 


235 

U absorption cross section in barns . 


SF(I) 


235 

U fission cross section in barns. 


SIGMAF(I) 


Fuel absorption cross section in cm 


FSIGMAF(I) 


2f(E i ). 

Fuel absorption cross section in cm"'*'. 


AA(I) 


1 + a(E^), a(E^) = fuel capture to fission 
ratio . 


SN(I) 


Slowing down source of neutrons, S(E), 
Eq. (7.1). 


10. E(37 ) 


E(37) = .0253 eV. This page shows E = kT 
values for the variables used in the 


PR 


cladding calculations. 
P^, Eqs . (4.8), (4.19). 


FJ(j) 


J o , Eq. (4.22). 


TT 


Cladding thickness in cm. 


FO 


Solid angle factor, Eq. (4.26). 


FIN 


Fraction of neutrons scattered in, Eq. (4.48) 


SIGS 


Cladding scattering cross section in cm 


SIGT 


Cladding total cross section in cm 


G 


g, Eq. (4.45). 


S 


T, Eq. (4.35). 


PU 


P u , Eq. (4.43). 


PUNIF 


P unif ’ E< 5- < 4 - 44) ' 
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RSUM 

TAU 

RH(^SUM 

CTE 



r in + r out 



11. BLANK 

12. BLANK 

13. 



, Eq . (4.46). 
t, Eq . (4.40). 

P. 



+ Ppil# ., Eq. (4.47). 
idinj 

absorption rate ratio, Eq . (4.55). 



in K out J 

Cladding absorption rate to fuel element 



W0 


Input 


item 10 . 


T 


Input 


item 2 . 


TITLE 


Input 


item 9 . 


E(I) 


Energy, E ± . 


FLUX (I) 


^(E.) 


, Eq . (6.1). 


FC(I) 


*c< E i> 


, Eq. (6.4). 


FS (I) 


? f (E i ) 


, Eq. (6.3). 


AB(I) 


Z“(E) 


1 + 1-£(e) 


W(j) 


Input 


item 10 . 


T 


Input 


item 2 . 


TITLE 


Input 


item 9 . 


E(I) 


Energy, E ^ . 


PS (I) 


W 


, Eq. (3.4). 


PM(I) 


V E i> 


, Eq. (3.6). 


PR(I) 


P rr<V 


, Eq. (4.8), 
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PE(I) 

TU(I) 

RFA(I) 


P e (E.), Eq. (4.21). 
f(E t ), Eq. (3.8). 
f/(l-f), Eq. (3.8). 


RTE(I) 


Rod abs. rate/element abs . rate, Eq. (4.53) 


FJ(]) 


J c (E i ), Eq. (4.22). 


15 . W(J) 


Input item 10 . 


T 


Input item 2. 


TITLE 


Input item 9 . 


F 


f, Eq . (6.15). 


ETA 


v- 


ETAF 


r)£, Eq . (6.16). 


DIS 

DISC 

DIFFC 


5®, Eq. (6.8). 
6^ Eq. (6.9). 
D, Eq. (6.38). 


DIFFL 


L, Eq. (6.39). 


AL 

ALl 

VM 

VF 

VC 

HV 

FN25 

FN28 


H t , Eq. (6.17). 

£ = £ /r)£. 

m t 

v m , Eq. (6.6). 
v f , Eq . (6.6). 
v c , Eq. (6.6). 

Eq . (6 .33) . 

?3S 

U number density. 

238 

U number density. 


FN16 


Oxygen in UO 2 number density. 
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FNC 


Cladding number density, input item 8. 


FNM 


Moderator number density. 


N 


Input item 2 . 


K 


Number of iterations required to solve 
Eq. (6.1). 


EPSI 


Input item 2 . 


RH(j> 


Input item 2 . 


ENR 


Fuel enrichment in atom per cent . 


AP 


Input item 10. 


Rl 


Equivalent cell radius, cm. 


RO 


Fuel element radius, cm. 


VOLM 


2 

Moderator area, cm . 


VOLC 


2 

Cladding area, cm . 


VOLF 


2 

Fuel rod area, cm . 


R 


Actual water to oxide volume ratio. 


TITLE 


Input item 10 . 


16. 


Page 16 lists the effective cross sections 
for each region of the cell, Eq . (6.23), 
and for the homogenized cell, Eq . (6.31); 
and the mean cross sections for each 
region, Eq . (6.24), and for the homogenized 
cell, Eq. (6.32). 
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Pages 12 through 16 are repeated for each lattice con- 
sidered, while consideration of a fuel element with new 
specifications will cause pages 10 and 11 to be printed 
with the subsequent printing of pages 12 through 16 for 
each lattice considered. 
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FORTRAN LISTING OF COUTH 

This appendix lists the Fortran statement of the COUTH 
source program as instrumented for calculation with the 
Nelkin water kernel as the scattering model. SUBROUTINE 
FKERN, which calculates the monatomic gas model for the 
scattering kernel, is also listed. For calculations with 
the monatomic gas model, SUBROUTINE WATER, FUNCTION ONE and 
FUNCTION TWO are replaced by SUBROUTINE FKERN . In addition 
the main program of COUTH must be modified to read in the 
values of the transport cross section of water as cal- 
culated with the Nelkin model for use in the calculation 
of the moderator escape probabilities . 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



Appendix E 



PROGRAM COUTH 

D I MFNS I ON U <55 ) 9 E <55 ) 9 PS (55 ) ?PM(5F) ,Pf(55) * R FA < 55 ) * TU (55 ) » AB ( 55 ) « 
1 SN ( 55 ; , Sl< <55 <>55 ) * SX ( 55 ) , aL-PhA (55 ) , Ph I N (55 ) 9 DU ( 55 ) * A A (55 ) i FEA (55 ) 9 
2GAMMa ( 55 ) , OFLTa <55 ) , SaUPhA (55)i SPh IN ( 55 > * eLuX ( 55 ) »ZZ (55 ) 9 SXA (55 ) » 
3XA ( 55 ) •> y.F <55 ) * FS ( 55 ) » ST ( 55 ) 9 aMURAR ( 55 ) , KN ( 5 ) , 1 M ( 5 ) , CF ( 55 ) , 

4S IG^AF ( 55 ) , FSIGMAF ( 55 5 9 OF ( 55 ) , oV ( 55 ) 9 V ( 55 ) , SS25 (55 ) * SS28 <55 ) , 

5SS 1 6 ( 55 ) 4 0 <55 ) , DO <55 ) < T I TLp ( 1 0 ) •> PR ( 55 ) 9 F J 0 (55 ) » R TE ( 55 ) 9 FC ( 55 ) * 



6CT E ( 55 ) •; F8 ( 55 ) <TRB(55) 

ORDER OF DATA INPUT IS AS FOLLOWS 0 



T , RHO , N , EPS I 

XF 

s K 

AMijBAR 

CF 

DE n 9 ENR 9 R 0 9 Rc I 9 Rf 9 AMc 9 SSC 9 S AC 9 FNc 

title 



14 

2(FlO»6)9 I 5 9 2X 9 ^ 1 0 0 2 
1 3 ( E6 0 1 ) 

I3(E6o 1 ) 

7 (El l 03 ) 

20(E4»3> 

1 1 (F7«4 > 

8 (F8.4 ) 9 4X 9 El 2 0 4 
10 <A 8 ) 



U/0 , AP 2 ( E 1 0 o 5 ) 

IND = 1 9 NeLKIN WATER KERNEL FOR 20 DEG. C READ IN. 

IND = 29 NELKIN WATER KERNEL CALCULATED AT To 
T = MODERATOR TEMPERATURE IN DEGREES KELVIN, 

RHO IS THE MODERATOR DENSITY IN GM/CC. 

N = NUMBER OF LETHARGY INTERVALS, 

EPSI = CONVERGENCE CRITERION FOR THE ITERATION. 

XA IS h'ONECKS U235 ABSORPTION CROSS SECTION IN BARNS. 

XF IS HONECKS U235 FISSION CROSS SECTION IN BARNS. 

SK IS THE SCATTERING MATRIX. 

SK IS SIG m A(U(J) TO U ( i 5 ) /'UN 1 T LfTHARGY OR f( I )*SIGMa(F(U) TO 
E < 1 ) ) /UN IT EV . SI GMA IS IN BARNS . 

AMUBAR IS THE AVERAGE COSINE OF THF SCATTERING ANGLE. 

CF IS THE CORRECTION FACTOR TO THE SK DIAGONALS. 

DEN IS THE FUEL DENS I T Y IN GM/CC. 

ENR IS THE FUEL ENRICHMENT EXPRESSED AS WEIGHT PER CENj. 

RO IS THE FUEL ELEMENT DIAMETER IN INCHES. 

RCI 7 S THE FUEL CLADDING INNER DIAMETER IN INCHES. 

RF IS THE FUEL ROD DIAMETER IN INCHES. 

AMc IS THE CLADDING SCATTERING MASS IN AMU. 

SSC IS T HE CLADDING FREE ATOM SCATTERING CROSS SECTION IN BARNS. 

sac is the cladding 2200 m/sfc absorption cross section in barns, 
fnc is the cladding concentration in ATOMS/RARN-CM, 

WO IS THE NOMINAL WATER TO URANIUM OXIDE VOLUME RATIO. 

AP IS THE LATTICE PITCH IN INCHES. 

19 FORMAT I 1 H 1 ; 

1 99 format (t>o H i scattering kernel 9 sigma (u<jj to u<i>)/unit lethargy 9 

1 KT = F7 . 5 v 6H 9 J - ! 2 9 4 H TO 129/) 

200 FORMAT ( I 4 s 2X 9 1 O ( E U o 3 ) ) 

203 FORMAT (37 !l WATER TO U r ANiuM OXIDE VOLUME RaT I O , p5„ 1 , 7H TO 1.0 9 / 9 2 
14H ABSOLUTE TEMPERATURE = F8.29/.27H THERMAL UTILIZATION (F) = F8. 

25 9 / 9 7H ETA = F8.5</o9h ETA*F = F8.5./.33H MOdE r ATOR DISADVANTAGE P 
3ACT0R - F6 i 5 s / 9 32H CLADDING DISADVANTAGE R ACTOR = F8.5,/,30H DIEFU 
4S I ON COEFFICIENT (CM) = F8.49 /925H DIFFUSION LENGTH (CM) = F8.4./t 

525 h lifetimes (Microseconds) 9/9 iox. i ih thermal = F 8 . 59 /, iox.8h mfa 

6N = f8,5(/< 27H AVERAGE VELOCITIES ( M/SFC ) 9 / 9 1 OX 9 1 3H MODERATOR = f8 
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7 o 2./,I0X.8h FUfL = F8 0 2./, 10X, 12h C l ADD! n G = f8 . 2 • / . 10Xi20H hOMOge 
8NIZED CELL = F8.2) 

204 FORMAT ( / / , 3 0 H N|jMbe r Op E N E R GY INje r VALs = !4,/,24H N|jMBE R OF I TE 

1RATI0NS = !4t/t 25H CONVERGENCE CRITERION = El 1 «3,/i29H MODERATOR D 

2ENSITY (GM/CC) = F8.5./,24H FUEL DENSITY (GM/CC) = F8„5,/,34H FUEL 
3 ENRICHMENT (ATOM PERCENT) = F8o5./.26H LATTICE PITCH (INCHES) = F 
48 o 5 . // t 1 1H RADII (CM)'/' 10x, 19 h EQUIVALENT CELL = F8.5./, 10X. 16H f 
5UFL ELEMENT = F8 » 5 « / . lOXt I2 h FUEL ROD = F8.5.//.14H A R EAS <CM**2), 
6/ ♦ lOXi 1 3 h MODERATOR = f8 „ 5 , / . 1 OX , 1 2h CLaDDI n G = f8 „ 5 , / , 1 OX * 1 2h FUE 
7L ROD = F8,5.//t46H ACTUAL WATER TO URANIUM OXIDE VOLUME RATIO = F 
88. 5. / ) 

206 FORMAT ( 4X i 3 1 HEFFECT I VE CROSS SECTIONS ( 1 /CM ) , 29X . 26HMEAN C R OSS SE 
ICTI ONS ( 1 /CM ) , / ,4X » 31 ( 1 H“ ) . 29X , 26 < 1 H“ ) ./« 2 (8X. 1 OH ABSORPT I ON , 4 2X ) . / 
2.2(1 2X. 9hM0d E RaT0R,f 1 Oo 5.29X ) . /.2 ( 1 2X . 4 HFUeL . F 1 5 . 5 , 29X ) , /, 2 ( 1 2x . 8 h 
3CLADD ING.F1 1 .5 . 29X ) . / . 2 (8X . 1 OhSCATTfR I Ng . 4 2X ) . / . 2 ( 1 2x . 9 hM0dE r AT0R , 
4f10 o 5.29X) i/i 2 ( 12X. 4hfUF l .F l5.5.29X)./.2( 1 2X. 8 HCLADD I Ng , F 1 1 . 5 . 29x ) 
5 * / . 2 ( 8X . 9 HT r A n S r ORt . 4 3X ) . / . 2 ( 1 2X . 9HM0dE r AT0R , f 10.5.29X),/«2(12X.4H 
6FUFL. F 15.5, 29X ) , /, 2 ( 1 2X , SHCLADD I N G , F 1 1 .5, 29X ) , / , 2 ( 8X , 5hT0t aL , 4 7X ) , 
7/ » 2 ( 12X,9hMOoeRATOR,f 10 0 5,29X ) ,/,2( 1 2X , 4HFUFL , F 1 5 . 5 , 29X ) ,/,2 ( 1 2X.8 
8HCLADD I NG. FI 1 .5.29X ) . / .2 <8X , 7HF I SS I ON , 4 5X ) ♦ / . 2 ( 1 2X , 4HFUEL . F 1 5 . 5 . 29 
9X ) ) 

207 FORMAT (//.2(8X.31HH0M0geNiZED CFLL c r OSS SEC T 1 ONS , 2 1 X ) . / ♦ 2 ( 1 2X . 1 0 
lHARSORPT iON.f9o5.29X) ./,2( 1 2X , 1 OhSc ATTe r I Ng . F9 . 5 . 29X ) ./.2( 12x.9hTR 
2 AN SPORT . Fl 0 0 5, 29X ) ,/ , 2 ( 1 2X . 5 hT 0T AL . F 1 4 , 5 . 29X ) . / . 2 ( 1 2X . 7HF ! SS I ON , p l 
32 o 5 . 29X ) ) 

208 FORMAT </»34H NuMrE r DENSITIFS ( NuMppR/p a r N _ cM ) . / . 1 OX . 8H U235 = El 
12.4./. 10X.8 h U238 = p 1 2 . 4 , / , 1 OX , 1 Oh OXYGEN = e 1 2 . 4 , / , 1 OX ♦ 1 2h CLADD 
2 1 NG = E12.4./* lOX. 1 3H MODERATOR = E12.4) 

334 FORMAT ( 1 1 < F7 « 4 ) ) 

815 FORMAT (6H I . 3X .6HFNeRgV , 7X . 7HDELTA E.6X.8hVEL0cITY.2X.7hDELta 

1 V»3X.8hLfTHA r GV. 8H DF l TA U , /// . ( 2X • I 4 . 2X . 2 ( p 1 1 . 3 , 2X ) . 2 ( f 8 . 1 , 2x ) . 2 
2 ( F7 . 4 . 2X ) ) 

1000 FORMAT (/. ( IX. !4.7(E15o4) ) ) 

1001 FORMAT ( 1 3 ( F6 ° 1 ) ) 

1004 FORMAT (37h1WATE r TO URANIUM OXIDF vOLuMp R AT I 0 , F5 . 1 . 8H TO 1.0.. 

1 24 H ABSOLUTE TEMPERATURE = F6.1) 

1005 FORMAT (2(F10.R>) 

1006 FORMAT (2 <F1 0.6 ) . I5.2X.E1 0.2 ) 

1007 FORMAT ( 8 < F8 . 4 ) . 4X . E 1 2« 4 ) 

1008 FORMAT ( 1 0 ( AS ) ) 

1009 FORMAT ( / . 1 0 ( A8 ) ) 

1100 FORMAT ( 5H 1 I .4X.6HFNE r GV. 7X. 15HU235 ABSORPT I ON , 2X . 1 2HU235 FlSSl 

ION, IX ♦ 15HFUFL ABSORPT I ON , 1 X , 1 2HFUEL F I SS I ON , 5X . 9h 1 + ALPHA , 4X , 1 2hS 

2L0w INg D°WN, / ,23X , 4 ( 1 3HC R °SS SECT I ON, 2X ) . 1 9X , 6HS0u R CE . / . 26X . 2 ( 7H (B 
BARNS ) . 8X ) . 2 (6H ( 1 /CM ) , 9X ) , /, 39 ( 3H- ) ) 

1101 FORMAT ( 5h 1 I , 4X . 6HFNERGY . 7X , 1 4HH20 SCATTE R I N G , 2X , 1 3HH20 TRANsPO 
1RT.2X. 12HU235 SCATTER , 3X . 1 2HU238 SC AT TE r . 4X , 1 1 HO 1 6 Sc ATTE R . 7x . 6 hMu 

2 RA r . /.23X.5 ( 1 3HC r 0SS SfCTiON, 2X ) , / , 26X , 5 ( 7H ( B A R NS ) .8X ) , / . 39 (3H~ 

3 ) ) 

1102 FORMAT ( 3H0 I ,8X.6HENe r GY,9X,9hM0dE r AT0R,6X,8HCLaDDINg,7X.4HFUEL, 1 

1 3X . 1 6HARS0RPTI ON C R OSS . / . 27X . 3 ( 8HSPFCTRUM , 7X ) . 2X , 1 8HSECT I ON (1 /cM 
?H20 ) . /. 32 < 3H- ).//. < I4.6X,4(E12.4.7X),Fl0o4) ) 

1103 FORMAT OHO I , 8X ,6HE n E r GY , 20X , 2H r S . 8x , 2h r M , 8x . 2hP r , 8x , 2HPE . 8x . 2HTU 

1 ,7X.3h r FA.7X,3h r TF.6X.5h j 0uT./.38(3h - ) , // . ( I 4 , 6X . F 1 2 . 4 , 8x . 8 ( F 1 0 
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1212 
4000 
4444 
9108 
9 1 99 
9298 



2.4 )) ) 

FORMAT (14) 

FORMAT ( 7 ( E 1 1 .7) ) 

FORMAT (20 (F4.3) ) 

FORMAT <FX,4HJ = * 1 0 ( 12 . 9X ) . /, 3H 

FORMAT ( I 4 » 2X . F ( E 1 1 o3) ) 

FORMAT (5X»4HJ = . 5 ( I 2 * 9X ) . / OH 

READ 1212. IND 

READ 1 006. T.RHO.N.EPSI 

NM=N+ 1 



I ) 
1 ) 



READ 1001 , (XA( 1 ) . 1=1 ,NM ) 

Read 1001 . (xf< i ) * i=i .nm> 

Read 4000 , ( (sk ( j, 1 ) , j = 1 , 1 ) , 1 = 1 , nm > 
read 4444 . ( amubar ( 1 ) . 1 = 1 .nm ) 

READ 334 . ( CF ( I ) . I = 1 , NM ) 



T = 8 « 6 1 67E— 5*T 
OM 1 =0.205 



OM2=0 .481 



OMR=0 .06 
EM= 1 8 o 



EMR=2.32 

EMV =3./( 1 o-l o/18«-l o/2.32 ) 

FA 25= 1 0.0 
F A28 = 8 o 3 
FA 1 6 = 3 o 76 

AM25=235« 117/1 .00896 
AM28 = 238 . 125/1 .00896 
AM 1 6= 16.0/1 .00896 
RM=238. 1 25/235 .117 
ANU=2. 43 

START CALCULATION Of ENERGY. VELOCITY. LFTHARGY NESH. 
ENDS WITH STATEMENT 94. 

8=0 • 1 *LOGF (10.) 

U( 1 ) = 0.0 

DV ( 1 ) =0 o i *2200 . 

V ( 1 ) = 44 0 o 
V ( 2 ) = V ( 1 )+DV( 1 ) 

DV(NM )=0o0 
DU ( NM ) =0. 0 
DE (NM ) =0.0 
DO 90 1=2.34 

DV ( I ) =DV ( 1 ) 

90 V ( I + 1 ) =V ( I ) +DV ( I ) 

DO 91 1=35. N 

DV ( I ) =DV ( 1 )*EXPF(B*FLOATF( 1-34 ) ) 

91 V ( 1 +1 ) =V( I )+DV ( I ) 

DO 71 1 I = 1 . NM 

D ( I ) = V ( I ) 

711 DD ( I ) =DV ( I ) 

DO 92 I = 1 .NM 
V( I ) =D ( NM+ 1 — I ) 

92 E ( I ) = ( V ( I )**2 )*. 0252977/ ( 2200. **? ) 

DO 93 I = 1 . N 

93 DV ( I ) =DD (NM- I ) 



CALCULAT 1 ON 
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DO 94 1=2, NM 

U( I ) = LOGF ( E ( 1 )/E( I ) ) 

DU ( 1 - 1 ) = U ( I ) -U ( I - 1 ) 

94 DE ( I - 1 ) =E ( I - 1 ) -E ( I ) 

PRINT 19 

PR I NT 8 1 5 , < I * E ( I ) , DE ( I ) » V < I ) , D V < I ) , U ( I ) , DU ( I ) , I = 1 , NM ) 

PRINT 19 
DO 712 1=1,5 
KJ=10o*FL0ATF< I ) 

IM ( I ) =K J 
712 KN ( I ) =KJ— 9 

I F ( I ND — 2 ) 222 , 223 ,223 

223 CALL WATER(E,T,N,SK ) 

DO 224 1=1, NM 

224 AMUBAR ( I )=AMURAR( I )-( 1 0 . *T- o 2526 ) /8 • 6 1 637 
222 DO 110 J=2,NM 

JK= J- 1 

DO ! 1 0 I = 1 , JK 

88= ( ( E ( I ) /F ( J ) )**2 )*EXPF ( (E ( J )~E ( I ) )/T ) 

1 1 0 SK( J, I >=SK< I , J)*BB 
DO 951 L=1 ,4 
Nl =KN (L ) 

M I = I M ( L ) 

PRINT 1 99 » T , N I , M I 
PRINT 9198, ( I , I=N1 ,MI ) 

951 PRINT 200 , ( I , ( SK( J, I ) , J=N I , M I ) , I = j ,NM ) 

Nl =4 1 
Ml =45 

PRINT 1 99 , T , N I , M I 
PRINT 9298 , ( I , I =NI , M I ) 

PRINT 9199,(1, (SK< J, I ) , J = NI ,MI ) , I = 1 , NM ) 

DO 445 I = 1 , NM 
445 SK< I , I >=SK( I , I )*CF< I ) 

DO 119 J= 1 , NM 
DO 116 I = 1 , NM 
1 16 SXA ( I ) =SK ( J , I ) 

119 SX< J)=TRAP(SXA,DU,N) 

C SX IS THE SCATTERING CROSS SECTION OF H20 IN BARNS, 

DO 447 1=1, NM 

SS25( I ) =FA25*SCR0SS ( E ( I ) ,T, AM25) 

SS28 ( I ) =FA28*SCR0SS ( E < I ) , T , AM28 ) 

SS 1 6 < I >=FAl6*SCROSS(E( I ) , T , AM l 6 ) 

SS25 , 28 ,16 ARE THE SCATTERING CROSS SECTIONS Of U235 , U238 , O 1 6 
IN BARNS* 

4 47 ST ( I ) = .656*SQRTF< ,0253/E( I ) )+SX( I )*( 1 0 -AMUBAR ( I ) ) 

C ST I S THE TRANSPORT CROSS SECTION OF H20 IN BARNS, 

DO 448 I = 1 , NM 
DO 448 J= 1 , NM 

448 SK ( J , I )=SK(J, I )*RH0*o60247/18»02 
PRINT 1101 

PRINT l 000, ( j , f( i ) , SX ( I ) ,ST ( I ) ,SS25 < I ) » SSP8< I ) , SS 16(1), AMUBAR ( I ) , 
1 I = 1 , NM ) 

TE = T/E |V1 +OM j / ( p „ *E m v ) -t-OM2/E M V+OMR* ( j „/ (ExPF ( OMR/j ) “1 . ) + ,5 ) /E^R 

c te is the effective temperature o 
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DO 8 I = 1 , NM 

ST ( I ) =ST ( I ) * RHO* o 6024 7/ 1 8 » 02 
SX( I )=SX( I )*RHO*o60247/lSo0? 

C Sl< iST i AND SX ARE NOW IN (I /CM)® 

A ! = SQRTF ( E ( I )/TE) 

8 SN( I )=E( I ) *ERF ( A 1 ) 

BETA=TRAP(SN, DUiN ) 

DO 9 1=1 i NM 

9 SN ( I )=SN( I )/BETA 

C SN IS THE NORMALIZED SOURCE 

T=T/8o6167E-5 

714 Read ioo7,oen,fnr,ro,rci .rf> amc,ssc»sac»fnc 

IF (DEN ) 41 , 4 1 , 71 8 

718 READ 1 008 , T ! TLE 
AMC = AMC/ 1 o 00096 
R0=R0*2 o 84/2 o 
RF=RF*2® 84/2 o 
RCI=RCI*2®54/2o 

ENR=lOO»*ENR*RM/(ENR*RM+< 100®-ENR) ) 

FN25 = ENR*.6 0247*DEN* 1 ® 0E~2/?70 . 26 
FN28=FN25*( 1 00,-ENR l/ENR 
FN l 6 = 2 »0 l * .60247*DEN/270 c 26 
FNM=.60247*RHO/18«02 
DO 713 I = 1 , NM 

S IGMAF ( I ) =FN25*XA ( I ) +FN28*2® 71 *SQRTF ( .0253/E ( I ) ) 

FS I GMAF ( I ) =FN28*XF ( I ) 

PR ( I ) =PRODF (FN25 .FN28® S IG^AF (I ) , FS I GMAF ( ! ) , QFN, Rf ,RC I «R0 « 

! SS25 ( I ) i SS28 ( I ) ,SSl6< I ) ) 

713 AA ( I ) =S IGMAF ( I l/FSIGMAF ( I ) 

PRINT 1100 

PRINT l 000 , ( i , p ( i ) , XA ( I ) * XF ( I ) ♦ S igmaf ( I ) ♦ FS I r,MAF ( ! ) ♦ AA ( 1 ) ♦ SN ( i > , 
1 I = 1 (NM ) 

CALL PESCF ( E* FNc» SSC .SAC* AMc .RC I .RO.NM.T , PR <FOO<PE<CTE ) 

DO 47 I = 1 »NM 

47 RTE ( I )=FJ0 ( I )* ( i .-PR ( I ) )/ (PR ( I )* ( 1 .-PE ( I ) ) ) 

RTE IS THE RATIO OF THE NUMBER OF NEUTRONS ABSORBED IN THE FUEL 
ROD TO THE NUMBER A8S0RBED IN THE FUEL ELEMENT. 

301 READ 1005. WO. AP 

IF(WO) 714.714,141 
141 PRINT 19 

PRINT 1004, WO, T 
PRINT 1009, TITLE 

R i = SQRTF ( ( ( AP*2 o 84 ) **2 )*1 . 7 3208/ ( 2 « *3 . l 4 1 fo )) 

Y=R1 /RO 

V0LM = 3 o 1 4 1 59* ( (Rl**2)-<R0**2) ) 

VOLF = 3 • 1 4 1 89* ( RF**2 ) 

VOLC = 3 • 1 4 1 89* ( ( R0**2 ) - ( RC 1**2) ) 

V0LCELL = 3 • 1 4 159* ( Rl **2 > 

R= VOLM/VOLF 
RC=VOLM/VOLC 
2 DO 7 I = 1 , NM 

PS ( I ) =PSLOF (E ( I > .RO , Y . RHO, ST ( I ) ) 

PM ( I ) =PMODF ( E < I ) , RO , Y , RHO , ST ( I ) ) 

TU< ! ) = ( ( 1 » — PE ( I ) ) *PS ( I ) )/( 1 « -PM ( ! )*PE ( I ) ) 



no n o o o o 



Appendix E 176 

tu is the Ratio of the number of nfutRons absorbed in the fuel 
ELEMENT TO THE number ABSORBED in THE CELLo 
Rea ( i ) = tu ( i )/< 1 9 -tu < i ) ) 

REA IS THE RATIO OE THE NUMBER OF NEUTRONS ABSORBED IN THE FUEL 
ELEMENT TO THE NUMBER ABSORBED IN THE MODERATOR. 

FEA ( I ) = REA( I )*o656*Rh 0*.6024 7*SQRTE ( e 0253/E ( I ) 1/18.02 
FEA IS THE FUEL ELEMENT ABSORPTION CROSS SECTION IN 1 /CM 0 
7 AB<I) = U.+REA<!)) * o656*RhO*„ 6024 7*SQRTF ( .0253/E < I ) ) /I 8 . 02 

AB IS THE TOTAL ABSORPTION CROSS SECTION JN THE CELL PER CM H20. 
BEGIN ITERATIVE SOLUTION Of INTEGRAL EQUATION, SOLUTION ENDS ON 36 
DO 112 1=1. NM 

112 PHI N ( I ) = 1 . 000 
C INITIAL GUESS AT SPECTRUM IS A CONSTANT 

K =1 

DO 13 I = 1 , NM 

1 3 ALPHA ( I ) =PH I N ( I ) *AB ( I ) 

8 = TRAP ( ALPHA, DU. N ) 

DO 23 I = 1 , NM 
23 PH I N ( I )=PHIN< I )/B 

14 DO IB I = 1 , NM 
DO 16 J= 1 .NM 

16 GAMMA ( J ) =PH I N ( J )*S< ( J . I ) 

15 DELTA ( I )=TRAP(GAMMA,DU.N) 

DO 27 1=1, NM 

27 SPH I N ( I ) = ( SN ( I )+DELTA( I) ) / ( AB ( I ) +SX ( I ) ) 

DO 28 I = ] , NM 

28 SALPHA ( I )=SPH1N( I )*AB ( I ) 

C=TRAP (S ALPHA, DU. N) 

DO 29 1=1, NM 

29 SPH IN ( I )=SPHIN( I ) /C 

DO 32 1=1, NM , 2 

RA = SPH I N ( I )/PHIN( I ) 

A BSD I F = A 8 SE (RA- 1 . > 

IF ( ABSDIF-EPS I ) 32,32,38 

32 CONTINUE 
GO TO 36 

38 K=K+1 

IF (K-lOO) 39,39,36 

39 DO 60 J= 1 , NM 

60 PH I N ( J)=SPHIN( J) 

GO TO 14 

36 DO 37 1=1, NM 

37 FLUX ( I ) =SPH I N ( I ) /E ( I ) 

DO 70 1=1 , NM 

70 ZZ ( I ) =FEA ( I )*RTE ( I ) *SPH IN ( I ) 

Cl=TRAP(ZZ, DU , N ) 

DO 71 1=1, NM 

71 ZZ ( I ) = AB < I )*SPHIN < I ) 

C 2 = TRAP ( ZZ , DU * N ) 

F=C1/C2 

DO 72 1=1 , NM 

FS ( I ) =SPH I N ( I ) *R*EEA ( I ) *RTE C I ) /S I GMAE ( I ) 

FC < I )=FEA< I >*CTF < I )*SPHIN( i ) *Rc*SQRTF (F ( I ) ) / ( SQ r TF ( . 0253 ) *FNc*SAC ) 

72 ZZ< I )=FEA< I )*RTE< I )*SPHIN< I )/AA< I ) 
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C 3 =TRAP(ZZ,DU.N) 

ETAF=ANU*C3/C2 
ET A = ET AF/F 
DO 73 I = 1 ,NM 

73 ZZ ( I )=SPHIN( I )/V< I ) 

C 4 = TRAP ( ZZ i DU . N ) 

DO 74 I = 1 i NM 

74 ZZ ( I ) =FS ( I )/V ( I ) 

C 5 = TRAP ( ZZ . DU » N ) 

DO 75 1=1 .NM 

75 ZZ ( ! ) =FC ( I )/V ( I ) 

C 6 =TRAP (ZZ.DU.N) 

DI S=C4/C5 
D I SC =C6/C5 
C7 = TRAP ( SPH I N . DU . N ) 

C 8 =TRAP(FS,DU.N) 

C9=TRAP (FC, DU .N) 

VM=C7/C4 
VF=C8/C5 
VC=C9/C6 
DO 77 I = 1 .NM 

77 ZZ ( I )=SPHIN( I ) *<>656*RH0* e 60247*SQRTF ( . 0253/E < I ) )/ 18.02 
ClO=TRAP(ZZ« DU « N ) 

DO 78 I =1 .NM 

78 ZZ< I )=FS( I )*SIGMAF< I ) 

Cl i=TRAP(ZZ.DU.N) 

DO 79 I = 1 .NM 

79 ZZ ( I )=FC( I ) *SAC*FNC*SQRTF ( .0253/E< I ) ) 

C12=TRAP(ZZ.DU.N) 

DO 80 1=1. NM 

80 ZZ( I )=SPHIN( I )*ST< I ) 

C13=TRAP(ZZ.DU.N) 

DO 81 1=1 .NM 

81 ZZ ( I >=SPHIN( I )*SX( I ) 

C 14=TRAP(ZZ.DU.N) 

DO 82 I = 1 .NM 

82 ZZ ( I ) =FS ( I )*(SS25( I )*FN25 + SS28 < I )*FN28+SSl6< I )*FNl 6 ) 
C!5=TRAP(ZZ.DU.N) 

DO 83 I = 1 .NM 

8 3 ZZ ( I )=FS( I ) * ( SS25 ( l )*FN25*( 1 . “2 . / ( 3 .*aM25 ) )+SS28( ! )*fN28*( 1 .“2./ 

1 ( 3 «*AM28 ) )+SS 1 6 < I )*FNi 6 * <lo-?o/(?«*AMli 5 ))) 

C16=TRAP(ZZ.DU,N) 

T=T*8o61 67E— 5 
DO 86 ! = 1 .NM 

F 8 ( I ) = ( SPH I N ( I )*VOLM+FS ( 1 )*VOLF + FC ( I ) *VOLC ) /VOLCELL 

TRp ( I ) = (ST ( 1 )*S p HI N ( I ) *vOLM+ (ps ( j j *S !G M AF ( I ) +Z.Z ( ! ) )*V 0 l F+ (FC (!)*(( 
1 SAC*fNc*SQRTf ( e 0253/E ( I ) ) ) + ( FNc*SSC*SC R OSS ( F ( 1 ) . T. aMc )*(l*“2./(3.* 
2 AMC ) ) ) ) ) )/ ( VOLCELL*FB ( I) ) 

86 ZZ ( I ) =FB ( I ) /TRR ( I ) 

C20=TRAP(ZZ.DU,N) 

C2l=TRAP(FB.DU.N) 

D!FFC = C20/(3.*C21 ) 

DO 84 I = 1 .NM 

84 ZZ < ! ) =FC < I )*FNC*SSC*SCR0SS (E<! ) . T . AMC ) 
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C17=TRAP(ZZ*DU*N) 

T=T/8»61 67E~5 
C 1 R=C 1 7* ( 1 o —2 o / ( 3 o * AMC ) ) 

DO 85 I = 1 *NM 

85 ZZ< I )=FS( I )*FSIGMAF ( I ) 

C1<9 = TRAP ( ZZ * DU * N ) 

SBAM=C 1 0/C7 
SBAF = C 1 1 /C8 
SBAC=C 1 2/C9 
SBSM=Cl 4/C7 
SBSF = C 1 5/C8 
SBSC=C 1 7/C9 
SBTRM=C 1 3/C7 
S8TRF = (Cl 1+C16 )/C8 
SBTRC = (C12 + C18 )/C9 

sbtm=sbam+sbsm 

SBTF =SBAF+SBSF 

SBTC=SBAC+SBSC 

SBFF = C 1 9/C8 

V0 = 2200 o 0 

SHAM=SBAM*VM/VO 

SHAF=SBAF*VF/VO 

SHAC = SB AC*VC/VO 

SHSM=SBSM*vM/VO 

SHSF=SBSF*VF/VO 

SHSC = SBSC*VC/V 0 

SHTRM=SBTRM*VM/V0 

SHTRF=SBTRF*VF/VO 

SHTRC=SBTRC*VC/VO 

SHTM=SBTM*VM/VO 

SHTF=SBTF*VF/VO 

SHTC=SBTC*VC/VO 

SHFF=SBFF*VF/VO 

DZ=C7*V0LM+C8*V0LF+C9*V0LC 

SBA= ( Cl 0#V0LM+C1 l*VOLF+Cl 2*VOLC )/DZ 

SBS= ( C 1 4*V0LM+C 1 5*V0LF+C 1 7*V0LC ) /DZ 

SBTR = ( C 1 3*V0LM+< Cl 1 +C 1 6 ) *VOLF+ (C 1 P+C 1 8 )*VOLC )/DZ 

SBT=SBA+SBS 

SBF=C 1 9*V0LF/DZ 

HV=DZ/ ( C4*V0LM+C5*V0LF+C6*V0LC ) 

SHA=SBA*HV/VO 

SHS=SBS*HV/VO 

SHTR=SBTR*HV/VO 

SHT=SBT*HV/V0 

SHF=SBF*HV/V0 

D I FFL =SQRTF ( D I FFC/SBA ) 

AL= ( ( C4*V0LM+C5*V0LF+C6*V0LC )*1 oOE+4 )/ (C2*V0LM ) 

AL 1 =AL/ETAF 
DO 76 ! = 1 *NM 
FC( I )=FC( I )/E( I ) 

76 FS ( 1 ) =FS ( I ) /F ( I ) 

40 PR 1 NT 1 102* < I ) *FLUX( I )»FC< I ) «FS( ! ).AR< I ) » I = 1«NM) 

PRINT 1 004 * W O * T 
PRINT 1009, TITLE 



b \ x 



B x.*c i ■ r: x/l •-} 
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PRiNt 1 1 03 « ( I »F( I ) 1 PS ( I ) ,PM ( I ) , PR( I ) , P E ( I ) , TU ( I ) , r FA ( I) » r TE ( I ) * 

1 F JO ( I M = 1 ,NM) 

PRINT 203, WO, T , p, ETA , ETAF , D I S , D I SC » 0 IFFC , D I FF>- , AL * A*- 1 , VM , VF, VC » HV 
PRINT 208»FN25iFN28«FN16'FNC ,FNM 

PRINT 2 04 iNtK, EPS I , RhO , DFN , ENR , AP , R l , Ro , R F * VOLM , yOLc , VOLp , R 
PRINT 1009, TITLE 
PRINT 1004, WO, T 

PRINT 206«SHAM i SBAM * SHAF«SBAF« SHAC ’ SBAC * ShSM, SbSM , ShSF * SpSF< SpSc * 

1 SBSC , SHTRN , SBTRM , SHTRf , SBTRf , SHTRc 1 SRT r C. , SHTM , SBTM , SHTF, SBTF , SHTC * 
2SBTC,SHFF,SBFF 

PRlNT 207* SHA ’ SBA i ShS , SpS , SpTR , SrTR, ShT , Sb"1" »SHF , SBF 
PRINT 1009, TITLE 
PRINT 19 
GO TO 301 
41 CONTINUE 
END 



FUNCTION PROoF ( FN25 , pN28 , S I G^A A , S I GM AF , DE N , Rf , R c I , R 0 , SS25 , SS28 , 

1 SSI 6 ) 

THIS CALCULATES THE ESCAPE PROBABILITY FROM THE FUEL ROD USING 

THE PARABOLIC FLUX APPROXIMATION, SEE K«B. CADY THESIS, PG 118, 

FN 16=2®01#« 60247*DEN/270 .26 

SS=SS2B*FN2B+SS28*FN28+SSl 6*FN 1 6 

SA=S I GMAA 

ST = SS + S A 

SF = S I GM AF 

C=SS/ST 

A=ST*RF 

B= 1 o-C 

IF <A-.5 ) 1,1,2 

1 CALL BESS ( A , F I 0 *F I 1 ,F I 2 ,FKO ,FK1 , FK2 ) 

F 1 0 = A* < F I 1*FK0 + 1.-A*(FI 1 *FK i + f I 0*FK0 ) ) 

FI 1=(FI 1 *FK 1+1 .-2.*F10)/3. 

FMOO= 1 o —2 o *F 1 1 

PR = 1 o -2 o*A*B* ( 1 . — FMOO )/ ( 1 C*FMOO ) 

GO TO 3 

2 CALL BESS(A,FI0,FI 1 ,FI2,FK0,FKi ,FK2) 

F 1 0 = A* ( F I 1*FK0+1 .-A* (FI 1 *FK 1+FI0*FK0) ) 

FI 1 = (FI 1 *FK 1 + 1 .-2.*F1 0 )/3. 

F20=<FI2*FK0-1 o+2 o*Fl 0+2.*F 1 1 )/3, 

F2 1 = ( ( F I 2*FK l /A ) +F 1 1-F201/5. 

F2?= ( (F I2*FK2/ ( A**2 )) + (l./(3.*< A**2 ) ) )“2.*F21 )/7. 

FMOO= 1 »-2.*F 1 1 
FM02= 1 .-4.*Fl 1 +8.*F21 

FM22= 1 .-6«*F1 1 +24.*F21+24.*F22“2./ ( A**2 ) 

FMO= 2 o * A* ( l .-FMOO ) 

FM2=A* ( 1 ® -FM02 ) 

G= 1 .-C*(4.*FMOO+4.*FM22~6o*FM02) + (C**2 ) * ( 4 . *FM00*FM22“3 . * 

1 (FM02**2 ) ) 

PR=1 o-B*(FMO-C*(FMO*(4«.*FM22 - 3o*FM02)+6**FM2* (FMoO“FMo2) ) ) /G 

3 PRODF= ( 1 o-RF/RCI )+RF*PR/RC I 
END 
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SURROUT I Ne PrScF ( E * F n C * ssc * SAC * AM C * Re * RO * NM , T * PR * F-JO ,Pe,CTE) 

DI MENS I ON E (55 ) *PR (55 ) .FJO (55 ) »PE ( 55 ) * CTE ( 55 ) 

C THIS CALCULATES THE ESCAPE PROBAB1LTY FROM THE FUEL ELEMENT GIVEN 

C THE ESCAPE PROBABILITY FROM THE FUEL ROD. THIN REGION THEORy 

C IS USED AS DEVELOPED IN K„B. CADY THESIS* APPENDIX J. 

19 FORMAT < 1 H 1 ) 

50 FORMAT <10h ENERGY = E12.4*/*31H FUEL ROD ESCAPE PROBABILITY = F8. 
1 5 * / * 1 9H ROq CURRENT OuT = E8 C 5*/*35H FUEL ELEMENT ESCAPE PROBABlLl 
2TY = F8 e 5 * / * 27h CLADDING THICKNESS (CM) = F8.5*/,5H FO =F8.5,/,25H 
3 FRACTION SCATTERED IN = F8.5*/*20H SIGMA SCATTERING = E8.5,/,15H 
4SIGMA TOTAL = E8.5*/,5H G = E8.5,/*5H T = E8,5*/,6H PU = E8.5,/*9h 
5 PuNlF = F8.5*/»14H RlN + ROuT = F8.5*/,7H TAU = F8 0 5./*18H RhOiN 
6+ RHOOUT = F8°5*/*7H CTE = F8.5) 

PI =3 o 14159 
GAm=«5772157 
T=T*8o6167E-5 
ZE=RF/RO 

tt=ro-rf 

FK = 1 o- <ZE**2 ) 

F0=2.*ASINF(ZE )/PI 
F 1 = ZE 

F2=F0+2o*ZE*SQRTF (FK)/PI 
F3--3 o*ZE*< 1 ,-(ZE**2 )/3o )/2o 
Dl =PI* (ZE **2 ) /4 « 

Fl.AM = LOGF( 4 o/SQRTF(FK ) ) 

D2= (2o/3o )*( 1 .-o 75*FK-3«* (FLA M-.25)*(FI<**2)/1 6.-3 «*< FLAM- ( i i ./ 

1 1 2 o ) ) * ( r K**3 ) /64 o ) 

FLl=(F2-4e*Dl/PI )/Fl 

FL2= ( 8o*F3/3.-2.*FK*Fl-4.*D2 )/F 1 

Cl = (P I *F0/2 o+ZE*SORTF (FK ) )/ ( 2.*ZE ) 

C2=D2/ ( ZE**2 ) 

FLL1=4.*C1/PI-ZE 

FLL 2 = 2 .*FK+Bo* (ZE**2 ) /3 . -4 o * ZE*C2 
FIN=ASINF<2«*ZE/( 1 °+ZE) )/PI 
DO 1 0 I = 1 , NM 

SIGS = SSC*FNC*SCROSS <E ( I ) *T* AMC ) 

S I GT=S I GS+ENc*SAC*SQRTF ( . 0253/E ( I ) ) 

S 1 GA = S I GT— S I GS 

S=F 1 ~S I GT*R0*F l *FLl +. 5 * ( <SIGT*R0 )**2 ) *F 1 *FL2 
X=SI GT*7T 

P|J= .5*X* ( GA m +LOgF (X )“1 »5) + l •“( X**2 )/6 • + ( X**3 )/48 ( X**4 )/360 • + 

1 (X**5)/2880. 

PUNIF=PU/( 1 » — ( S I GS* ( 1 o — PU ) /S I GT ) ) 

I F ( X— o 05 ) 31*31*32 

3 1 G= ( 1 o ~ F 1 ) ~ 2 o * < 1 o-F2)*SIGT*RO+Bo *( 1 .~F3)*( ( S I GT*RO ) ** 2 ) /3o 
GO TO 33 

32 G=— .021 2“o0266*L0GF (X )/ (ZE**2.2 ) 

33 CONTINUE 

RSUM = SIGS*( 1 o — G — S ) *PUN I F/S I GT 

T A U = i .-SIGT*RO*FLLl +. 5 * ( (SIGT*R0 )**2)*FLL2 

RHOSUM=SIGS* ( 1 .-TAU )*PUNIF/S IGT 

R I N=F I N*RSUM 

ROUT=RSUM-R IN 

RHO I N=F I N*RHOSUM 
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RHOOUT=RHOSUM— RHO IN 

F JO ( I )=PRC I ) * ( S+R I N ) / ( 1 o —RHO I N*PR ( I ) ) 

PE ( I ) =G+R0UT+F JO ( I ) * ( TAU+RHOOUT ) 

CTF ( I ) =S 1 GA* ( ( 1 - ~G~S )+FJO ( I ) * ( i . ~T AU ) ) / ( ( 1 <. _P F (I ) )* ( S I GA + < S I GS* 

1 PU ) ) ) 

IF ( 1-37) 10,40, 10 

40 PRINT 19 

PRINT 50 , E < I ) ,PR ( I ) ,F JO ( i ) ,Pe ( I ) , TT, F0,FlN, S I GS , S I GT , G , S , P U • p U N I F , 
1 RSUM , TAU , RHOSUM , CTE ( I ) 

PRINT 19 
10 CONTINUE 

T=T/8»6167E~5 

END 



FUNCTION ERF (2 ) 

ERF(Z) CALCULATES the ERROR FUNCTION, 

AZ= ABSF (Z ) 

IF ( AZ-4o ) 85,86,86 

86 ERF = A Z/Z 
RETURN 

85 IF (AZ-2. ) 87,88,88 

87 T 1 = 1 o 

T 2 = T 1 * (Z**2)/3o 

T3=T2* ( Z**2 ) *3 • / ( 2 • *5 - ) 

T4 = T 3* (Z**2 )*5 ®/2 1 • 

T5 = T4* (Z**2 ) *7-/36- 
T 6 = T5* (Z**2 ) *9-/55 - 
T7=T6* (Z**2 )*1 1 o/78. 

T8=T7* (Z**2 )*1 3»/l 05- 
T9=T8* (Z**2 )*1 5-/1 36 o 
TlO = T9*(Z**2)* 17-/1 71 - 
T l l =T1 0* ( Z**2 )*1 9-/21 0. 

T l 2 = T 1 1 * ( Z**2 ) *2 1 • /253 - 
T 13 =Tl 2 *(Z** 2 )*23-/300o 
T l 4 =T 1 3* < Z**2 ) *25-/35 1 - 
Tl5=T14* (Z**2 ) *27-/406- 
T 1 6=T 1 5* (Z**2 ) *29-/465- 
T 1 7 = T 1 6* (Z**2 ) *3 1 -/528- 
T 1 8=T 1 7* ( Z**2 )*33-/595« 

T 1 9=T 1 8* ( Z**2 )*35-/666® 

T20 = T 19* (Z**2 ) *37-/74 1 - 
T21 =T20* (Z**2 ) *19 -5/820 - 

ERF = 2 - * » 564 1 896*Z* ( T 1 ~T2+T3~T4+T5-T6+T7-T8+T9~T 1 0+T 1 1 — T 1 2+T 1 3 
1 — T 1 4 + T 1 5-T1 6+T1 7-T1 8 + T1 9-T20+T21 ) 

RETURN 

88 ZA = 2 - * ( Z**2 ) 

Si = 1 - 

S2= 1 • /ZA 

S3 = 3 o / ( Z A**2 ) 

S4= 1 5 - / ( Z A**3 ) 

S5=52 -5/ (ZA**4 ) 

fRf=( 1 o ~ - 564 1 896* ( S 1 ~ S2+S3 - S4 + S5 ) / ( AZ*FX p F ( Z**2 ) ) )*(AZ/Z) 

END 
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SUBROUTINE BESS(Z,C,D,E,F.G,H) 

C BESS CALCULATES MODIFIED BESSEL FUNCTIONS OF THE FIRST AND 

C SECOND KIND. 

C C= 10 ( Z ) ,D = I 1 (Z ) ,F = 12 (Z ) , F = KO (Z ) , G = K I (Z ) ,H=K2<Z ) 

DIMENSION B ( 1 00 ) ,BI< ( 3 ) 

I F ( Z~ 45 « ) 1.1.2 

1 N=XF1XF(2.*Z + 10. ) 

GO TO 3 

2 N = l 0 0 

3 S ( N- 1 ) = 1 .0E-40 
DO 8 I =N, 1 00 

8 B ( I ) =0 oO 
SUM=2 • 0E-40 
NJ=N- 1 
DO 4 L=2.NJ 
NL=N-L 

B ( NL ) =FLOATF (NL )*2« *B ( NL+ 1 )/Z+B(NL+2 ) 

4 SUM--SUM+2. *B (NL ) 

SUM = SUM— B < 1 ) 

FACTOR= ( EXPF ( Z ) )/SUM 
DO 5 L= 1 . N 
B(L)=FACTOR*B(L) 

IF ( B ( L ) — 1 « OE— 307 ) 6.6.5 
6 B(L)=0.0 

5 CONTINUE 
C=8< 1 ) 

D=B ( 2 ) 

E= B ( 3 ) 

IF (Z-5. ) 81,82.82 
62 IF (Z-6« ) 110,111,111 

1 1 0 KL= 1 0 

GO TO 140 

111 IF ( Z— 7 o ) 112,113,113 

1 12 KL= 1 2 

GO TO 140 

1 13 1 F ( Z— 8 ® ) 114,115,115 

114 KL= 1 4 

GO TO 140 

1 IB IF CZ-9. ) 116,117,117 

116 KL =16 

GO TO 140 

117 IF(Z-lOo) 118,119,119 

1 1 8 KL= 1 8 

GO TO 140 

119 KL= 10 

GO TO 140 
140 W= 1 oO 

SUM2= 1 . 0 
DO 83 1=1, KL 

W=-W*( (2«*FL0ATF( I )-l o )**2)/( ( Z*8 . )*FLOATF( I ) ) 

83 SUM2= SUM2+W 

SUM2 = SUM2— W/2 » 

Rl< ( i ) = ( SUMp*SQRTF <3.141 592654/ ( 2 • *Z ) ) ) /FXPF ( Z ) 

GO TO 9 
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81 SUM 1=0.0 

DO 7 !=3i 99 . 2 

7 SUM 1 = SU M 1 +B ( I ) /FLOATF ( I — 1 ) 

SUM 1 = 4 • *SUM 1 

BK ( 1 ) =SUM 1 - ( .57721 5665+LOGF ( Z/2. )) *B < 1 ) 
9 BK < 2 ) = < < 1 • /Z ) -BK ( 1 > *B < 2 ) > /B < 1 ) 

F=BK ( 1 > 

G=BK ( 2 ) 

BK ( 3 ) = ( ( 1 o /Z ) -BK ( 2 ) *B ( 3 ) ) /B ( 2 ) 

H=BK( 3) 

END 



FUNCTION PMODF ( E . RO . Y ♦ RHO . S I GM ATR ) 

c this calculates the escape probability of neutrons incident on 

c THE MODERATOR FROM THE FUEL. IT IS CALCULATED FROM PSLOF BY A 

C RECIPROCITY RELATIONSHIP. SEE K . B • CADY THESIS. PG 88. 

C SIGMATR IS The transport CROSS section OF THE MODERATOR 

c flambda is the linear extrapolation length 

C SIGMAMA is the absorbtion CROSS section in the MODERATOR 

R 1 =RO*Y 

FLAMBDA= (.623/(1 »+2o48*SIGMATR*R0 ) )+. 71 04 
S IGMAMA=.656*Rh0*.60247*SQRTF ( • 0253/E ) / 1 8 . 02 
D=1 «/(3.*SIGMATR) 

DELTA= FLAMBDA/ SI GMATR 
FK=SQRTF (S I GMAMA/D ) 

RK1 =FK*R1 
RKO=FK*RO 

CALL BESS (RKO.FIOO.FI 1 0 , F I 20 , FKOO . FK 1 0 . FK20 ) 

CALL BESS(RKi ,FI01 ,FI 1 I .FI21 «FKOl . FK 1 1 ,FK21 ) 

B= (FK i l *F I OO+F I 1 1 *FKOO )/ (FKi 0*F I 1 I “FK \ ] *F I 10 ) 

C = B+FK*DELT A 

XE = 2./ (C*RO* ( ( Y **2 ) — 1 . )*FK) 

PMODF = 1 * - 2 « *S I GMAMA*RO* ( (Y**2 )~1 . ) *XE 
END 



C 

C 

C 

C 



FUNCTION PSLOF ( E ♦ RO , Y , RHO . S I GM ATR ) 

this calculates the escape probability of neutrons 

AT ENERGY (E) DUE TO COLLISIONS AT OTHER ENERGIES. 
THEORY is used WITH a SPATIALLY UNIFORM SOURCE and 
DEPENDENT LINEAR EXTRAPOLATION LENGTH. 

FLAMBDA = ( .623/ ( 1 c+2.48*SIGMATR*R0 ) ) + 0 7l04 
SIGMAMA=»656 *Rh0*.60247*SQRTF( .0253/E 1/18.02 
R 1 =RO*Y 

D= 1 . / (3 .*S I GMATR ) 

delta=flambda/sigmatr 

FK=SQRTF (SI GMAMA/D) 

RK 1 =FK*R1 
RKO=FK*RO 

CALL BESS (RKO.FIOO.FI 1 0 . F I 20 , FKOO , FK i 0 ♦ FK20 ) 

CALL BESS(RKi .F l0l.FIll.FI21 .FKOl . FK i l , fK21 ) 

B= (FK l i*F I OO + F I 1 1*FK00 )/ (FK i 0*FI 1 l -FK 1 1 *F I 1 0 ) 

C = B+FK*DELT A 
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PSL0F = 2»/ (C*RO* ( ( Y**2 ) - 1 . )*FI<) 
END 



FUNCTION TRAP ( S « DU « N ) 

DIMENSION S(56),DU(56) 

USES THE TRAPEZOIDAL RULE TO INTEGRATE THE N+ 1 DIMENSIONAL 
VECTOR So DU IS THE LETHARGY INCREMENTo 
C = 0.0 

DO 2 1=1 iN 

2 C=C+.5*(S( I )+S( 1+1 ) )*DU( I ) 

TRAP=C 

END 



FUNCTION SCROSS (E. T , A ) 

calculates scattering cross section/fRee atom cross section from 

THE MONATOMIC GAS MODEL. E AND T ARE IN EV. A IS IN NEUTRON 
MASSES. 

Y=E*A/T 

IF (Y-700. ) 1,1,2 

1 X= SQRTF ( Y ) 

SCROSS = ( l . + .5/Y )*FRF ( X )+l ./ ( X*EXPF ( Y ) *SQRTF (3.14159) ) 

RETURN 

2 SCROSS = ( 1 • +c5/Y ) 

END 



subroutine wateR(E,t,n,sk> 

water CALCULATES THE NELKlN WATER KERNEL. FUNCTIONS ONE AND TWO 
ARE NECESSARY. 

D I Mens I on p (P 5 ) , SK ( 55 ,55 ) , WGT ( 1 2 ) ♦ SKa (55 ) ,KN ( 5 ) , I M (5 ) » A M u (12) 

OM1 =0.205 
0M2=0 .48 1 
0MR=0 .06 
EM= 18. 

EMR = 2 ® 32 

EMV =3o/( 1 «“1 o/18.-l o /2 • 32 ) 

19 FORMAT ( 1H1 ) 

NM=N+ 1 

KN ( 1 ) = 1 

KN ( 2 ) = 1 1 

KN (3 ) =2 1 

KN (4 ) =31 

I M ( 1 )=10 

I M ( 2 ) =20 

IM(3) =30 

I M (4 ) =40 

AMU ( 1 )=. 9602899 

AMU ( 2 >=.7966665 

AMU (3) =.5255324 

AMU ( 4 ) = . 1 834346 

AMU ( 5 >=.98940093 

AMU ( 6 ) = .94457502 
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AMU ( 7 ) = « 86563 1 20 
AMU<8)=. 75540441 
AMU (9 ) =o61787624 
AMU ( 10)=. 45801678 
AMU ( 1 1 >=.28160355 
AMu ( 12)=. 09501251 
WGT ( 1 ) = . 101228=: 

WGT <2 >=.2223810 
WGT ( 3 >=.3137066 
WGT (4 > = .3626838 
WGT(5>=. 02715246 
WGT (6 ) = .06225352 
WGT(7)=. 09515851 
WGT ( 8 ) = . 12462897 
WGT ( 9 )=. 14959599 
WGT ( 1 0 )=. 16915652 
WGT (11)=. 18260342 
WGT ( 1 2 >= . 1 8945061 
DO 102 1=1 . NM 
DO 101 J=1,I 

BEE=E( I >+E( J) 

C= SQRTE ( E ( I )*E ( J) ) 

D=0. 

H= 0 . 

N J=NM+ J 

NA=.5*FL0ATF(NJ) 

IF(I-NA) 150,150.80 
150 KM=5 
KL= 1 2 
GO TO 306 
80 KMi 1 
KL = 4 

306 IF ( E < O ) “ . 33 ) 300,300.301 

301 IF(E(J)-.34> 302.302,303 
300 M J= 1 

GO TO 304 

302 M J=3 

304 0M=0MR 
AM=EMR 
L J = 3 

CM= 1 . /EM 
GX=T /EM 

B=EXPF (OM/ (2.*T ) ) 

A 1 = ( ( B**2 ) + 1 . ) / ( < < B**2 ) ~1 . >*EMR*OMR ) + i . / ( fM\/*OM } ) +2 ./ (EMv*0l v l2 > 
GO TO 305 

303 OM=OM 1 
AM=EMV 
M J = 3 

L J= 1 

CM=1 ./EM+1 «/EMR 

GX = T/EM + OMR*( i e / ( EXPF (OMR/T)— l o ) + o 5 ) /EMR 
B=EXPF(0M/(2.*T) ) 

AI = 1 o /(EMV*0Ml ) +2 » / ( EMV*OM2 ) 

305 IF(I-J) 901,902, 901 
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90? ALPHA = E < J )/SQRTF(3 « I 4 1 59*GX ) 

BETA=E( J)*SQRTF( 16«/<3* 14159*T ) ) 

GO TO 903 
901 ALPHA=Oo 
BETA=0 . 

9 O 3 20X=2o*B/(AM*OM*<B**2-l • ) ) 

DO 151 K=KM,KL 
X= BEE+2 0 *C*AMU ( K ) 

XX=SQRTF (X ) 

Z=X*ZOX 

F = ON^ (E(J),E(I ) 1 X. GX ,OM,cM, MJ,LJ, A ! , B*Zi EMv ) 

F= <2o*81 « 5/4 0 )*(F~ALPHA/XX) + (4o24/4o )* 

1 (TWO(E( J) » E ( I ) tXtT ) “BET A/XX ) 

X=BEE-2.*C*AMU (K) 

XX=SQRTF (X ) 

Z=X*ZOX 

G = 0 N£<E( 0 ) , F < I > «X.GX,OM,cM,MJ,LJ,AI , B.Z,EMV) 

G= (2. *81 .5/4. )*CG-ALPHA/XX>+ ( 4 * 24/4. )* 
1<TW0(E(J)«E(I)»X.T) -BETA/XX ) 

D=D+WGT (K )* (F+G ) 

1F1 H=H+WGT (K ) * AMU ( K ) * ( G— F ) 

IF(I-J) 70,802.70 

802 FA = ( (2 o*81 .5/4. ) * ALPHA+ ( 4 . 24/4 . ) *BET A ) *2 . /SQRTF ( E ( -> ) ) 
D= D+F A 
H= H+F A/3 0 
70 SK< J, I )=D 

101 SKA( J)=SK(J, I ) 

PRINT 1004,1 

1004 FORMAT (4HI = 13) 

1005 FORMAT (19H SCATTERING KERNEL / ( 7 ( I 3 , E 1 3 • 4 ) / ) ) 

PRINT 1 005, ( J, SKA ( J ) , J=i , I ) 

102 CONTINUE 

DO 110 J=2«NM 
JK=J— 1 

DO 110 I = 1 , JK 

BB=< <E< I >/E< J) )**2)*EXPF( ( E ( J ) ~E ( I ) ) /T ) 

110 SKI J, I )=SK( I , J)*BB 
DO 951 L=1 ,4 
N I =KN ( L ) 

Ml =1 M (L ) 

PRINT 199.T 

199 FORMAT ( 25H 1 SCATTER I NG KERNEL, KT = F7.4) 

PRINT 9199, ( I , I =N I ,MI ) 

9199 FORMAT ( 7X , 1 0 ( I 2 ♦ 9X ) ) 

951 PRINT 200, ( J , (SK( J, I ) ,J=NI ,MI ), 1 = 1 ,NM) 

200 FORMAT ( I 4 , 2X . 1 0 ( El 1 . 3 ) ) 

Nl =41 

M I =45 

PRINT 199, T 

PRINT 9199, ( I , I=NI ,MI ) 

PRINT 220« ( I . (SK( J, I ) , J=NI , M I ) , I = 1 ,NM ) 

220 FORMAT ( I 4 , 2X , 5 ( E 1 1 . 3 ) ) 

9876 CONTINUE 
END 
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FUNCTION ONe(EP«E«X,GXiOM,cM,MJ,LJ,AI ,B,Z*EMV) 

this calculates <sigma<up to u.mu>/(Sigma round atom/4>) f or 
NELKIN PROTON. EP 1 E 1 X 1 and t are in ev. 

DIMENSION B I (200 ) . AA < 3 ) .CC < 3 ) 

OM 1 =0 .205 
OM2=0 .48 1 
AA ( 1 ) = 1 « 

AA ( 2 ) = X/ (EMV*0M1 ) 

AA ( 3 ) = ( AA (2 )**2 )/2» 

CC ( 1 > = 1 c 

CC(2) =2.*X/ (EMV*0M2 ) 

CC(3)=(CC(2)**2)/2. 

PROD 1 =0.0 
SUM=0 . 

D=E*SQRTF(E/(EP*X*GX*3. 14159 ) ) 

N= 1 . 5*Z+ 1 0 o 
BI (N) =0. 

SUM= 1 oOE-20 
B I (N- 1 ) =SUM 
NN=N-1 

DO 11 J=2.NN 
I =N-J 

BI ( I ) =2.*FL0ATF< I )*BI ( 1 + 1 >/Z+BI (1+2) 

1 1 SUM=SUM+BI ( I ) 

SUM=2**SUM-BI ( 1 ) 

C=EXPF (Z )/SUM 

DO 12 J=1 » NN 

12 BI ( J)=C*B! ( J) 

DO 1050 M=l ,MJ 
FM = M— 1 

PR0D = 0 « 

DO 17 L= 1 iLJ 
EL=L-1 
SUM=0 . 

DO 13 N= 1 ,NN 
EN=N— 1 

BB=E-EP+X*CM+EN*0M+EL*0M1 +FM*0M2 
AZ= ( BB**2 ) / ( 4 . *X*GX ) 

I F ( AZ — 705 . ) 112.112.113 

112 P=( (B**EN)*BI (N) )/EXPF(AZ) 

GO TO 114 

113 P= 0 . 0 

114 SUM= SUM + P 
W= P/SUM 

I F ( W— 1 0 OE — 4 ) 14.14.13 

13 CONTINUE 
PRINT 500, N 

500 FORMAT (5H N = 13) 

14 DO 15 N=2 . NN 
EN=N- 1 

BB=E-EP+X*CM-EN*0M+EL*0M1+FM*0M2 

AZ=(BB**2)/(4.*X*GX) 

I F ( AZ — 7 O 5 . ) 115.115.116 

1 1-5' P= (B I ( N)/(B**EN ) )/EXPF ( AZ ) 
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GO TO 117 
1 16 P=0®0 
1 1 7 SUM= SUM+P 
W=P/SUM 

I F ( W— 1 • OE — 4 ) 16,16,15 

15 CONTINUE 
PRINT 500, N 

16 PROD=AA (L )*SUM+PROD 

17 CONTINUE 

PROD 1 = CC (M )*PROD + PRODl 
1050 CONTINUE 

ONE = D*PRODl /EXPF ( X*A I ) 
END 



FUNCTION TWO ( EP , E , X , T ) 

this calculates (sigma(up to u,mu>/(sigma bound atom/4) > 

FOR A FREE GAS OF MASS 16* EP,E,X,AND T ARE IN EV* 
B=E-EP+X/16» 

D=E*SQRTF(E*16./<X*T*3. 1415P*EP) ) 

TWO = D/EXPF( (B**2 )*1 6®/ (4**X*T ) 1 
END 

END COUTH 



SUBROUT I NE FKERN ( NM , SFA , T, A, E, SK ) 

subroutine fkern calculates the scattering matrix using THE 

MONATOMIC GAS MODEL. 

D I MENS I ON E <55 ) *SK (55 ,55 ) , W (55 ) ,RX (55 ) , EX (55 ) 

ET A = < A+l * )/ ( 2**SQRTF ( A ) ) 

R 10 = < A-l * )/ ( 2«*SQRTF ( A ) ) 

DO 43 1=1, NM 

W ( 1 ) = SQRTF ( E ( I )/T ) 

RX ( I ) =R I 0*W ( I ) 

43 EX ( I ) =ETA*W ( I ) 

DO 10 1=1, NM 

DO 1 0 J= 1 , I 

1 0 SK ( J, I ) = ,5*SFA* (FTA**2)*(F( I )/F(J) )*( fRf ( FX ( I )“Rx < J ) )+E R F ( EX ( I ) + 
1 RX ( U ) ) + ( EXPF ( ( W ( J ) **2 )“ ( W ( I ) ) ) ) * (E R F (EX ( J > ~Rx ( I ) )-E R F ( R X ( I ) + 
2EX ( J ) ) ) ) 

END 
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NOTATION 

Roman Letters 

f 

a Fuel rod radius in mean free path, a = S r^ 

A Fuel element area, Eq . (2.12) 

A Ratio of scatterer mass to neutron mass, 

Eq. (5.1) 

A Parameter for substitution into the Nelkin 

Water Kernel, Eqs . (5.17) and (5.18) 

A Mass number, Eq . (7.3) 

c Ratio of the scattering to the total cross 

section, p. 27 

C(E) Correction factors to the diagonal terms in 

the Nelkin Water Matrix, p. 59 

C-,^ Eq. (4.36) 

C 2 Eq. (4.37) 

d(E) Linear extrapolation length into the fuel 

element, Eq . (2.18) 

d Q (E) Linear extrapolation length into a black hole, 

Eq. (C .4) 

D Diffusion coefficient for the homogenized cell, 

Eq. (6.38) 

D x Eq. (4.30) 

D 2 Eq. (4.32) 
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D (E) 
m 


Diffusion coefficient for the moderator, Eq . (2.13) 


D f (E) 


Diffusion coefficient for the fuel, p. 27 


E 


Neutron energy in eV 


E c 


Thermal cutoff energy 


f 


Thermal utilization, Eqs . (6.13) and (6.15) 


f (M-) 


Eq. (5.21) 


f(E) 


Eqs. (2.7), (2.28), and (3.8) 


F o 


Eq. (4.26) 


F 1 


Eq. (4.27) 


F 2 


Eq. (4.28) 


F 3 


Eq. (4.29) 


8 


Sec. 4.3, Eq. (4.45) 


G(E) 


Eq. (2.17) 


h 


Planck's Constant 


V z > 


Modified Bessel function of the first kind 


J_ (E) 


Current into the fuel element, Eq . (2.24) 


V E > 


Net current entering the cladding from the 
fuel rod, Eq . (4.22) 


J n (E) 


Net current entering the fuel element, 
Eq. (2.12) 


K n< z > 


Modified Bessel function of the second kind 




Thermal lifetime, Eq . (6.17) 


L 


Diffusion length for the homogenized cell, 
Eq. (6.39) 




Eq. (4.33) 
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"1 

3t\ 

Xl 

M 

M(E) 



N(E) 

N(r,E) 

N^E) 

N 1 

P e (E) 

P m< E > 

P r (E) 

P rv< E > 

p s (E) 

P 

u 

P unif 

r o 

r l 

r 

c 

r f 



Eq. (4.34) 

Eq. (4.38) 

Eq. (4.39) 
p. 47 

Maxwellian distribution, Eq . (2.6) 
Cylindrical moments. Sec. 4.1 
p. 47 
p . 47 

Spatially averaged source, Eq . (2.11) 

Neutron source, Eq . (2.11) 

Spatially averaged energy dependent neutron 
density in region i, Eq . (6.5) 

Total neutron density in region i, Eq . (6.7) 
Fuel element escape probability. Chap. 4 
Moderator escape probability. Sec. 3.2 
Fuel rod escape probability. Sec. 4.1 
Escape probability for the fuel rod and void 
region, Sec. 4.2 

Moderator escape probability. Sec. 3.1 
Eq. (4.43) 

Eq. (4.44) 

Fuel element radius 
Equivalent cell radius 
Inner radius of the cladding 
Fuel rod radius 
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r . 
in 


Sec. 4.3, Eq. (4.49) 


r* 

out 

RRj 


Sec. 4.3, Eq. (4.50) 

Reaction rate of type j in region i, Eqs . (6.18) 
and (6.19) 


S(E) 


Neutron slowing down source, Eq . (2.3) 


S(E) 


Spatially averaged slowing down source of 
neutrons, Eqs. (2.10), (5.27), and (7.1) 


S(r,E) 


Spatially dependent slowing down source of 
neutrons, Eq . (2.10) 


t 


Eq . (4.24) 


T 


Moderator temperature in electron volts 


T 


Sec. 4.3, Eq . (4.35) 


u 


Lethargy 


v(E) 


Neutron velocity 


v 1 


Average neutron velocity in region i, Eq . (6.6) 
The volume of region i 


w k 


Weighting coefficients for the Gaussian 
quadrature integration formula, p. 50 


X 


Eq . (4.42) 


X 


Eq. (5.2) 


y 


p. 15, y = r 1 /r Q 


y 


Eq. (5.14) 


z 


Eqs. (5.17) and (5.18) 


z 


Eq. (5.5) 
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Greek Letters 



a(E) Fuel capture to fission ratio 

P(E) Eq. (4.2) 

7 Euler's Constant 

6™ Moderator neutron density disadvantage factor, 

Eq. (6.8) 

6^ Cladding neutron density disadvantage factor, 

Eq. (6.9) 

6™ Moderator flux disadvantage factor, Eq . (6.10) 

6^ Cladding flux disadvantage factor, Eq . (6.11) 

e Convergence criterion, Eq . (C.24) 

C Eq. (4.20) 

7) Number of fission neutrons produced per thermal 

absorption in the fuel, Eq . (6.16) 
r) Eq. (5.7) 

0 Eq. (4.25) 

K p . 15 

/Cf p. 28 

A Eq . ( C . 4 ) 

TV Eq. (4.31) 

p Cosine of the scattering angle 

p(E) Average cosine of the scattering angle, Eq . (5.24) 

p^ A zero of the Legendre polynomials, p. 50 

v Number of fission neutrons produced per thermal 

fission 

1 Eq. (4.23) 

£ Eqs . (5.17) and (5.18) 
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P 

Pin 

Pout 

a b 

a f 

°J< E > 

a(E 

a(E 

a*(E- 

a 1 (E' 

^(E) 

Z. 

J 

I 1 : 

J 

Z(E ' — 

T 

<p m (E) 



= 3.14159 

Eq. (5.8) 

Sec. 4.3, Eq . (4.51) 

Sec. 4.3, Eq. (4.52) 

Bound atom microscopic cross section 
Free atom microscopic cross section 
Microscopic cross section expressed in barns 
for interaction type j in material i 
E,|_i) Differential scattering kernel, Eqs . (5.1) 
and (5.15) 

E) Scattering kernel, Eqs. (5.4) and (5.19) 

E) Modified diagonals for the Nelkin water 

kernel, Eq . (5.26) 

^E) First moment of the scattering kernel, 

Eq. (5.23) 

Macroscopic cross section expressed in cm ^ 
for interaction type j in material i 
Macroscopic effective cross section for inter- 
action type j in material i, Eqs. (6.23) and 
(6.31) 

Macroscopic mean cross section for interaction 
type j in material i, Eqs. (6.24) and (6.32) 

E) Macroscopic scattering kernel, Eq. (2.1) 

Sec. 4.3, Eq . (4.40) 

Energy dependent flux in an infinite homogeneous 



moderator, Eq . (2.1) 
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cp x (E) Spatially averaged energy dependent flux in 

region i 

cpg Spatially averaged 2200 m/sec flux in region i 

cp^ Spatially averaged thermal flux in region i 

^(u.) k L iteration for the moderator flux at 
lethargy u^, Eq. (C.22) 
cp£(u.) cp^(u^) normalized, Eq . (C.23) 

9 u (r) Eq. (4.5) 

X Eqs . (5.17) and (5.18) 

X (E) Eq. (2.21) 

<x> r p. 47 

od^ p. 47 

o>2 p. 47 



Subscripts and Superscripts 
a absorption 

c cladding 

f fuel 

f fission 

h homogenized cell 

m moderator 

s scattering 

t total 



tr 



transport 
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